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El cacique de Iraca y su sobrino Ramiquin 
gobernaban sobre la tierra en una noche absoluta. 
Para resolver la situacion, el cacique de Iraca 
decidio que su sobrino ascendiera sobre los cielos 
y trajera la luz. Este se dirigio vertiginosamente 
hacia las alturas y de pronto se transformo en 
un astro incandescente y luminoso: Ramiquin' 
se habfa convertido en el Sol. Pero su tib 
no estaba satisfecho del todo pues una parte 
del dfa se hallaba aun en tinieblas y esto le 
recordaba a la humanidad, con miedo y tris- 
teza, la epoca en que todo era tinieblas. Fue 
entonces cuando el cacique de Iraca resolvio 
hacer lo mismo que su sobrino, perdiendose 
en la boveda celestial. Y se convirtio en un 
astro de luz mas tenue: la Luna. Su luz servia 
para alegrar a la gente durante la ausencia del Sol. 

Leyenda muisca o chibcha. Colombia. 



Resumen 



Este trabajo de tesis se enmarca en un concepto mas general denominado "Alta 
resolucion en ffsica solar" . El trabajo consiste de dos partes claramente definidas. 
La primera parte trata sobre el desarrollo instrumental para observaciones solares y 
la segunda parte esta dedicada a la explotacion cientifica de datos solares obtenidos 
con intrumentacion solar puntera. 

En la primera parte de la tesis se trabaja el tema de la alta resolucion y la restau- 
racion de imagenes para la obtencion de una alta calidad de imagen. Se comienza 
con una revision teorica del problema que representa la turbulencia atmosferica y las 
aberraciones instrumentales en la calidad de las imagenes. Esto plantea la necesi- 
dad de implantar tecnicas de restauracion post-facto de las imagenes, que sumadas 
a las correcciones en tiempo real de la Optica Adaptativa, nos den imagenes cada 
vez mas cercanas a la realidad, es decir al objeto verdadero, que en nuestro caso 
es la region del Sol que queremos estudiar. 

La forma evidente, aunque no por ello la mas sencilla, de evitar el efecto negativo 
de la turbulencia atmosferica sobre la calidad de nuestras imagenes es el uso de 
telescopios espaciales. Fuera de la atmosfera terrestre, las observaciones no estan'an 
afectadas por aberraciones atmosfericas. Sin embargo, aun seguinan existiendo 
aberraciones instrumentales degradando las imagenes, aunque con menos intensi- 
dad. El problema principal de una mision espacial es su elevado costo de puesta en 
orbita, mantenimiento y actualizacion. 

Un esfuerzo por tener una observacion solar, sin el efecto contra prod ucente de 
la atmosfera terrestre, es el desarrollo de la mision SUNRISE, una colaboracion 
entre la Agencia Espacial Alemana, DLR, la Estadounidense NASA y el Programa 
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Nacional Espanol del Espacio. Este proyecto lanzara un globo aerostatico con un 
telescopio de 1 metro de apertura que tendra durante 15 di'as la posibilidad de ob- 
servar ininterrumpidamente al Sol con alta resolucion espacial, temporal y espectral 
sin precedentes. El objetivo principal de SUNRISE es el estudio de la formacion de 
estructuras magneticas en la atmosfera solar y su interaccion con los flujos convec- 
tivos de plasma. Para cumplir este objetivo se cuenta con el instrumento Imaging 
Magnetograph experiment (IMaX), un magnetografo desarrollado enteramente por 
instituciones espanolas, lideradas por el Institute de Astrofi'sica de Canarias en el 
cual he realizado el presente trabajo. Este instrumento sera capaz de producir ma- 
pas del campo magnetico de regiones extensas de la superficie solar midiendo la 
polarizacion de la luz en determinadas Imeas espectrales. Como miembro del equipo 
de IMaX, he desarrollado un metodo de calibracion en vuelo para caracterizar las 
aberraciones que afectaran las imagenes en IMaX. La descripcion del metodo de 
calibracion, asi" como tambien las pruebas de su robustez, constituyen el nucleo de 
la primera parte de esta tesis doctoral. 

En la segunda parte de la tesis nos centramos en el tema del estudio de flujos 
horizontales en regiones solares activas. Se utilizan dates de observaciones so- 
lares desde Tierra y desde el espacio y se aplica el metodo de reconstruccion de 
imagenes expuesto en la primera parte para restaurar el material observado. Estu- 
diamos los movimientos propios de estructuras dentro y fuera de regiones solares 
activas. A traves de tecnicas de correlacion local y los subsiguientes mapas de flujo 
que generamos, podemos cuantificar los flujos horizontales en las regiones obser- 
vadas. 

La primera region activa estudiada corresponde a un complejo grupo de man- 
chas solares de configuracion 5. Se infiere el campo de velocidades horizontales 
sobre de una serie temporal de alta resolucion, y a partir de este mapa de flujos se 
encuentra una correlacion entre la presencia de flujos de gran velocidad alrededor 
de las manchas solares hacia afuera y la existencia de penumbra. La zona afectada 
por estos flujos se denomina foso (en ingles moat). Se sugiere una relacion entre 
flujos radiales hacia afuera a lo largo de los filamentos penumbrales (flujo Evershed) 
y los flujos fotosfericos, tambien radiales, en la granulacion circundante a las man- 
chas solares. Para confirmar este resultado, se estudia una muestra mas amplia de 
manchas solares con gran variedad de configuraciones penumbrales, y nuevamente 
se encuentra la misma dependencia foso-penumbra. En las areas donde las umbras 
son adjacentes a la granulacion circundante, no hay evidencias de presencia de estos 
flujos a gran escala. 
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Finalmente se estudia el campo de velocidades horizontales alrededor de poros 
(manchas solares sin penumbra). Se trabaja con imagenes restauradas de alta 
resolucion, que conforman una serie temporal estable y de larga duracion, sobre la 
cual se estudian nuevamente los movimientos propios en todo el campo de obser- 
vacion y sus propiedades alrededor de poros solares. Como resultado relevante, no 
se encuentra ninguna evidencia de flujo de foso alrededor de los poros. 



Abstract 



The present thesis work can be framed in a more general concept designated as 
"High resolution in solar physics". The work consists of two clearly defined 
parts. The first part concerning instrumental development for solar observations 
and the second one devoted to the scientific exploitation of solar data acquired with 
cutting edge solar instrumentation. 

The first part of this thesis is dedicated to the topic of high-resolution ob- 
servations and image restoration to obtain high-quality images. It begins with a 
theoretical reviewing of the problem that represents the atmospheric turbulence 
and the instrumental aberrations on the image quality. This problem force us to 
implement post-facto image restoration techniques that, added to the real-time 
corrections performed by the Adaptive Optics, gives us images closer to reality, i.e. 
to the true object, that in our case is the region of the Sun we want to study. 

The most evident solution, although not the easiest one, to get rid of the neg- 
ative effect of the atmospheric turbulence on the image quality, is the use of space 
telescopes. Out from the Earth's atmosphere, the observations would not be af- 
fected by atmospheric aberrations. Nevertheless, the instrumental errors would still 
be present degrading the images although with less strength. The counter part 
of this solution is the elevated costs of launching, maintenance and updating of a 
space instrument. 

To have good solar observations overcoming the negative influence of the Earth' 
s atmosphere, one effort is being made with the development of the so-called SUN- 
RISE mission, a collaboration between the German Space Agency DLR, the north 
American NASA and the National Spanish Space Program. This project consists 
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in a balloon-borne mission that will launch a 1-meter telescope to the stratosphere 
and will record data uninterruptedly during 15 days, with unprecedented temporal, 
spatial and spectral resolution. The main aim of SUNRISE is to study the for- 
mation of magnetic structures in the solar atmosphere and their interaction with 
the convective plasma flows. In order to do so, the on-board instrument called 
Imaging Magnetograph experiment (IMaX) has been designed and implemented 
entirely by Spanish institutions leaded by the Instituto de Astrofisica de Canarias in 
which I have developed the present work. This instrument will be able to produce 
magnetic field maps of extensive solar regions by measuring the light polarization 
in certain spectral lines. As a member of the IMaX team, I have been in charge 
of performing numerical simulations to identify and evaluate possible optical error 
sources. My main contribution to the project has been the development of an 
in-flight calibration method to characterize the aberrations affecting the images in 
IMaX. The description of this calibration method as well as the test to prove its 
robustness make up the core of the first part of this thesis. 

The second part of the thesis is centered on the study of horizontal flows in 
solar active regions. Data from ground-based and space observations are used and 
the reconstruction techniques explained in the first part are successfully employed 
to restore the images. We focus on the proper motions of structures in and around 
solar active regions. The way to quantify the horizontal flows in the field-of-view 
consist of using local correlation tracking techniques that generate flow maps. 

The first active region studied corresponds to a very complex sunspot group 
with J-configuration. The horizontal velocity field is inferred from a high-quality 
time series and, from these flow maps, it is found a correlation between the presence 
of strong velocity flows (moats) surrounding sunspots and the existence of penum- 
bra. A relation between radial outflows along the penumbral filaments (Evershed 
flow) and the photospheric outflows in the granulation surrounding the sunspot is 
then suggested. To confirm this result, a larger sample of sunspots with a variety 
of penumbral configurations is studied and once again the same moat-penumbra 
dependence is found. In the areas where umbrae are adjacent to the surrounding 
photospheric granulation there are not evidence of these large-scale outflows. 

Finally, the horizontal velocity field is studied around pores (sunspots lacking 
penumbrae). Working with restored high-resolution images conforming a stable and 
long-duration time series, the analysis of proper motions is performed again within 
the field-of-view and the properties of the motions around solar pores are studied 
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and quantified in detail. As a relevant result, there are no evidences of moat-like 
flows around the pores. 



Introduction 



The Sun is undoubtedly the most important astronomical object for the human 
kind. Although its proximity and the fact that it has been largely studied, our 
knowledge is yet quite poor in understanding how it really works: which are the 
mechanisms that rule its behaviour, with the eleven years solar cycle, and further 
more facts concerning the magnetic nature and its implications in the structure of 
our closest star. 

The study of the Sun embraces nowadays multiple branches: from the study of 
the solar interior through helioseismology techniques up to the study of the solar 
wind that impacts the Earth after traveling a long distance across the emptied 
space. The majority of advances in solar physics have come thanks to the new 
telescopes and observational techniques. The quality of the images taken by tele- 
scopes has improved considerably in the last decades and nowadays it has almost 
reached the diffraction limit imposed by the instrument. Computational tools have 
also taken part in the development of solar physics by improving modelling and 
simulations which are in their turn corroborated by observations, hereby the impor- 
tance of reaching high-resolution observing levels. 

The Sun evolves and changes in different time and spatial scales. Added to the 
11 years cycle that marks the biggest changes in the solar magnetic field structure, 
it also presents more rapid changes that can be seen from Earth. Appearance of 
sunspots, for instance, was first observed in the XVII century with the use of the 
telescope in astronomy though there were many records of sunspots observations 
since antiquity. 

The aim of this thesis associated to the Imaging Magnetograph experiment 
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(IMaX) project developed by several Spanish institutions 1 has two main branches. 
Firstly, the initial part of the thesis deals with instrumental aspects, reviewing the 
latest restoration techniques and concerning the design of a method to calibrate 
the in-flight instrumental aberrations of IMaX. The second part of the thesis is 
dedicated to the study of the dynamics of structures in the solar photosphere, in 
particular in solar active regions. A more detailed information of the contents for 
each chapter is described below: 

In Chapter §1 we start reviewing the theory concerning the image formation 
which is crucial to understand many concepts developed through all the thesis, 
such as the degradation of images and the restoration techniques applied to recon- 
struct them. 

Chapter §2 is devoted to the instrument IMaX, describing first the project it 
is associated to, the objectives of the mission and some important characteristics 
about it. Moreover, the main concern of this section is the definition of the method 
for the IMaX in-flight calibration and the subsequent test of the robustness of the 
method taking into account all possible error sources. 

Chapter §3 is the first chapter within the second part of the thesis. It intro- 
duces the theory on solar active regions and the current knowledge in the field, i.e. 
the description of the structures conforming a solar active region, ranging from the 
widely studied umbra and penumbra of sunspots to the more recently discovered 
fine structures. 

In Chapter §4 we employ the restoration techniques to correct solar data and 
show the first results yield by the study of horizontal proper motions in and around 
sunspots. Maps of horizontal velocities are presented for a complex solar active 
region with a ^-configuration. Some of the results in this chapter have already 
been published in the The Astrophysical Journal Letters (ApJL). 

Chapter §5 extends the study of proper motions to a larger sample of sunspots 
displaying various penumbral configurations. Most of the results presented in this 
chapter have already been published in the The Astrophysical Journal (ApJ). 

Chapter §6 extends the study of solar active regions done in the previous two 

1 lnstituto de Astroffsica de Canarias (IAC), Instituto de Astroffsica de Andalucfa (IAA), Instituto 
Nacional de Tecnica Aeroespacial (INTA), Grupo Astronomfa y Ciencias del Espacio (GACE) 
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chapters by analyzing a sample of solar pores from ground-based and space obser- 
vations. This chapter is a extended version of the work to be published in the The 
Astrophysical Journal (ApJ). 



Chapter §7 summarizes the conclusions of the work and the final discussion. 
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Part I 

Defining a method for in-flight 
calibration of IMaX aberrations 
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1 

Foundations on image restoration 

n the first part of this thesis I concentrate on the techniques we employ to achieve 
high-resolution solar images that enable us to study the Sun at very small spatial 
scales (fine details). 

1.1 Introduction 

Current problems trying to explain the physics of the Sun require to resolve very 
tiny structures as the first step to be able to model what is actually happening in 
different layers and regions in the Sun. 

The Earth's atmosphere can be considered as an isotropic turbulent medium. 
The quality of the solar images captured by ground-based telescopes is severely af- 
fected by the atmospheric turbulence. The image degradation is generally described 
as the combination of three main contributions, as follows: 

• Structures smearing (blurring). 

• Global displacements of the image (image motion). 

• Distortion of the structures caused by the differential image motion of differ- 
ent patches in the FOV (stretching). 

All these degradation effects considered together are usually referred to with 
the term seeing and represent the first problem to face if we are interested in high 
resolution data. 
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There have been recent efforts trying to circumvent the atmospheric influence 
on the observed images. Space-based telescopes (Domingo, Fleck & Poland, 1995, 
SOHO) are undoubtedly the obvious best choice to get rid of the atmospheric ef- 
fects. We all have got fascinated with the images obtained by space facilities and 
nowadays we are still getting astonishing data from them (HINODE; Kosugi et al., 
2007). But the elevated cost of construction and operation of these space solar 
telescopes has made the scientific community to think about an alternative. 

In recent years and in order to improve the ground-based observations, the 
Adaptive Optics (AO Rimmele, 2000; Scharmer et al., 2000) has made possible to 
partially correct both, the instrumental and the atmospheric aberrations. The idea 
is conceptually simple and uses optical elements deforming in real-time to compen- 
sate the wavefront aberration induced by the atmosphere and the telescope. 

Due to the temporal scale in the evolution of the seeing, the AO only pur- 
sues low-order corrections that moreover are limited to an isoplanatic patch of a 
few arc seconds, so that it is compulsory to implement post-facto computational 
techniques to complement the real-time corrections. Powerful numerical codes for 
image restoration have been thus developed in the last decade, every single one 
requiring an especially designed observing strategy. 

In order to quantify how the aberrations affect the images, it is worth to dedicate 
the next section to the mathematical formalism describing the image formation in 
the telescope. 

1.2 Image formation 

The basic process describing the formation of images is shown in Figure 1.1a. The 
object of study is placed on the so-called object plane represented by the coordi- 
nates system (£,??), and another coordinate system (x, y) called the image plane is 
situated right in the focal plane of the telescope. The wave coming from a point 
source in the object plane passes through the atmosphere and part of it enters the 
telescope and forms the corresponding image (not a point anymore but a spot) 
centered at position (x',y') in the image plane. 

If the object source is the impulse unit function its image corresponds to the 
impulse response function of the optical transmission system (atmosphere + tele- 
scope). This response is also named Point Spread Function (PSF) of the system. 
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FIGURE 1.1 — a) Skecth representing the image formation process, and b) Detail of the image 
plane showing that the intensity at each point has a contribution from the distribution of intensities 
centered at points in the neighbourhood. 



The PSF can be viewed as the normalized distribution of intensity in the image 
of a point source and can be expressed as s(x,y;x',y';t), where t represent the 
variability of the transmission system (e.g. atmospheric turbulence evolves at a 
fast rate) and (x',y') reflects that the PSF is, in general, space variant, i.e. the 
distribution of intensity in the image of a point source changes with its location 
(£, rj) in the object plane. Note that (a/, y') is the image conjugated point of (£, rf). 

Considering an extended object like the Sun as composed by multitude of in- 
coherent point sources, with individual intensities i (^,7?), we can estimate the 
resulting intensity distribution in the focal plane of the telescope. Assuming a 
linear optical system and incoherent illumination, this intensity can be expressed as 

i(x,y,t)= // i (x',y')s(x,y;x',y';t)dx' dy', (1.1) 



where i (x',y') represents the distribution of intensity in the ideal image (i.e. 
also in the object) that would produce a perfect system free from aberrations 
and with infinite aperture. The intensity at each point (x,y) of the image has 
a contribution from the images centered at points in the neighbourhood, e.g. 

(x'i,y'i),(x' 2 ,y' 2 ),(x' 3 ,y' 3 ) in Figure 1.1b . 
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For an isoplanatic system having a spatially invariant PSF, i.e. 

s(x, y; x', y'\ t) = s(x — x',y — y', t), the equation 1.1 can be written as 




(1.2) 



where * stands for convolution and q is the vectorial notation for the coordinates in 
the image points. Using the convolution theorem 1 , the intensity can be expressed 
in the Fourier domain as the product 



where capital letters stand for the Fourier transforms of the functions in lowercase, 
u is the frequency vector in the Fourier domain and, S is the so-called Optical 
Transfer Function (OTF) of the system. The amplitude of the OTF is defined as 
the Modulation Transfer Function (MTF): MTF(u, t) = |OTF(u,i)|. According 
to equation 1.3 the MTF is a filter that attenuates the amplitudes of the Fourier 
components of the ideal image to form the observed image. 

1.3 Image restoration as a particular case of the Inverse Problem 
in Physics 

Image restoration fits into the format of the Inverse Problem in Physics 2 which, in 
general, can be considered to as the solution f(q') of the Fredholm inhomogeneous 
integral equation of the I st kind 



where JC(q,q') is known as the kernel of the integral equation. Equation 1.2 is a 
particular case of equation 1.4, where the kernel is the PSF. 

The inverse problem in our particular case is usually called image reconstruction 
or image restoration since the target is to achieve an estimate of the ideal image 

1 The convolution theorem states that the Fourier transform of the convolution of two functions 
is the product of their respective Fourier transforms. 

2 The inverse problem frequently occurs in different branches of science where the values of 
some model parameters must be obtained from the observed data. 



I(u,t) =I (u)S(u,t) 



(1.3) 




(1.4) 
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or equivalently of the true object starting from a degraded image. Because in our 
particular case the inverse problem can be formulated as a convolution equation 
(eq. 1.2), image restoration can also be referred to as a deconvolution problem, the 
formal solution of which can be expressed from equation 1.3 as follows, 



Io(u) 



S(u,t)' 



or in the measuring domain as 



i (q) = T 



-i 



I(u,t) 
S(u,t) 



(1.5) 



(1.6) 



The symbol * over iq means that the restoration will not be complete in any 
case since the transmission system operates as a low-pass spatial frequency filter 
with a given cut-off. In other words, the restoration problem will render an estimate 
of the true object. 



In order to compute i {q) it is mandatory to characterize the PSF (or its Fourier 
transform the OTF) describing the optical system. Several approaches to the PSF 
determination lead to different numerical methods for the restoration of solar im- 
ages (Bonet, 1999). 



1.3.1 Noise contribution 

An additional difficulty in the inversion problem arises from the fact that in the real 
case, the observed image is affected by noise caused by different sources being the 
readout and the photon noise the more relevant components. Though the latter is 
proportional to the square root of the number of photons in the incoming signal, 
in most solar physics applications the assumption of uncorrelated signal and noise 
gives good results. In this manner, equation 1.2 can be completed including the 
noise as and additive contribution, and reformulated as 3 

i n (q) =io(q)*s(q)+n(q), (1.7) 

and in the Fourier transformed domain 

Note that the variable t is dropped from the formulae describing the image formation in order 
to shorten the notation. However, one has to keep in mind that the formulae describe instantaneous 
events. 
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I N {u) = I {u) S{u) + N{u) = I(u) + N(u). (1.8) 
The restoration based in formula 1.5 leads to 

In(u) _ I(u) N(u) 

S{u) S{u) + S(u) ' 1 ' 

where the term N(u)/S(u) represents a noise amplification, making compulsory 
to pursue a noise filtering, previous to the restoration process. 

1.3.2 Noise filtering 

Filtering of the noisy signal is standardly done by using the so-called optimum 
filter ®n(u), described for ID problems by Brault & White (1971), which is a 
real function that weights the diverse Fourier spectral components according to the 
noise level at each frequency. 

I F (u) = I N (u).<S> N (u) = [I(u) + N(u)] $ N (u). (1.10) 
The optimum filter is formulated as 

, s |J(tx)| 2 

W = MW' (L11) 

Combining the noise filtering together with the deconvolution, we eventually 
find the optimum filter for restoration <£^(«), better known as the Wiener-Helstrom 
filter 4 



R{ } |i»| 2 + |iV(«)| 2 5(«) MTF 2 ( W ) + l/SNR(w)' { ' 1 

where SNR(w) = |/ (it)| 2 /|A/"(u)| 2 is the signal-to-noise ratio (SNR). Since we 
do not know a priori the function \I (u)\ 2 , some models for SNR(u) are com- 
monly assumed (Collados, 1986). Hereafter the superscript * stands for complex 
conjugate. 



4 Named after the optimal estimation theory of Norman Wiener, this filter simply acts separating 
signals based on their frequency spectra. The gain of the filter at each frequency is determined by 
the OTF of the system and the relative amount of signal and noise at that frequency. 
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FIGURE 1.2 — Sketch of a possible setup for the implementation of the phase diversity technique. 

1.4 The Phase Diversity Technique for image restoration 

The Phase Diversity (PD) technique was first proposed by Gonsalves & Childlaw 
(1979) as a new method to infer phase aberrations working with images of extended 
incoherent objects formed through an optical system. This technique of image 
reconstruction (see also Gonsalves, 1982; Paxman et a I . , 1992) requires to use at 
least two images of the object we want to reconstruct. One of the images is the 
conventional one degraded by an unknown aberration (atmosphere and telescope) 
and the second one is a strictly simultaneous image of the same object affected 
by the same unknown aberration plus a known intentionally induced aberration. 
Figure 1.2 shows an optical setup where we induce a known amount of defocus by 
displacing backwards the camera with respect to the nominal focus. 

The system of equations that mathematically describes the formation of both 
images can be expressed as 




(1.13) 



where : 



*l (?) 1*2(9) ' s the focus-defocus image pair. 

•si(q), £2(9) are the corresponding PSF's for both channels. 

n i{q), ^2(9) a re the noise additive terms. 

i {q) is the so-called true object. 



If we do not consider the noise terms the system above is determined since 
it has two equations and two unknowns: i (q) and s\(q). Note that s\(q) and 
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S2(q) differ in a well-known defocus value and consequently are analytically related. 

Nevertheless, the noise terms force a statistical solution of the problem. Pax- 
man et al. (1992) propose a least squares solution for the system of equations. For 
Gaussian noise, they derived the following error metric in the Fourier domain to be 
minimized 

L(I , Si) = ]T{|/iH - / (w)5i(w)| 2 + 7 \l 2 (u) - I (u)S2(u)\ 2 }, (1.14) 
u 

where capital letters stand for the Fourier transforms of the functions in lowercase 
in the system of equations 1.13, 7 = crf/a^ is introduced by Lofdahl & Scharmer 
(1994) as a weighting factor to equalize the noise contributions for the case when 
the noise variances a\ and a\ are not the same in both images. S represents the 
OTF of the system which can be derived as the auto-correlation of the so-called 
generalized pupil function H, 

S(u) oc J J H(r).H*(r - Xfu)dr, (1.15) 

where A and / stand for the working wavelength and the effective focal length 
of the optical system, respectively, and u is a vector with dimensions of spatial 
frequency. H, in turn, can be expressed in terms of the joint phase aberration <fi(r) 
caused by the telescope and the turbulence as 

H{r) = |lT(r)|e^. (1.16) 

The phase aberration can be parametrized by using a Zernike polynomials ex- 
pansion (see Appendix A) 

4(r) = <f>(Rp) = 1 £ i a j Z j (p), (1.17) 

where R is the pupil radius, p is the vectorial notation for (p, 9) with p defined in 
a unit circle (0 < p < 1) and 9 being the azimuthal angle. 

Thus, equation 1.16 can be then written as 

H(Rp) = \H{Rp)\A S{R P^U a ^P)], (1.18) 

where a = {ctj,j = 1,2, J} are the coefficients of the terms in the Zernike 
expansion and S(Rp) is the diverse phase, i.e. the induced defocus in our case. 
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The final result is that the OTF (equation 1.15) and therefore the error metric L 
(equation 1.14) are parametrized by the expansion in Zernike polynomials, so that 
one can write: L(I ,a). The diverse phase S(Rp) will be zero in the case of Si 
that corresponds to the focus image. 

Part of the minimization of the equation 1.14 can be performed analytically. 
The solution of the equation dL/dI Q = yields an object estimate I (u) that 
minimizes the equation 1.14, 

f (,,\ _ J i(") g iK a )+7^(u)5 2 *(u,a) 

IoW ~ |5i(« ) a)|2 + 7 |5 2 (« ) a)|2 ' ^> 

Replacing I (u) in the equation 1.14, we obtain a so-called modified error 
metric 

T i \ ^ |/i(u)S 2 (u,q) - ; 2 (ii)5i(u,Q)| 2 

Lm(ol) = > 1 r——, 7775 rm ttt — ■ (l-^ u ) 

V \Si(u,a)\ 2 + 1 \S 2 (u,cx)\ 2 

Note that this modified metric is not explicitly depending on the Fourier trans- 
form of the object I (u) but only includes the unknown a vector. By means of 
non-linear optimization techniques, we find the a vector, characterizing the aber- 
ration components, that minimizes the equation 1.20. Once these components are 
determined, Si and S 2 can be calculated from equations 1.18 and 1.15 and the 
object estimate can be eventually derived from equation 1.19, completing in this 
way the restoration process. 

In conditions of poor seeing, the restored images are sometimes contaminated by 
some artifacts like periodic strips or other regular patterns. This is a consequence of 
having zeros or quasi-zeros in the OTFs at some specific spatial frequencies which 
produce a poor SNR at these frequencies. Thus, the division in equation 1.19 
by the squares of these nearly zero values cause excessive amplification of some 
spectral components. To circumvent this drawback, a sort of speckle summation 
of various realizations close in time has been successfully used. This summation 
carries out compensation for spectral information gaps in the Fourier domain. From 
equation 1.19 and summing up K realizations, one can easily derive 

K / \ K 



^2[Iki(u)Sl 1 (u,ex)+-fI k2 {u)Sl 2 {u,cx) j = f (w)^( \S k i{u,cx)\ 2 +-f\S k2 {u,cx) 

(1.21) 



k=l v ' k=l 
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where I (u) is not affected by the fc-subscript because it must be unique, i.e. the 
estimate of the true signal that we are seeking 5 . 

1.5 Multi-Object Multi- Frame Blind Deconvolution 

The MFBD (Multi-Frame Blind Deconvolution) is a restoration method that uses 
multiple frames, bringing in such a way complementary information to recover the 
aberrations affecting the images (MFBD; Lofdahl, 2002, 1996). The method gen- 
erally works in a better way when the contrast is high, the exposure time is short 
and the noise is low. 

An extension of the MFBD, named Multi-Object Multi-Frame Blind Deconvo- 
lution (MOMFBD; Van Noort, Rouppe van der Voort & Lofdahl, 2005), has made 
possible to use, apart from multiple frames, also multiple objects simultaneously ob- 
served, to restore the images in a more efficient way. By different objects we mean 
the same field-of-view in the Sun but observed in different wavelengths (within a 
rather narrow spectral range). For each object, PD focus-defocus pairs can also be 
included. The observations in all channels must be simultaneous so that we can 
assure a common atmospheric aberration. The joint restoration of several objects 
has also the advantage that almost perfect alignment can be achieved between all 
of them. 

Figure 1.3 shows a sketch of the optical setup at the Swedish Solar Telescope 
(Observatorio Roque de los Muchachos). The MOMFBD technique is applied inde- 
pendently to the blue and the red channels. The G-band and G-cont beams in the 
figure are splitted into two channels (focus-defocus) that produce PD image-pairs. 
Red channel has also its own PD-pair 6 . 

Throughout this thesis work we employ all the restoration methods described 
above (i.e. PD, MFBD and MOMFBD) depending on the particular circumstances 
for each situation. The high-spatial-resolution achieved by using these methods is 
evident from the results shown in all the different chapters. 



J For more detailed information on the PD-method see e.g. Bonet & Marquez (2003); Criscuoli 
et al (2005) and references therein. 

G The detailed description of the setup will be commented in the next chapter when dealing 
with the observations we pursued in different campaigns. 
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FIGURE 1.3 — Optical setup for the red (left) and the blue (right) beam at the SST showing the 
configuration of the optical elements and cameras in the proper disposition to pursue corrections 
with MOMFBD. Arrows mark the direction of the incoming light. Cameras marked as 1.6, 4.2, 
and 6,3 indicate the number of mega-pixels. ND stands for a neutral density filter placed along 
the red beam to balance the exposure level. The SOUP filter is a tunable Lyot filter with 72 
mA bandpass for the Fe I 6302 A line and 128 mA bandpass for Ha. Taken from Van Noort, 
Rouppe van der Voort & Lofdahl (2005). 
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In-flight calibration of IMaX 
aberrations 

One of the main concerns of this thesis has an instrumental nature. This chap- 
ter is devoted to present the assignments I have developed in this field as a 
member of the Imaging Magnetograph experiment (IMaX) project at the IAC in 
Tenerife. Initially, a brief preliminary part is presented to introduce the balloon- 
borne SUNRISE project conceived to pursue high-resolution solar observations and 
from which IMaX is part of. Afterwards, I will introduce the IMaX concept. Fi- 
nally, I will explain with great detail the specific part I have performed within the 
IMaX project that consists in designing a robust method to calibrate the in-flight 
instrumental aberrations of IMaX. 



2.1 SUNRISE 

The SUNRISE project consists of a balloon-borne 1-m aperture solar telescope 
(Figure 2.1) that aims at high-resolution spectro-polarimetric observations of the 
solar atmosphere, to be flown in the framework of NASA's LDB (Long Duration 
Balloon) program in 2009, in a series of flights on circumpolar trajectories at a float 
altitude of 35 - 40 km. The main goal of the project is understanding the formation 
of magnetic structures in the solar atmosphere and their interaction with the plasma 
convective flows. The SUNRISE equipment consist of the main telescope feeding 
two focal-plane instruments through a light distribution system. These instruments 
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have been developed by several PI institutions : 

• Telescope: Gregorian design. Carbon fiber based telescope structure with lm 
Schott Zerodur lightweighted primary mirror. 

• ISLiD: Image Stabilisation and Light Distribution System that ensures ca- 
pability of simultaneous observations with all science instruments, based on 
all-dielectric dichroic beam splitters (see Figure 2.2). An important part of 
the ISLID is the so-called Correlator and Wavefront Sensor (CWS) system 
in charge of image stabilisation and correction of optics misalignments and 
defocus. 

• SUFI: The Sunrise Filter Imager is a filtergraph for high-resolution images in 
the visible and UV spectral ranges. 

• IMaX: Imaging Magnetograph experiment, is a magnetograph providing fast- 
cadence two-dimensional maps of the complete magnetic field vector and the 
line-of-sight velocity as well as white-light images with high spatial resolution. 

With all this equipment, SUNRISE will provide spectra and images resolving spatial 
scales down to 35 and 70 km on the Sun in the ultra-violet and visible ranges, 
respectively. Figure 2.2 sketches the concept of the SUNRISE project. The main 
telescope in panel (b) puts the light beam into the optical bench in panel (a) 
being M4 the interface between both subsystems. The light beam is folded by 
the tip/tilt mirror and afterwards splitted into two channels, one feeding the CWS 
for wavefront sensing and the other reaching a dichroic plate that separates the 
specific wavelength for SUFI and IMaX, respectively. The CWS generates electric 
signals: SI for steering the tip/tilt mirror and S2 to control M2 in order to correct 
for defocus and optical misalignments. 

2.2 The Imaging Magnetograph experiment 

As briefly explained in the last section IMaX is one of the instruments part of the 
payload of the SUNRISE balloon. It will allow to study the dynamics and evolution 
of the solar magnetic field as well as its interaction with the plasma, with high 
temporal, spatial and spectral resolutions and unprecedented polarimetric sensitiv- 
ity. This instrument will then provide magnetograms of extended solar regions by 

1 Max Planck Institute for Solar System Research (MPS), Instituto de Astroffsica de Canarias 
(IAC), Kiepenheuer-lnstitut fur Sonnenphysik (KIS), High Altitude Observatory (HAO), Lockheed- 
Martin Solar and Astrophysics Laboratory (LMSAL). 
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FIGURE 2.2 — Box diagram showing the SUNRISE configuration, a) Light distribution system 
(ISLiD) and the two focal instruments (IMaX and SUFI), b) Telescope and folding mirrors. 
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combining high temporal cadence and polarimetric precision while preserving the 
bidimensional integrity of the images. To meet this goal IMaX has to work as a: 

• High-efficient image acquisition system. 

• Near diffraction limited imager. 

• High resolving power spectrograph. 

• High sensitivity polarimeter. 

Figure 2.3 shows the final optical and mechanical design of the instrument in a 
3D representation; the main optical components will be described in detail in the 
next section. Table 2.1 lists the main parameters of the instrument. 

2.2.1 Optical description 

The box diagram in Figure 2.4 shows a concise description of the IMaX optical 
configuration 2 (see figure caption). 

Figure 2.4 (lower panel) shows a more detailed view of the IMaX optical layout 
that will be briefly commented in the next lines. The optical interface with SUNRISE 
is the ISLiD Focus F4. Next to F4 we locate the prefilter (PF) set (bandwidth lA) 
and the modulator based on Liquid Cristal Variable Retarders (LCVRs or ROCLIs). 
After suffering the selected polarization, the beam passes through a collimator 
system consisting of 2 lenses and a doublet. In the collimated space we locate one 
solid Fabry-Perot Interferometer (which cavity is LiNbC>3 3 ) in the position of a pupil 
image of the SUNRISE system. The etalon operates in double pass which means 
that the light passes once through it, then reflects back by means of two folding 
mirrors and crosses for the second time the etalon. The beam is then focused by 
an imaging optical system (camera optics) that consist of a doublet and two lenses. 
A cubic polarizing beamsplitter divides the light beam into two branches, to form 
in the CCDs respective images with orthogonal polarization states. A parallel glass 
plate can be optionally inserted in the light path between the BS and one of the 
CCDs for calibration purposes by applying the post-facto phase diversity technique. 
We will extensively detail the phase diversity mechanism in section §2.3. 

2 From the IMaX Final Optical Design document SUN-IMaX-RP-IX200-023. 
^Due to its unique electro-optical, photoelastic, piezoelectric and non-linear properties Lithium 
Niobate (LiNbOs) is widely used in a variety of integrated and active acousto-optical devices. 
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TABLE 2.1 — Main IMaX parameters. 



PAR AMFTFR 


\/A IMF 


Avpercure 


1 m 

i m 


cTieciive TOCdi lengin 


ar nn m 
^fo.uu m 


PlXPl ^I7P 

1 1 y\ \_ 1 Jl 


12 ^/m 


Spatial sampling 


0.055 arcsec/pix 


Working wavelength A 


525.02 nm (Fel) 


Pixels 


1024 x 1024 


FOV 


50 x 50 arcsec 


Image exposure time 


< 200 ms 


IMaX weight 


< 80 Kg 


IMaX data rate 


< 850 KB/s 


SPECTROPOLARI METRIC MODE 




Number of wavelengths 


4 + 1 continuum 


Time for l,Q,U,V 


< 30 sec SNR 1000 


Time for l,V 


< 15 sec SNR 1000 


MAGNETOGRAPH MODE 




Number of wavelengths 


4 + 1 continuum 


Time for l,Q,U,V 


< 20 sec SNR 1000 


Time for l,V 


< 10 sec SNR 1000 


DEEP MAGNETOGRAPH MODE 




Number of wavelengths 


1 


Time for l,V 


< 160 sec for SNR 1000 
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FIGURE 2.3 — Sketch showing the IMaX design with all its components, including among others 
the etalon working in double-pass, the two science cameras and the box devoted to the proximity 
electronics. The two CCD cameras will also be employed in combination with a plane-parallel 
glass plate to obtain bursts of PD image-pairs to calibrate the in-flight optical aberrations. All the 
information about the project can be found in the website http://www.iac.es/proyect/IMaX/. 



2.3 The Phase Diversity plate 

As commented in the last section, the IMaX design includes a mechanism to apply 
the Phase Diversity technique for calibration purposes (wavefront sensing). Ac- 
cording to section §1.4 where the PD method is described, we need to combine 
information included in at least two simultaneous images of the object, one be- 
ing the conventional focal-plane image that is degraded by the unknown system 
aberrations we are interested in, and the other one affected by the same unknown 
aberrations plus an extra known aberration intentionally induced. 

For the sake of simplicity this extra aberration is commonly chosen as a defo- 
cus (see Appendix A) induced by simply displacing one of the cameras out of focus 
by a certain known distance that we will also term diversity. For IMaX another 
alternative was implemented since the displacement of the camera could generate 
inertial problems caused by redistribution of masses or vibrations when moving a 
heavy CCD. One alternative option to induce defocus in one of the images with- 
out moving the CCD camera to avoid the above mentioned problems, consist on 
intercalating a light plane-parallel glass plate in front of one the CCDs to induce a 
known displacement of the image focus as will be described below. 

The defocus aberration can be described as an excess/defect in the radius of 
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FIGURE 2.4 — Upper panel: IMaX box diagram showing the direction of the incoming light 
(IL) from the Telescope and ISLiD and the focus position (F4) just before entering IMaX. The 
main IMaX optical components are: the prefilter and the Liquid Cristal Variable Retarders (PF & 
LCVRs), the collimator (COLL), the etalon (E) located in a conjugate pupil and working in double 
pass, the camera optics (CAM), the cubic polarizing beamsplitter (BS) that divides the light beam 
into two branches, to form in the CCD respective images with orthogonal polarization states. 
In one of the branches, a phase diversity (PD) device can be optionally inserted for calibration. 
Lower panel: IMaX optical layout. The light goes from focus F4 to the CCD cameras after 
passing through all the optical elements drawn in the re-scaled diagram. ROCLIs are the same as 
LCVRs. Based in the document IMaX Final Optical Design, SUN-IMaX-RP-IX200-023. 
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FIGURE 2.5 — Geometrical configuration for a defocus system. The sketch shows the corre- 
spondence between the PV defocus aberration ds and the longitudinal displacement of the image 
plane, Az, along the optical axis. 

the spherical wavefront with respect to the value that would be required to form the 
image of a point source at the nominal focus F (see Figure 2.5). The subtraction 
between both wavefronts give us the defocus aberration function (f> that turns out 
to be a paraboloid function. 

Geometric considerations lead to the following relationship between the axial 
displacement of the image plane (Az) and the defocus aberration function 



where D is the aperture, < p < 1, pD/2 is the radial distance from the pupil 
center, A is the working wavelength and / the effective focal length of the system. 

On the other hand, the defocus aberration function expressed as a Zernike 
polynomial in the Noll's basis (see Appendix A) is 




(2.1) 



a4\/3(2p 2 — 1) = 2a^V3p 2 — a^V3, 



(2.2) 



where is the defocus weighting coefficient (in radian). The last term a4\/3 in 
equation 2.2 represents an offset that is not affecting the shape of the wavefront 
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but only an axial shift (piston). Thus we can relate expressions in equations 2.1 
and 2.2 as 



ttAz /L> x 



\ p 2 = 2a 4 VV , (2.3) 



4A V/ 

and finally obtain the displacement in terms of the coefficient 04 (in radian): 

Az = -—{d) " 4 - (2 - 4) 

By replacing the IMaX working parameters (Table 2.1) in the last expression 
we obtain Az as proportional to <x 4 as follows: Az = - 4.68923 0,4 mm. 

This relation will be used later on to have an intuitive and also quantitative 
feeling about how significant is a defocus aberration expressed as a certain a± value. 

Apart from the axial displacement Az, another way to characterize the amount 
of defocus is by means of the Peak-to-Valley (PV) optical path difference (OPD) 
which obviously corresponds to the value of the defocus aberration function at the 
pupil edge. According to expression 2.1 this value is 

2ttAs hAz { D x ' 2 



or in wavelength units 



a saw;- (2 - 6) 

Table 2.2 presents the defocus aberration in terms of the PV optical path differ- 
ence corresponding to the focus displacement along the optical axis for a telescope 
with /Number = f/D = 45.0 and A525.02 nm, which are the adopted values in 
IMaX. 



pv = — = ~sr I7' ' (2 - 5) 



As Az fD^ 



As a rule of thumb, an amount of defocus corresponding to As = 1 A at the 
pupil edge has proved to be satisfactory for the phase diversity inversions. Using 
the IMaX parameters (see Table 2.1), and by setting As = 1 A, we obtain a value 
for the focus displacement along the axis Az = 8.51 mm. 
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Optical Path Difference (OPD) 
at the pupil edge (As/A) [number of waves] 


Focus displacement along the optical axis 
Az [mm] for /|\|umber = 45 ar| d A525.02 nm 


1/18 


0.472518 


1/15 


0.567022 


1/14 


0.607523 


1/13 


0.654256 


1/12 


0.708777 


1/11 


0.773211 


1/10 


0.850532 


1/9 


0.945036 


1/8 


1.06317 


1/7 


1.21505 


1/6 


1.41755 


1/5 


1.70106 


1/4 


2.12633 


1/3 


2.83511 


1/2 


4.25266 


1 


8.50532 



TABLE 2.2 — Relation between the optical path difference (OPD) at the pupil edge in wave- 
length units and the focus displacement (Az) along the optical axis for a particular telescope and 
wavelength A. 
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As mentioned above, we will induce the image plane displacement by employing 
a glass plate of thickness d, as shown in Figure 2.6. This displacement can be 
mathematically formulated by the expression 

Az = d(l-^\, (2.7) 

where n is the refractive index of the material. The thickness of the plate re- 
quired to produce a desired displacement depends on the refractive index of the 
optical glass commercially available, and for IMaX the selected material was Fused 
Silica 4 with n=1.461. Assuming a PV aberration of 1 A or equivalently Az=8.51 
mm as reported above, the required thickness for our PD plate would be d=27 mm. 



FIGURE 2.6 — Displacement on 
axis produced by a plane parallel 
plate. An incident ray forming an 
angle a with the normal to the 
surface is deflected when entering 
the plate and eventually leaves the 
plate with the same angle but in- 
ducing an OPD as shown in the 
figure. 





n 










-H — — ^ 




d 


47 


< > 



Defocus tolerance 

The image blurring permitted for an instrument can be specified by the diameter 
of the blur spot or as the angle (/3) subtended by such a diameter. For a given 
(3 of acceptable blurring, the focus deepness is lower towards the lens that in the 
backward direction. 



4 Fused Silica is a high purity synthetic amorphous silicon dioxide. This noncrystalline, colorless, 
silica glass combines a very low thermal expansion coefficient with excellent optical qualities and 
exceptional transmittance over a wide spectral range and is also resistant to scratching and thermal 
shock. 
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A criterion for the defocus tolerance can be made by fixing the blur angle (3 per- 
mitted for our purposes. For instance, we can select as the value of the diffraction 
cutoff wavelength, that is slightly greater than the Full-Width at Half-Maximum 
(FWHM) of the Airy spot. According to this criterion the defocus tolerance Az t 
can be expressed as 



Nevertheless, this criterion is quite severe and there is another less-strict and 
more commonly used to establish the limit of the defocus tolerance in an optical 
system. This criterion is based on the loss of intensity in the central part of the PSF. 

Rayleigh 5 first (for the spherical aberration) and other authors afterwards (for 
other aberrations), demonstrated that when the PV value of the error in the wave- 
front is < A/4, the central intensity of the PSF of the system reduces by less than 
20% with respect to the central intensity of the Airy function. 

Marechal demonstrated that when the rms of the wavefront aberration in an 
optical system is < A/14, the central intensity of the PSF also reduces by less than 
20%, or in other words, the Strehl ratio 6 is > 0.8. 

Both limits described above, A/4 for the PV value or A/14 for the rms, are com- 
monly accepted as tolerable when establishing whether or not a system achieves 
the resolution limit for practical purposes. These limits are referred to as the A/4- 
Rayleigh and the Marechal image quality criteria, respectively. 

Born & Wolf (1999) study the distribution of intensity in a volume around the 
focus in a diffraction limited optical system and, in particular, the intensity in the 
central part of the PSF along the optical axis. In this study it is demonstrated that 
to reduce the central intensity of the PSF by a 20%, the displacement as measured 
from the nominal focus is given by the equation 



5 Lord Rayleigh (John William Strutt), 1842-1919, Nobel prize in Physics in 1904. 

6 The Strehl ratio of an optical system is defined as the quotient between the central intensities 
of the real PSF and the theoretical one assuming a diffraction limited system, the latter being the 
Airy function. SR =PSF(0, 0) /Airy(0, 0). This ratio is closely related to the sharpness criteria 
for optics defined by Karl Strehl. 




(2.8) 
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) 



2 



Az t = ± 2A 



(2.9) 



This is the classical limit adopted as the tolerance for the defocus, an it is based 
on the above described Rayleigh and Marechal criteria. Note that the difference in 
both defocus tolerances in equations 2.8 and 2.9 is a factor of 2, meaning that the 
last one is more permissive as a quality criterion. 

With the IMaX parameters (in Table 2.1) the defocus tolerances Azt for the 
instrument according to both defocus criteria are, respectively, 



Nevertheless, the modelled defocus error in IMaX, stemming from optical tol- 
erancing and thermal effects, remains within the tolerance range even in the more 
restrictive case (± 1.06 mm). 

2.4 Defining the method for the in-flight calibration of the 
IMaX image aberrations 

The goal of this section is to define the procedure that will be employed in IMaX to 
calibrate the image degradation during the flight. The proposed calibration proce- 
dure is based on the PD-speckle technique for image reconstruction and wavefront 
sensing (see section §1.4). To that aim IMaX has been provided by a PD-plate to 
intentionally induce a controlled defocus. 

The PD-plate can be optionally intercalated in one of the IMaX imaging chan- 
nels (see Figure 2.7 and section 2.3) so that a simultaneous focus-defocus image- 
pair, i.e. a PD image-pair, can be recorded. From this pair an estimate of the 
aberrations will be possible in post-processing by means of a PD inversion code. 
PD-inversions throughout this chapter have been performed with an IDL code de- 
veloped by Bonet & Marquez (2003). 

Assuming a long-term variation in the instrumental aberrations, the image ac- 
quisition for calibration could be performed with a cadence of one hour. A burst 




(2.10) 




(2.11) 



28 



CHAPTER 2. In-flight calibration of IMaX aberrations 




of 25-30 PD-pairs taken in a short time interval in the spectral continuum would 
be enough each time. This way we get a large amount of information 7 collected 
within a sufficiently short time period so that it can be assumed negligible evolution 
of the solar structures in the FOV. The averaged results from the PD-inversions of 
these calibration image-pairs in a post-facto process will provide a description of the 
aberrations also affecting our science spectral images and, in turn, the maps of the 
full-Stokes vector. Subsequently, the science images will be reconstructed by us- 
ing a standard deconvolution code to nearly reach the diffraction limit of SUNRISE. 

The organogram of Figure 2.8 shows a schematic representation of the sequen- 
tial steps in the post-processing calibration procedure. 



2.5 Testing the robustness of the calibration method 

The purpose of this section is to evaluate by means of numerical simulations, the 
robustness of the method we have implemented for the wavefront error calibra- 
tion in IMaX, versus a variety of aberration assumptions. To pursue the numerical 
experiments we select an isoplanatic patch consisting of a portion of one of the 
synthesized images of solar granulation, from the MHD model developed by Vogler 
et al. (2005); Schussler et al (2003). The spatial resolution of this synthetic image 
is much higher than that corresponding to the cutoff frequency in SUNRISE 8 , and 

improving the SNR. 

8 ln order to be as realistic as possible constructing the theoretical images one has to consider 
the spectral contamination by the wings and lobes of the filter transmission curve T(X). For this 
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FIGURE 2.8 — Organogram explaining the strategy for the calibration procedure that will be 
applied to the IMaX data. Within a time interval of ~ 10 s a burst of focus-defocus image-pairs 
are recorded and labeled from 1 to K. Every single image is divided in a total of N x M square 
isoplanatic patches. Note in the figure that the patches are overlapped since their borders will be 
corrupted by mathematical manipulations in the Fourier domain during the PD-inversion process. 
For all the focus-defocus image-pairs the PD inversions are independently performed for every 
single patch-pair (corresponding to a common focus-defocus portion in the FOV) to obtain a set 
of L — 1 Zernike coefficients (from 02 to q_l) describing the wavefront aberration affecting such a 
portion of the FOV. Thus, for instance, the nomenclature af nm stands for the Zernike coefficient 
of index 4 (defocus) computed from the patch centered at coordinates (n,m) in the last image-pair 
recorded which is labeled by K. The coefficients extracted from every patch are then averaged for 
all the K realizations and hereby we are able to describe the wavefront aberration in any given 
patch of the FOV. Interpolation over the whole FOV will result in continuous maps for every 
aberration coefficient ai(x,y), with 2< I < L. The OTF of the system for every location on the 
image is straightforward computed from these cti(x, y) functions (eqs.1.18 and 1.15). 
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represents what will be termed hereafter the true object. 

The testing procedure will consist in simulating the formation of PD image-pairs 
as produced by a telescope with the SUNRISE aperture and a given set of aberra- 
tions (model of aberrations). A total of 30 image-pairs for different photon noise 
realizations is hereby obtained. The resulting degraded image-pairs are re-sampled 
to the IMaX pixel size and inverted with the PD-code. This set of images will be 
processed according to the diagram sketched in Figure 2.8 applied to the particular 
case of only one isoplanatic patch, and the set of averaged aberrations retrieved 
from the inversions will be compared to the input model of aberrations. 

In our simulations we intend to reproduce the global aberration as close as pos- 
sible to the real one affecting the image quality when the instrument is acquiring 
data in real-time. For this purpose we need to, firstly, identify the potential sources 
of error and, secondly, make realistic assumptions about the quantitative contribu- 
tion from each of them. 



2.5.1 Identifying the error sources 

The contribution from the different error sources can be mathematically represented 
through the generalized pupil function (see equation 1.18) as follows, 



H(p, 9) = \H(p, 9)\ exp i [5(p, 9) + <j> Hp pie(p, 0) + Mp, &) + <t>(P, 0) + <f> atm (p, 9)] 



where the terms in the exponential stand for: 



purpose a sequence of images are constructed from the theoretical model for different A's around 
ourworkingA: i(x,y,X). This sequence has been convolved with T(A): i'(x, y, X)=i(x, y, A)*T(A) 
and from the resulting sequence i'(x,y,\) the central image in our working A is extracted. These 
calculations have been done by D. Orozco Suarez. 



(2.12) 
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Differential aberration between focus/defocus channels = phase diverse. 
Main mirror polishing error. 
Phase error from the etalon. 

Low-order optical system aberrations common in focus/defocus channels. 
Atmospheric aberration which is negligible in IMaX. 
Transmission function over the pupil. 

In addition we have to include in the error budget the pixel integration effect 
in the CCD and the noise. 

6(p,0) and 4>{p,0) will be expressed as a polynomial expansion in a basis of 
Zernike as explained in section §1.4. The rest of the contributions will be intro- 
duced as amplitude- or phase-screens constructed according to a model or obtained 
from direct laboratory measurements. 

We emphasize that the formulation of the Zernike's basis we use corresponds 
in all cases to that introduced by Noll (1976) as stated in section §1.4. 



6(p,9) 

4>ripple(P, 0) = 

4>atm(p,0) = 

\H(p,0)\ = 



2.5.2 Quantifying the contribution of the error sources 

The opto-mechanical elements in SUNRISE which are affecting the formation of our 
images are: Telescope, ISLiD, mechanical interface with IMaX and IMaX 
itself (see Figure 2.2 for a schematic view of the optical design including all these 
components). The first step is, consequently, the compilation of data from the 
design and specifications of all the different components 9 . 

Table 2.3 lists a compendium of some nominal values for the rms wavefront 
error (rms-WFE) derived from the optical design including the statistical modelling 
for the optical tolerancing and the thermal behaviour of the systems 10 . 

9 Most of these data have already been published in the IMaX technical document: SUN-IMaX- 
RP-IX200-023 IMaX Final Optical Design (2007). 

10 Thermal effects are simulated with an optical design software (CODE V) by varying the refrac- 
tive index and the geometry of the opto-mechanical components as a function of the temperature 
and the expansion coefficients. Temperature variations, for instance, are taken into account in 
the expected range except in those components having their own temperature control in a lower 
range (i.e Liquid Cristal Variable Retarders - LCVRs -, etalon, etc). Simulations are meant to 
study the aberration for three possible cases: operative, hot and cold. It is mandatory to know 
the thermal variation of the refractive index (dn/dT) for every optical material and the thermal 
expansion coefficient for every opto-mechanical element. The thermal effects have been studied 
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The rms values shown hereafter are computed excluding piston and image dis- 
placement (tip/tilt) from the error budget. 



Low-order optical system aberrations 

In addition to the above mentioned mean parameters describing the rms aberration 
in Table 2.3, we also have available the sets of Zernike coefficients corresponding 
to the nominal design for both, Telescope+ISLiD and IMaX, up to a total of 
36 modes (7 th order) starting from the term of index 4, i.e. defocus (see columns 
2 and 3 in Table 2.4). Respective rms-WFE are included in Table 2.3 (1/26 and 
1/140 waves). 

Initially, at a first stage we represented the model for low-order aberrations 
in our numerical experiments by re-scaling the coefficients of the nominal design 
by a factor accounting for the thermal effects and optical tolerancing 11 , excluding 
the etalon (the contribution of which is explicitly introduced in the simulations 
by means of amplitude- and phase-screens supplied by the manufacturer) and the 
ISLiD-IMaX interface contribution. In Table 2.3 the re-scaling factors correspond 
to: Telescope+ISLiD (with optical tolerancing and thermal effects, item 2) and 
IMaX (with optical tolerancing and thermal effects, item 7), so that their values 
are 1/8.25 and 1/19, respectively. The resulting sets of re-scaled Zernike coeffi- 
cients were summed up to give the desired low-order aberrations model (column 4 

for IMaX as a whole and also for its main components independently (e.g. the PD-plate, etc). 
The calculations were made under the assumption of a temperature variation range of 15±7.5°C. 
However, further improvements in the thermal system (May 2008) have provided a more restrictive 
range for the temperature excursions: 24.5±0.5°C. This implies that the rms-WFE values in the 
table, derived from the former thermal model, are rather pessimistic estimations. Nevertheless, 
we preserved them because new estimations are not available for all cases yet. Moreover, testing 
the robustness of our calibration method with over-degraded images will reinforce the conclusions 
about the reliability of the calibration method. Optical tolerancing is a critical step in the design 
of an optical system. The objective is to define a fabrication and assembly tolerance budget and to 
accurately predict the resulting as-built performance, including the effects of compensation. Also 
part of the study is determining the best set of compensators. 

n Let {cij,j = 1, 2, J} be the set of Zernike coefficients of the expantion in the Noll's basis, 
approaching a wavefront aberration. The orthogonality properties of these basis functions facilitate 

the calculation of the rms-WFE as \JY1 a |- Re-scaling the coefficients to provide a wavefront 

aberration with a given {rms-WFE) 1 consists of simply computing a new set of coefficients a' ; such 
that a- = a^rms-WFE)' /{rms-WFE). 
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TABLE 2.3 — Aberrations budget from the modelling of different optical systems in SUNRISE. 
Polishing errors in the main mirror are not included. 







rms-WFE 


ITEM 


SYSTEM 1 " 


(no tip/tilt) 
[waves] 


1 


Telescope+ISLiD (nom. design) averaging 9 points in FOV 


1/26 


2 


Telescope+ISLiD (nom. design+op.tolerancing+therm.eff.)* 


1/8.25 


3 


Interface ISLiD-IMaX (misalignment) 


1/42 


4 


IMaX without etalon (nominal design) 


1/140 


5 


IMaX without etalon (nominal design + optical tolerancing) 


1/22 


6 


IMaX without etalon (nominal design + thermal effects) 


1/38 


7 


IMaX without etalon (nom. design+ op.tolerancing+therm. eff.) A 


1/19 


8 


Etalon (lab. calibration CSIRO ) 


1/23 


9 


IMaX (nom. design + op. tolerancing + thermal eff. + etalon) 


1/14.6 


10 


Total: Teles.+ISLiD+Interface+IMaX® 


1/7.1 



f Symbol "+" is employed to include another element. Measured by the etalon's manufacturer. 

| RSS = Root Sum Square: □ Equivalent to the RSS of items 7 and 8. 

rms-WFE = ^/(rms-WFEi) 2 + (rms-WFE 2 ) 2 . © Equivalent to RSS of items 2, 3 and 9. 
A Equivalent to RSS of items 5 and 6. 
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TABLE 2.4 — Error budget for different sub-systems in SUNRISE. Zernike coefficients in Noll's 
basis. 







IMaX 


Tel-|-ISLiE> 


T -|- 1 S -|- 1 1V1 ^ 


IMaX 


Xel.+ISLiD 


Sum prev. 


Etalon aberr. 


TOTAL 


Zj 


no etalon 




Sum of re-scaled 


no etalon 


re-scaled 


two columns 


Low-order modes 


Sum prev. 




nom.aberr. 


nom.aberr. 


nom.aberr. 


Zygo 


nom.aberr. 




Zygo 


two col. 




[rad] 


[rad] 


[rad] 


[rad] 


[rad] 


[rad] 


[rad] 


[rad] 


4 


0.0445282 


-0.0714395 


0.0493187 


0.0411549 


-0.278658 


-0.237503 


0.647168* 


-0.237503 


5 


-1.24691e-09 


2.51327e-07 


9.71148e-07 


0.0157080 


9.80332e-07 


0.0157089 


-0.0565487 


-0.0408397 


6 


-0.00246226 


0.180514 


0.685980 


0.279602 


0.704116 


0.983718 


0.175929 


1.15965 


7 


-1.11203e-09 


-0.000384531 


-0.00149992 


0.00314159 


-0.00149991 


0.00164168 


0.0314159 


0.033057 


8 


0.000775724 


0.0186024 


0.0782744 


-0.0408407 


0.0725608 


0.0317201 


-0.0753982 


-0.0436782 


9 


1.52930e-09 


0.000275518 


0.00107470 


0.00314159 


0.00107469 


0.00421628 


-0.0502655 


-0.0460492 


10 


3.28740e-07 


-0.000634476 


-0.00247243 


0.0219912 


-0.00247485 


0.0195163 


-0.0188496 


6.66745e-4 


11 


-0.00423551 


0.00933983 


0.00523408 


0.0219912 


0.0364311 


0.0584223 


-0.0188496 


0.0395727 


12 


-6.63747e-08 


-0.000165059 


-0.000644322 


0.0251327 


-0.000643833 


0.0244889 


-0.0125664 


0.0119225 


13 


6.54083e-10 


6.28319e-08 


2.49901e-07 


0.00000 


2.45083e-07 


2.45083e-07 


0.0314159 


0.0314162 


14 


-6.65045e-07 


2.51327e-06 


4.90487e-06 


-0.00628319 


9.80332e-06 


-0.00627338 


0.0125664 


0.00629299 


15 


1.23446e-09 


0.00000 


9.09256e-09 


-0.0125664 


0.00000 


-0.0125664 


0.00628319 


-0.00628319 


16 


2.13090e-09 


-1.54566e-05 


-6.02747e-05 


0.0219912 


-6.02904e-05 


0.0219309 


0.0125664 


0.0344972 


17 


4.98823e-10 


-4.39823e-07 


-1.71191e-06 


0.0785398 


-1.71558e-06 


0.0785381 


-0.00628319 


0.0722549 


18 


7.93558e-10 


1.88496e-07 


7.41094e-07 


-0.0219912 


7.35249e-07 


-0.0219904 


0.00628319 


-0.0157072 


19 


-7.84298e-10 


1.25664e-07 


4.84389e-07 


0.00628319 


4.90166e-07 


0.00628368 


-0.0125664 


-0.00628270 


20 


5.43063e-09 


0.00000 


3.99998e-08 


0.0282743 


0.00000 


0.0282743 


0.00628319 


0.0345575 


21 


1.04887e-09 


0.00000 


7.72552e-09 


-0.0125664 


0.00000 


-0.0125664 


0.0314159 


0.0188496 


22 


0.000234998 


5.34071e-06 


0.00175173 


0.185354 


2.08321e-05 


0.185375 


0.0314159 


0.216791 


23 


-6.92688e-ll 


0.00000 


-5.10206e-10 


0.0691150 


0.00000 


0.0691150 


-0.00628319 


0.0628319 


24 


6.84204e-09 


0.00000 


5.03957e-08 


0.0282743 


0.00000 


0.0282743 


0.0188496 


0.0471239 


25 


1.05472e-09 


0.00000 


7.76868e-09 


0.0188496 


0.00000 


0.0188496 


-0.00628319 


0.0125664 


26 


-7.45863e-07 


0.00000 


-5.49373e-06 


0.00314159 


0.00000 


0.00314159 


0.00000 


0.00314159 


27 


-1.08106e-09 


0.00000 


-7.96267e-09 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 


28 


5.19823e-09 


0.00000 


3.82881e-08 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 


29 


-6.84537e-10 


0.00000 


-5.04203e-09 


0.0691150 


0.00000 


0.0691150 


0.0188496 


0.0879646 


30 


-5.44505e-09 


0.00000 


-4.01060e-08 


0.0471239 


0.00000 


0.0471239 


0.00628319 


0.0534071 


31 


-3.69741e-10 


0.00000 


-2.72336e-09 


0.0188496 


0.00000 


0.0188496 


0.00628319 


0.0251327 


32 


-1.26343e-08 


0.00000 


-9.30591e-08 


0.0125664 


0.00000 


0.0125664 


0.00628319 


0.0188496 


33 


1.23635e-09 


0.00000 


9.10646e-09 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 


34 


5.01911e-09 


0.00000 


3.69687e-08 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 


35 


1.78661e-09 


0.00000 


1.31595e-08 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 


36 


-8.78520e-09 


0.00000 


-6.47082e-08 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 



A Telescope+ISLiD+IMaX. 

4b This contribution is set to zero (— * 0.00000) for all our calculations hereafter. See text for detailed information. 
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TABLE 2.5 — Aberrations budget (rms-WFE) from laboratory calibrations. 







rms-WFE 


ITEM 


SYSTEM^ 


(no tip/tilt) 






[waves] 


1 


IMaX excluding etalon and PD-plate D . From Zernike coeffs. 


1/16.8 


2 


IMaX (item 1) + (Telesc.+ISLiD) A . From Zernike coeffs. 


1/6 


3 


Etalon in double-pass: From the CSIRO phase-screen* 


1/26 


4 


Etalon in double-pass within the oven^. From Z coeffs. 


1/9.2 


5 


Etalon contribution in item 4 after removing defocus 


1/28.5 


6 


Total: Low-order aberration terms' 8 . From Zernike coeffs. 


1/5.2 



f Symbol "+" is employed to include another element. 
□ Laboratory calibration with Zygo at INTA. 

A Nominal coefficients of Teles. +ISLiD, re-scaled with the factor in item 2 of Table 2.3. 
4k Optimal selection of two sub-apertures in the phase-screen provided by the manufacturer. 

Tip and tilt have been removed. 
At operative temp. (23°C) and voltage (2000 v). Lab. calib. Zygo, INTA. 
<8> (re-scaled Teles. +ISLiD)+(IMaX excluding PD-plate and considering in the etalon 

only the low-order aberration terms except for the defocus; Zygo). 
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in Table 2.4). 

A criticism to this procedure is that the simple re-scaling preserves the distri- 
bution curve of the nominal Zernike coefficients, which might be a rather arbitrary 
statement because the contribution from optical tolerancing and thermal effects 
could modify in a non-proportional way the different aberration terms. Having this 
in mind, we changed, according to our possibilities, the criterion to construct the 
model of aberrations affecting the images in IMaX. 

Table 2.5 complements Table 2.3 by adding new empirical measurements and 
other considerations regarding the wavefront errors. Once IMaX has been manufac- 
tured and assembled (excluding etalon and PD-plate), the aberrations have been 
calibrated in the laboratory by using a Zygo Interferometer at INTA. Figure 2.9 
shows the arrangement in the laboratory including IMaX and Zygo, sketching the 
optical elements and the light paths along the whole system when calibrating the 
optical aberrations of IMaX. The etalon and PD-plate drawn in the sketch have 
been actually removed for this calibration as mentioned above. 

From these measurements we obtained the Zernike coefficients approaching the 
real IMaX aberrations (excluding etalon and PD-plate) at laboratory conditions 12 
(Column 5 in Table 2.4). The corresponding value for the rms-WFE is 1/16.8 
waves (item 1 in Table 2.5). 

Concerning the Telescope+ISLiD contribution, we do not have the laboratory 
measurements and therefore we have preserved the procedure of using the coeffi- 
cients from the nominal design after rescaling by a factor accounting for the average 
sources of error (optical tolerancing and thermal effects; item 2 in Table 2.3) numer- 
ically modelled (column 6 in Table 2.4). Figure 2.10 (left panel) plots the Zernike 
coefficients for IMaX as measured at the laboratory and the Telescope+ISLiD co- 
efficients as re-scaled from the nominal ones. The total added contribution of both 
(Column 7 in Table 2.4) is plotted in Figure 2.10 (right panel) and represents the 
set of coefficients we are going to use in the simulations. In all cases we present the 
Zernike coefficients corresponding to the Noll's basis and starting from coefficient 
«4 (defocus). The terms for piston, tip and tilt are excluded from our model since 
they do not represent figure errors. 



12 Atmospheric pressure and room temperature about 20°C. At a first stage, laboratory experi- 
ments for checking the thermal behaviour of the system have not been performed. 
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FIGURE 2.9 — Optical configuration of the system including IMaX and the Zygo interferometer 
when calibrating the optical aberrations of the IMaX instrument (upper figure) and a snapshot of 
the Zygo software interface (lower figure). 
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FIGURE 2.10 — Plot of the Zernike coefficients representing the IMaX (excluding etalon and PD- 
plate) and Telescope+ISLiD degradation independently (left panel), and their added contribution 
(right panel). The coefficients correspond to the Noll's basis, starting from the term of index 4 
(defocus). 



The calculation of the rms-WFE from the quadratic sum of coefficients in our 
model produce a joint aberration of 1/6 waves (item 2 in Table 2.5) to be com- 
pared with the value RSS (=1/7.6 waves) derived from items 2 and 7 in Table 2.3 
that represent the global error effect (excluding the ISLiD-IMaX interface and the 
etalon) as modelled with CODE V 13 . The difference might be partially ascribed to 
the approximate behavior of the RSS parameter representing the joint aberration 
of various optical systems. 



Phase and transmission errors from the etalon 

As described before, the etalon is placed onto a conjugate pupil of the SUNRISE 
system. Nevertheless, up to this step, we are not including in the model of aber- 
rations the etalon contribution. The effects induced by the etalon in amplitude 
(transmission) \H(p,8)\ and phase (j) e (p,6) will be explicitly included in our simu- 



13 Software by Optical Research Associates used at INTA for designing and analysing the IMaX 
optical configuration. 
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lations. The manufacturer, CSIRO , provided computer files mapping the errors in 
amplitude and phase all over the etalon circular surface of 70 mm of diameter (this 
information is what we call amplitude- and phase-screens). Based on these maps 
we selected two optimal circular areas of 25 mm of diameter each. The criterion for 
this selection was to achieve the best compromise to minimize both the amplitude 
and phase errors. The etalon surface was oriented so that the selected circular areas 
fitted the position of the pupil image on the etalon operating in double pass mode. 
Considering the double pass through the selected areas we derived from the error 
maps supplied by the manufacturer, the amplitude- and phase-screens displayed in 
Figure 2.11 that exhibit predominately high-spatial frequency structures (high-order 
Zernike modes in the phase-screen). The resulting rms-WFE value is as small as 
1/26 waves (Table 2.5, item 3). The residual piston and tip/tilt were removed from 
the phase-screen. 

Nevertheless, the problem gets complicated when enclosing the etalon in the 
pressurized and thermalized oven (the thermally controlled enclosure). Interfero- 
metric measurements (with Zygo at INTA) of the wavefront error in double-pass 
and parallel-beam configuration, unveiled in the joint etalon+oven system a sig- 
nificant extra contribution of low-order Zernike modes (column 8 in Table 2.4 and 
Figure 2.12) with a rms-WFE of 1/9.2 waves (item 4 in Table 2.5). Since this value 
comes from the measured Zernike coefficients, that do not fundamentally describe 
the high-spatial frequencies detected in the calibration maps from CSIRO, we will 
assume that both contributions are complementary and additive and so we will 
treat them in our simulations. The origin of low modes is not fully-characterized 
yet though we speculate they are originated by deformations of the oven windows 
caused by mechanical stresses or by the bulge of the etalon itself when applying 
high voltages (ranging from to 2000 volts) while scanning our working spectral 
line. 

Note that the major contribution to low-order modes in the etalon corresponds 
to the defocus coefficient (~0.65 rad or 3 mm of displacement), and we are con- 
fident it can be compensated when coupling the etalon+oven system into IMaX 
by optimizing the position of the image focal plane 15 . For this reason this defocus 
contribution will be considered as null in our simulations (first value in column 8 
in Table 2.4), so that the rms-WFE ascribed to the etalon+oven in low modes 

14 Australian Commonwealth Scientific and Industrial Research Organisation, 
website: http://www.csiro.au/ 

15 By means of an MTFs optical bench. 
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Amplitude in double-pass: Phase in double-pass: 

\ll(f},0) * (; (P,0) 




FIGURE 2.11 — Amplitude and phase aberration screens for the etalon (rms <f> e {p, 8) — 1/26 
waves). 

shrinks from 1/9.2 to 1/28.5 (item 5 in Table 2.5). The rest of the Zernike coeffi- 
cients will be added to the low-order Zernike coefficients characterizing the system 
Tel+ISLiD+IMaX (column 7 in Table 2.4) thus resulting the total budget of low- 
order aberration terms affecting the image formation in IMaX (column 9 in Table 2.4 
and Figure 2.13) with a rms- WFE= 1/5. 2 waves (item 6 in Table 2.5). The total 
phase error induced by low-order aberrations will be denoted hereafter as 4>z{p,&)- 

In the final array of Zernike coefficients remains a non-zero value for defocus 
stemming mainly from ISLiD. This coefficient can not be set to zero in the sim- 
ulations as we did in the case of the etalon defocus since the ISLiD system has 
been designed and integrated with total independence with respect to IMaX. Thus, 
we have not any chance to compensate the defocus contribution from ISLiD when 
fixing the optimum position of the image focal plane during the IMaX integration. 

Main mirror polishing errors 

Another effect we are meant to include in our numerical experiment is the so-called 
ripple. Commented briefly in section §2.5.1 this error is caused by the polishing 
tool on the main mirror surface of the SUNRISE telescope. The wavefront error 
induced by this effect is quantified by a certain rms-WFE that will be, from now 
on, represented as rms-ripple and modelled by a phase screen arbitrarily chosen out 
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FIGURE 2.12 — Plot of the Zernike coefficients representing the etalon degradation contribution 
from every sub-aperture (left panel) and for the double-pass configuration (right panel). The 
coefficients correspond to the Noll's basis, starting from the term of index 4 (defocus). Values for 
the rms of every single sub-aperture are 1/18.8 and 1/16.4 waves, respectively, and the one for 
the complete double-pass configuration 1/9.2 waves. 



of a sample of realizations which average power spectrum matches a von Karman 6 
power spectrum, for an outer-scale equal to the size of the polishing tool (hereafter 
ripple-scale). The size assigned to the ripple-scale in our simulations is 30 cm. 
Note that the amount of rms-ripple we are considering is derived from the wave- 
front and consequently it is twofold the ripple error in an optical surface working 
by reflection. The total phase error induced by ripple will be denoted hereafter as 

4>ripple(P, 0)- 

The ripple can be classified as a high-order modes contribution to the global 
WFE. Figure 2.14 shows a ripple screen realization. 



10 Theodore von Karman (1881-1963), originally from Hungary, is responsible for many advances 
in aerodynamics, supersonic and hypersonic airflow characterization, among others. The von Kar- 
man spectrum is a power spectrum of the refractive index fluctuations describing the atmospheric 
turbulence. Thus, phase errors caused by atmospheric turbulence are also standardly modelled to 
match, in average, a von Karman power spectrum for given values of the inner- and outer-scale 
of turbulence. Inner- and outer-scales mean the smallest/largest spatial scales of the fluctuations 
also referred to as smallest/largest eddies in the turbulent medium. 
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FIGURE 2.13 — Plot of the Zernike coefficients representing the low-order contribution from the 
two sub-systems: 1) Telescope+ISLiD+IMaX excluding etalon and PD-plate (also shown in the right 
panel of Figure 2.10), and 2) etalon with the defocus coefficient set to zero(compared with right 
panel of Figure 2.12), independently (left panel). The total added contribution is shown in the 
right panel. 



Phase Diversity plate 

As commented in section §2.3 the PD-plate will introduce a defocus in one of the 
cameras. Apart from that, we were interested in finding out other possible contri- 
butions from the PD-plate to the error budget of the global system. To that aim an 
interferometric calibration of the isolated glass plate was made at the IAC optical 
laboratory 17 by means of a Zygo interferometer. The aberrations were evaluated 
in 9 footprints of the IMaX telecentric beam (1.88 mm of diameter each) over the 
plate, arranged as shown in Figure 2.15 (right panel). A view of the optical setup 
for the calibration is shown in Figure 2.15 (left panel). 

The final results expressed by the Zernike coefficients for each single footprint 
are plotted in Figure 2.16, and the obtained rms-WFE is, in average, 1/21.3 waves. 
Although this is a small contribution, we suspect that it has been over-valuated (i.e. 
the aberration must be still smaller) because of some parasitic structures detected 
in the interferograms probably caused by reflections at the front and back faces 
of the plate. This reflections are the consequence of an inappropriate light source 



The laboratory conditions were: temperature 23°C and humidity 38%. 
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Figure 2.14 — Ripple 
screen randomly chosen 
out of a sample of many 
realizations (von Karman 
power spectrum). The 
figure displays an example 
illustrating the contribution 
of high-order aberrations 
induced by the polishing tool. 



TABLE 2.6 — Thermal model for the PD-plate (May 2008). 

PD INDUCED DEFOCUS 
THERMAL CASE T [°C] Focus displacement [mm] PV [waves] 

Hot 25l) 8A6 1.008 

Operative 24.5 8.51 1.015 

Cold 24.0 8.56 1.019 



wavelength differing from the IMaX working wavelength for which the surface coat- 
ing was designed and manufactured. In consequence, apart from its inherent task 
of displacing the image plane, we will assume a null contribution from the PD-plate 
to the error budget. 

Based on the IMaX thermal study 18 , Table 2.6 lists for three different cases 
within the expected temperature variations range, the amount of defocus (in mm) 
induced by the PD-plate. 

Detector contribution 

The CCD mainly produces three effects that degrade the image quality (Boreman, 
2001), namely: 

• Detector footprint effect: It consist of the integration of the image informa- 

18 Performed in May 2008 by Carmen Pastor Santos at INTA on the basis of the new thermal 
model calculations. 
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FIGURE 2.15 — PD-plate during the interferometric calibration at the IAC. Left: Zygo inter- 
ferometer and PD-plate optical configuration. Right: Cross-section of the PD-plate showing the 
positions and areas (footprints) where the WFEs are evaluated. Courtesy of the IAC optical 
engineer Felix Gracia Temich. 
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FIGURE 2.16 — Zernike coefficients of the PD-plate WFE for 9 footprints over the plate (dif- 
ferent colors). 
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tion over the surface of the detector elements since they have a physical size. 

A detector element with dimensions a x and a y performs a spatial averaging 
of the irradiance falling onto its surface, that in the frequency domain is 
equivalent to a spatial filtering represented by a transfer function OTFf ootpnnt ., 
mathematically expressed as 

U I rf 00 tprint(V x ,Vy) — • , (2.16) 

■KV x a x -KVyCLy 

where v x , u y represent the components of the spatial frequency. 

The smaller the pixel size the broader the resulting transfer function and, 
therefore, the lesser the smoothing in the image when recorded by the CCD. 

• Sampling effect: A sampled-imaging system is not shift-invariant and the 
position of the light reaching the detector, with respect to the pixels, will 
affect the final image. This effect is commonly treated as a statistical average 
of the relative image locations respect to the CCD. 

• Crosstalk effect: Charge-carrier diffusion 19 and charge-transfer inefficiency 20 
over the detector induce a spurious signal on the neighbourhood. 

In what follows, we only model in our simulations the first effect, which is the 
most significant. 

Noise 

Typically, the noise in the image recorded by a CCD has two components: the 
photon noise and the readout noise. Due to the high performance of the CCDs in 
IMaX we can neglect, for the sake of simplicity, the contribution from the readout 
noise. We will simulate the photon noise as having a Gaussian distribution with a 
given rms. As an IMaX requirement and in order to reach a spectral sensitivity of 
10~ 3 decisive to perform precise polarimetric measurements, the value we adopt in 
our simulations is rms-noise= 10~ 3 times the signal in the spectral continuum, i.e. 
a SNR=10 3 in the continuum. 

19 The absorption of photons in a semiconductor material is wavelength dependent: high for 
short-wavelength photons and decreases for longer-wavelength photons. With less absorption the 
long-wave photons penetrate deeper into the material and thus the charges generated should travel 
longer paths to be collected. 

20 lt is caused by incomplete transfer of charge packets along the CCD delay line. 
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2.5.3 Numerical simulations and results 

In order to do a systematical presentation of the numerical experiments we are going 
to describe hereafter for various aberration assumptions, and taking into account 
all the different error sources described above, we consider appropriate to classify 
these error contributions in three main groups, as follows: 

1. Low-order aberrations (LOA) caused by the following optical sub- 
systems: SUNRISE telescope (excluding polishing errors), ISLiD, Interface 
ISLiD-IMaX, and IMaX excluding high-spatial frequency inhomogeneities in 
the etalon but preserving figure errors caused by deformations in the etalon 
and in the windows of its thermally controlled enclosure. 

These aberrations, except for those from the ISLiD-IMaX Interface which will 
not be explicitly considered, are modelled by the set of Zernike coefficients 
listed in column 9 of Table 2.4 with a rms wavefront error rms-LOA = 1/5.2 
waves (Table 2.5). This will be our reference value for LOA in our simu- 
lations. So as to avoid the danger of going to one extreme or the other, 
whereby either excessively pessimistic considerations when evaluating the er- 
ror budget, or the possibility that unexpected error sources could arise (e.g 
uncontrolled spikes in the temperature or pressure variations), we will include 
other cases of aberration in our simulations, with rms-LOA values that vary 
from that of the reference. To do that as simple as possible, we will re-scale 
by a certain factor the list of coefficients in column 9 of Table 2.4. Thus, 
different cases of LOA will be simulated by assuming the following values for 
rms-LOA: 0, 1/12, 1/7, 1/5, 1/4, 1/3 and 1/2 waves. 

The efficiency of the PD-inversion is strongly dependent on the particular set 
of LOA terms included in the simulation. For instance, the problem becomes 
simpler when the main contribution to the aberration comes from very low 
order terms. Then, the inversion code gives good results by assuming a small 
number of unknowns (e.g. 15 or 21 Zernike coefficients or equivalently 15 or 
21 equations). This simplifies the iterative solution of the equations system 
and speed up the convergence. 

2. High-order aberrations (HOA) caused by ripple in the telescope main 
mirror and by high-spatial frequency inhomogeneities in the etalon surface as 
described by laboratory calibration maps. Figure errors in the etalon, induced 



2.5. Testing the robustness of the calibration method 



47 



by high voltages are not included in this budget but in the LOA. 

This HOA enter into the simulation code as respective phase-screens. The 
rms-WFE measured for the etalon in double-pass configuration (rms-etalon) 
is 1/26 waves (Table 2.5) and will be considered as a fixed value in all nu- 
merical experiments in which this effect is included. For the polishing errors 
three different cases will be considered: rms-ripple=0, 2/60 and 2/28 waves. 
The way we write down the precedent values as fractions with a nominator 2 
is meant to denote the effect of reflection of the wavefront in the main mirror 
(e.g. in the case 2/60 we are adding two contributions: 1/60 + 1/60, each 
of them corresponding to the rms-ripple characterizing the mirror surface). 
Apart from the etalon phase-screen, we also have available an amplitude- 
screen mapping its transmission inhomogeneities. 

3. Pixel integration in the CCD and Noise. The first effect is simu- 
lated by an analytical formula (see equation 2.13) as a function of the pixel 
size (fixed value equals 12 /im). The noise contribution is introduced into 
the code as noise-screens generated by an I DL 21 simulator of pseudo-random 
signals with Gaussian distribution for a given rms value. As justified above, 
we take for all considered cases a fixed value rms-noise= 10 -3 x (the signal 
in the spectral continuum). 

In what follows, a variety of aberration models will be made up by combining 
weighted contributions from different error sources. In an attempt to approach 
the real observations when recording bursts of 30 image-pairs for calibrations (see 
section §2.4), 30 realizations of noise will be systematically performed for every 
aberration model adopted. 

Table 2.7 reports on the different ingredients composing the aberration models 
considered in every numerical experiment. The results of these experiments are 
displayed in Figures from 2.17 to 2.28. 

Key notes for the interpretation of Figures 2.17 to 2.28 

Figures 2.17, 2.20, 2.23, 2.26 display PD-speckle reconstructed images from 30 
realizations for every aberration assumption. Lying below every restored scene we 
present one out of the 30 realizations of the in-focus degraded images. To facilitate 

21 The Interactive Data Language is a widely known programming language that is a popular 
data analysis language among scientists. 
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the visual comparison between different images, all of them are displayed by em- 
ploying a common gray-scale ranging from the minimum to the maximum intensity 
values of the true object shown at the upper right corner for reference. This way, 
real differences in contrast and resolution can easily be detected by eye. 

Apart from the visual comparison, we also calculate the values of different 
parameters to quantify the degree of success achieved in our numerical experiments 
basically consisting of: degradation followed by PD-inversion. Figures 2.18, 2.19, 
2.21, 2.22, 2.24, 2.25, 2.27 and 2.28 are arranged in three columns. The panels in 
left, middle and right columns of the figures correspond to different assumptions for 
the rms-ripple: 0, 2/60 and 2/28 waves, respectively. The abscissae in all panels 
stand for different cases of rms-LOA. Bullets linked by dotted lines show the mean 
values resulting from PD-inversions. Vertical bars represent the dispersion of these 
results caused by 30 different realizations of photon noise. Dashed lines represent 
the values of the different parameters in the degraded (smeared) images simulated. 

To evaluate the reliability of the results from the PD-process we need some 
expected reference values for comparison, e.g. the contrast of the true object that 
should ideally by reached after restoration, or the zero value that should ideally re- 
sult from the subtraction of certain wavefronts (see thick solid lines in Figures 2.18, 
2.21, 2.24 and 2.27). We also use as a reference the results of the restoration iz 
performed with expression 1.21, by only employing the input Zernike terms simu- 
lating the LOA of the system, 4>z (see thin solid lines in Figures 2.18, 2.21, 2.24 
and 2.27). 

Subscripts t,s,r and Z in Figures 2.18, 2.21, 2.24 and 2.27 stand for true ob- 
ject, smearing (simulated degradation), reconstruction with PD, and restoration 
performed with only the Zernikes representing the LOA in the input, respectively. 

Experiment 1 

Combined effects considered (see also Table 2.7). 

• rms-ripple = 0, 2/60 and 2/28 waves. 

• rms-LOA = 0, 1/12, 1/7, 1/5, 1/4, 1/3 and 1/2 waves. 

• rms-noise = l(P 3 x continuum signal (30 realizations per aberration model.) 

• No CCD effects. 

• No etalon effects: \H(p,8)\ =1 , e (p,0) =0. 

• PD-defocus in mm (degradation/inversion): 8.51 / 8.51. 

• Inversion with 25 coefficients. 
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TABLE 2.7 — Error contributions for diverse numerical simulations 



ERROR CONTRIBUTION EXPERIMENT 



SOURCE 


VALUE 


1 


2 


3 


4 


rms-ripple 


n 2 2 

°< 60' 28 waves 


V 


V 


V 


V 


rms-LOA 


n l lllll 

0, 12 , 7 , 5 , 4 , 3 , 2 waves 


V 


V 


V 


V 


rms- noise 

rms-etalon 


1CP 3 X continuum signal 

(30 realizations/aberr. model) 

1 

2g waves 


V 


V 
V 


V 
V 


V 
V 


Etalon amplitude 


\H(p,8)\^l 




V 


V 


V 


CCD 


12 /im/pix — > 0.055 "/pi x 






V 


V 


PD-defocus 


DEGRADATION / INVERSION 
8.51 mm (PV 1.00A) / 8.51 mm 
9.00 mm (PV 1.06A) / 8.51 mm 


v/ 









A glance at Figure 2.18, rows (a), (b), (c), reveals that the results derived from 
the application of PD to the simulated observations and the expected values for 
different parameters (thick solid line) are in good agreement when rms-ripple=0 
or 2/60 waves (left and middle columns) and rms-LOA=0, 1/12, 1/7, 1/5 or 1/4 
waves. Thus the retrieved rms contrast for i r ranges from 16% to 17.4% versus 
the reference rms(^)=17.5%; the mean absolute difference mean(|i t — i r \) remains 
always below 2.5% and rms(^ — i r ) below 3.5%. In many cases the results are also 
reasonably good for rms-L 0/4=1/3 waves but always noticeably worse for the case 
rms-L 0/4=1/2 waves. 

The figure shows that the results from the PD-inversions deviate progressively 
from the reference values as long as the amount of ripple and LOA augment in the 
simulated aberrations. This behavior reflects the limitations of the method stem- 
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ming from: 1) the unsensed aberrations due to the limited number of Zernikes used 
in the inversion (25 coefficients versus 36 employed to simulate LOA in the total 
system -see Table 2.4; 2) the influence of the also non-sensed HOA terms present in 
the ripple 22 ; 3) the noise in the data; and 4) the instabilities in the PD-method itself. 

The comparison of the bullet-dotted with the dashed lines evidences the enor- 
mous improvement obtained from PD-inversions with respect to the smeared im- 
ages. It is worth mentioning that Figure 2.18, rows (a), (b), (c), (d), show the good 
coincidence (within the above mentioned range rms-ripple < 2/60 and rms-LOA < 
1/4) between the bullet-dotted and the thin solid lines, the latter representing the 
results of the deconvolution performed by only employing the Zernike coefficients 
assumed as the input to simulate LOA in the total system. This coincidence, even 
extensive to the third column, reveals that the PD-algorithm compensate reason- 
ably well at least for the degradation caused in the simulation by LOA terms. Note 
also that, for a given parameter and for rms-LOA < 1/4, the variation trend in the 
bullet-dotted and thin solid lines is quite similar in the three columns, being the 
main difference between these variations a vertical offset which increases from the 
left to the right columns. We ascribe these offsets to a loss of contrast caused by 
the ripple. 

Ripple errors are high-frequency phase errors and can be viewed as noise in 
the wavefront since they can hardly be modelled even with very high order Zernike 
terms in a polynomial expansion. The consequence of these errors is a transfer of 
energy from the peak of the PSF to its far wings thus producing stray-light and 
consequently an unrecoverable loss of contrast as evidenced in Figure 2.18, row 
(a) - this degradation increases as the ripple-scale decreases (Bello, 2001). Thus, 
large values of rms-ripple have more impact on the image contrast than on the 
image resolution since the latter is related to the FWHM of the PSF that, for well 
corrected systems, changes rather slowly 23 . Note that the Strehl ratio is strongly 

22 ln further experiments also the inhomogeneities in the etalon phase will bring additional HOA 
contribution. 

23 According to equations 1.15, 1.16 and the auto-correlation theorem, the PSF of an optical 
system is proportional to the power spectrum of the generalized pupil function. Consequently, 
the presence of high-spatial frequencies in (f>(p,8) will cause a long tail in the power spectrum or 
equivalently long wings in the PSF of the system. Furthermore, the energy conservation forces 
the normalization of the volume enclosed by the PSF; thus the extension of its wings will be done 
at the expenses of an intensity decrease at shorter distances, in particular at its central point. 
If, apart from the ripple, the optical system is well corrected, this transfer of energy from short 
to long distances will provoke a clear decrease in the central intensity of the PSF as long as the 
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dependent on the amount of ripple (see panels in row (f) of Figure 2.18). The 
case rms-ripple=2/28 has only an academic interest in our simulations because we 
expect a better polishing quality in the SUNRISE main mirror surface 24 . 

Stray-light in solar photometry represents a problem that classically has been 
faced from various approaches, all of them based on modelling the far wings of the 
PSF by linear combinations of Gaussian and Lorentzian functions, the parameters 
of which are determined from solar aureole measurements (Martinez Pillet, 1992, 
and references therein). In spectro-polarimetry, being the case of IMaX, stray-light 
is a particularly delicate issue because it represents a contamination of the mag- 
netic signal at a given feature by the signal originated in other structures located 
all around the feature even at quite large distances (e.g. ~ arc minutes). 

The loss of contrast and spatial resolution, as long as the rms-ripple and 
rms-LOA increase, can be visually appreciated in the strips of restored images 
in Figure 2.17. Image resolution is particularly affected by rms-LOA > 1/4 waves: 
medium-high spatial frequency structures, grains-like or fringes-like patterns start 
arising. These artifacts originate because the OTF in these cases drops to close-to- 
zero values in a range at medium-high spatial frequencies 25 and are more evident in 
single-restored images before doing the speckle summation of equation 1.21 (see 
discussion at the end of section §1.4). Thus, the final product for single restorations 
is often arbitrarily affected by over restoration of some particular spectral compo- 
nents. 

In Figure 2.19, thick lines represent the variation of the Zernike coefficients (up 
to index 11) used as the input to simulate the LOA (i.e. the reference for compar- 
ison), whereas thin lines mark simply the zero level. Bullet-dotted lines represent 
again values retrieved from the PD-inversion process. In most cases the coincidence 
of the retrieved coefficients with the expected values is quite good even in the right 

ripple increases, but because of the very steep slopes in the PSF core, the FWHM of the PSF will 
scarcely be affected. 

24 ln fact the expected rms values for polishing errors as measured on the mirror surface range 
from 5 to 10 nm. For our working wavelength, A525.02 nm, this range corresponds to 1/100 
to 1/50 waves, which includes the value 1/60 assumed as the second case for rms-ripple in our 
simulations. 

2j This effect is favoured by the central obscuration in the SUNRISE'S pupil inherent to the 
particular design adopted for this telescope. The central obscuration depresses the central region of 
the MTF in comparison with its typical shape for diffraction-limited telescopes with clear aperture. 
See Figure 2.29 illustrating this depression. 
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column of the figure (we exclude from this comparison the case rms-L OA= 1/2 
waves). Nevertheless, the deviations encountered from the reference values could 
probably be ascribed to cross-talk between the limited number of Zernike terms 
employed in the PD-inversions (25 terms) in an attempt to compensate for the 
contribution of higher-order terms which are in fact present in the simulated wave- 
front: we have used 36 terms to simulate the system LOA and, furthermore, we 
can not disregard the contribution from polishing errors and noise to the signal. 
In other words, as argued by Lofdahl & Scharmer (1994) the inversion program 
compensates for the missing high-order terms by introducing errors in the lower 
Zernike coefficients. To understand such compensations, we note that the inver- 
sion algorithm does not attempt to derive the best fit to the wavefront. Because 
of the chosen error metric (see equation 1.19), the inversion attempts to find two 
internally consistent (i.e. derived from the same Zernike coefficients) transfer func- 
tions for the focus and defocus images such that the error metric is minimized. The 
value of the transfer function at a specific frequency is given by an average over the 
aperture of all differences in phase taken at a fixed separation (see equation 1.15). 
Thus, mainly at low and intermediate frequencies, different wavefronts can give 
similar transfer functions. This gives the possibilities for compensations and allows 
the inversion procedure to find the wavefront that gives the best fit to the transfer 
functions given the observed data and the number of Zernike coefficients used to 
represent the wavefront. Nevertheless the improvement given by such compen- 
sations is limited because of the absence of high-spatial frequencies in the lower 
aberrations. 



The panels in rows (d) and (e) in Figure 2.18 deserve a particular comment. 
The inversion process yields rms(</v) values that fit rms(ijl s ) very well, at least up to 
rms-LOA=l/4 waves. In spite of this coincidence, the deviation of the parameter 
rms((/> s — 4> r ) from zero increases as the ripple and the LOA augment. This variation 
preserves a similar trend in the three columns of row (e) although it is affected 
by a vertical offset in the middle and right columns. Note that the parameter 
rms(0 s — 4>z) is also affected by similar offsets. In fact these offsets coincide, as 
expected, with the rms-ripple value assigned to each case. The slope in the variation 
of rinses — 4> r ) could be justified in terms of the compensations mentioned in the 
previous paragraph (cross-talk between low order modes), which reproduce good 
transfer functions but not so good wavefront shapes. Apparently, the difference in 
the wavefront shape does not imply a substantial difference between rms((/> s ) and 
rms(0 r ) as shown in the panels of row (d). 
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Experiment 1 
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FIGURE 2.17 — Experiment 1: PD reconstruction for different amounts of aberration. Seven 
different contributions from low-order aberrations are adopted: rms-LOA = 0, A/12, A/7, A/5, 
A/4, A/3, A/2, and for each of them three possible cases of polishing errors are considered: rms- 
ripple = 0, 2A/60, 2A/28. The exam of the PD-speckle reconstructed images from 30 realizations 
allows us to check the robustness of the IMaX WFE calibration. The IMaX case is emphasized in 
yellow and the so-called true object representing the free-from-aberrations and -diffraction image 
is located in the upper right corner. All the images are displayed by employing a common gray 
scale ranging from the minimum to the maximum intensity values of the true object. Specific 
differences characterizing every Experiment are summarized at the bottom of the figure. 



54 



CHAPTER 2. In-flight calibration of IMaX aberrations 



EXPERIMENT 1 / Inversions evaluation 




C-.D D.1 {1.2 OS OA C-.E 



3fl 



2 



ID 



10 




d.d- n.i aj ou D-i 




D-.1 D.2 1 j (.4 U 




4,A t,1 U V Ai4 A,9 




-I 1 — J — ± L LJ 

fl.ii 1,1 U V U V 
m-LCw, [«mi) 




im-LCw, {vhwu) 




■0,0 fcl M CJ D-.+ D,D 0,0- D.1 US iJ SU DA Qui D-.1 D.Z flj «U ** 

imp-lAh (wh) ren-LAl C^n^wJ rmr-LO* 



FIGURE 2.18 — Comparative plot for the numerical Experiment 1. To understand the meaning of 
the different symbols and labels see section §2.5.3: Key notes for the interpretation of Figures 2.17 
to 2.28. 
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EXPERIMENT 1 / Inversions evaluation (cont) 
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EXPERIMENT 1 / Zernike coefficients evaluation 
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FIGURE 2.19 — Comparative plot of Zernike coefficients for the numerical Experiment 1. Bullets 
linked by dotted lines represent the variation of the mean Zernike coefficients retrieved from the 
PD-inversions of 30 realizations. Thick solid lines represent the variation of the original coefficients 
used as code inputs to simulate the degradation caused by low-order aberration (LOA) terms. Only 
the coefficients indexed from 4 to 11 are included in the figure. 
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EXPERIMENT 1 / Zernike coefficients evaluation (cont) 
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Experiment 2 

Combined effects considered (see also Table 2.7). 

• rms-ripple = 0, 2/60 and 2/28 waves. 

• rms-LOA = 0, 1/12, 1/7, 1/5, 1/4, 1/3 and 1/2 waves. 

• rms-noise = l(P 3 x continuum signal (30 realizations per aberration model.) 

• No CCD effects. 

• Etalon effects: \H(p,0)\ and <f) e (p,9) derived from lab. calibration maps. 

• PD-defocus in mm (degradation/inversion): 8.51 / 8.51. 

• Inversion with 25 coefficients. 



The only difference of this Experiment with respect to Experiment 1 is that 
now we include in the simulated degradation the contribution from \H{p,9)\ and 
4> r {p,d) mapping the transmission inhomogeneities and the HOA of the etalon, 
respectively. Note that the PD-code is only designed for the retrieval of the phase- 
error function in the pupil but not of the pupil transmission function. Thus, in 
the inversion process the function \H(p,9)\ cannot be considered as an unknown 
in any case but one has to presume a certain model for it. In a first attempt we 
used for the inversion the same amplitude-screen as the one applied to simulate 
the degradation but this way of proceeding generated a disastrous behaviour in the 
PD-code with strong fluctuations during the iterative process that minimizes the 
error metric, and eventually collapsed the calculations. So, we were compelled to 
use in the inversion the simple ideal model \H(p,9)\=l all over the pupil. 

Now we center our attention in the impact of these new complications in the 
efficiency of the PD-code to recover the global phase-error introduced as the input 
to simulate the degradation. Firstly, we concentrate in Figure 2.21, rows (a), (b), 
(c), for rms-ripple=0, 2/60 waves and for rms-LOA < 1/4 waves. As in Experiment 
1, a good agreement between the results of the inversion and the reference values 
is found. In most cases, the contrast rms(i r ) ranges from 16.5% to 18.4%, that 
is, close to the reference value 17.5%. The trends and ranges in the variations of 
the mean absolute differences mean(|z t — i r \) and of the rms(?t — i r ) remain quite 
similar to those in Figure 2.18. However, a significant difference with Experiment 1 
deserves special discussion. Now, for the cases with small rms-LOA < 1/7 waves, 
the inversion yields image contrasts rms(i r ) overvalued with respect to the reference 
one. 

We ascribe this result to the non-realistic assumption made when setting \H(p, 9)\=1 
in the inversions. The error metric to be minimized in the PD-inversions does not 
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FIGURE 2.20 — Experiment 2: PD reconstruction for different amounts of aberration. Seven 
different contributions from low-order aberrations are adopted: rms-LOA = 0, A/12, A/7, A/5, 
A/4, A/3, A/2, and for each of them three possible cases of polishing errors are considered: rms- 
ripple = 0, 2A/60, 2A/28. The exam of the PD-speckle reconstructed images from 30 realizations 
allows us to check the robustness of the IMaX WFE calibration. The IMaX case is emphasized in 
yellow and the so-called true object representing the free-from-aberrations and -diffraction image 
is located in the upper right corner. All the images are displayed by employing a common gray 
scale ranging from the minimum to the maximum intensity values of the true object. Specific 
differences characterizing every Experiment are summarized at the bottom of the figure. 
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EXPERIMENT 2 / Inversions evaluation 
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FIGURE 2.21 — Comparative plot for the numerical Experiment 2. To understand the meaning of 
the different symbols and labels see section §2.5.3: Key notes for the interpretation of Figures 2.17 
to 2.28. 
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EXPERIMENT 2 / Inversions evaluation (cont) 
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EXPERIMENT 2 / Zernike coefficients evaluation 
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FIGURE 2.22 — Comparative plot of Zernike coefficients for the numerical Experiment 2. Bullets 
linked by dotted lines represent the variation of the mean Zernike coefficients retrieved from the 
PD-inversions of 30 realizations. Thick solid lines represent the variation of the original coefficients 
used as code inputs to simulate the degradation caused by low-order aberration (LOA) terms. Only 
the coefficients indexed from 4 to 11 are included in the figure. 
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explicitly contain neither the phase, <fi, nor the amplitude, \H\, that define the gen- 
eralized pupil function (expression 1.16), but the OTFs that can be derived from 
the latter. Thus, the defective information assigned to the amplitude \H\ can be 
compensated by some extra errors - in comparison with Experiment 1 - in the 
phase functions provided that the resulting generalized pupil functions still yield 
appropriate OTFs that minimize the error metric. These extra errors are evident 
from the comparison of Figures 2.18 and 2.21, rows (d) and (e), and can be the 
cause of the over-restoration mentioned above. Figure 2.21, row (e), also shows 
the deviations of rms(4> s — <pz) hom zero. These vertical offsets are related to the 
HOA contribution from polishing errors in the main mirror and inhomogeneities in 
the etalon surface and, as expected, in the three panels they very closely amount 
to RSS=\/ rms-ripple 2 + rms-etalon 2 . 

Experiment 3 

Combined effects considered (see also Table 2.7). 

• rms-ripple = 0, 2/60 and 2/28 waves. 

• rms-LOA = 0, 1/12, 1/7, 1/5, 1/4, 1/3 and 1/2 waves. 

• rms-noise = 10 X continuum signal (30 realizations per aberration model.) 

• CCD effects. 

• Etalon effects: \H(p,6)\ and <j) e (p,8) derived from lab. calibration maps. 

• PD-defocus in mm (degradation/inversion): 8.51 / 8.51. 

• Inversion with 25 coefficients. 



In this Experiment we add to the image degradation one more contribution with 
respect to Experiment 2. It consist of the integration over the pixel area, of the 
information reaching the image plane as commented in section §2.5.2. Table 2.7 
ticks off the error contributions taken into account for this simulation. Dotted line 
in Figure 2.29 shows the OTFf ootpr ; nt representing the degradation effect which is 
incorporated to the total error budget by simple multiplication by the OTF describ- 
ing the global aberration in Experiment 2. The new ingredient in the simulated 
aberration has no significant impact in the results of the PD-inversion. Thus, for 
rms-ripple < 2/60 and rms-LOA < 1/4 waves the differential impact in the pa- 
rameters rms(v), mean(\i t — i r \) and rms(it — i r ) - Figure 2.24, rows (a), (b), 
(c) -, with respect to Experiment 2, is smaller than 2 tenths %. Similar coinci- 
dences between the results from Experiments 2 y 3 are achieved for the parameter 
rms(0 r ). However, the influence of the CCD pixel integration is remarkable in the 
Strehl ratio as evidenced by comparing Figures 2.21 and 2.24, rows (f). 
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FIGURE 2.23 — Experiment 3: PD reconstruction for different amounts of aberration. Seven 
different contributions from low-order aberrations are adopted: rms-LOA = 0, A/12, A/7, A/5, 
A/4, A/3, A/2, and for each of them three possible cases of polishing errors are considered: rms- 
ripple = 0, 2A/60, 2A/28. The exam of the PD-speckle reconstructed images from 30 realizations 
allows us to check the robustness of the IMaX WFE calibration. The IMaX case is emphasized in 
yellow and the so-called true object representing the free-from-aberrations and -diffraction image 
is located in the upper right corner. All the images are displayed by employing a common gray 
scale ranging from the minimum to the maximum intensity values of the true object. Specific 
differences characterizing every Experiment are summarized at the bottom of the figure. 
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FIGURE 2.24 — Comparative plot for the numerical Experiment 3. To understand the meaning of 
the different symbols and labels see section §2.5.3: Key notes for the interpretation of Figures 2.17 
to 2.28. 
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FIGURE 2.24 — (cont). 
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FIGURE 2.25 — Comparative plot of Zernike coefficients for the numerical Experiment 3. Bullets 
linked by dotted lines represent the variation of the mean Zernike coefficients retrieved from the 
PD-inversions of 30 realizations. Thick solid lines represent the variation of the original coefficients 
used as code inputs to simulate the degradation caused by low-order aberration (LOA) terms. Only 
the coefficients indexed from 4 to 11 are included in the figure. 
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Concerning the retrieved Zernike coefficients from the PD-process, they look 
very similar in Figures 2.22 and 2.25, within the range of aberrations mentioned 
above (rms-ripple < 2/60 and rms-LOA < 1/4 waves). Nevertheless, as discussed 
in previous paragraphs, small differences in the Zernike coefficients of both Exper- 
iments would not be crucial provided that the resulting generalized pupil functions 
yield appropriate OTFs to obtain a good minimization of the error metric. This 
seems to be the case in the present Experiment to judge by the good image restora- 
tions achieved. Therefore, we will consider the results from Experiment 3 as the 
reference for the next numerical test. 

Experiment 4 

Combined effects considered (see also Table 2.7). 

• rms-ripple = 0, 2/60 and 2/28 waves. 

• rms-LOA = 0, 1/12, 1/7, 1/5, 1/4, 1/3 and 1/2 waves. 

• rms-noise = l(P 3 x continuum signal (30 realizations per aberration model.) 

• CCD effects. 

• Etalon effects: \H(p,6)\ and <f) e (p,8) derived from lab. calibration maps. 

• PD-defocus in mm (degradation/inversion): 8.51 / 9.00. 

• Inversion with 25 coefficients. 



The recent IMaX thermal design (May 2008) predicts for the PD-plate, dur- 
ing the flight, a temperature variation range of 24.5 ± 0.5 °C. From this thermal 
model, simulations with CODE V estimate the defocus displacement induced by 
the PD-plate ranging from 8.46 mm (PV-WFE=1.008 waves; hot case) to 8.56 
mm (PV-WFE=1.019 waves; cold case), being 8.51 mm (PV-WFE=1.015 waves) 
the displacement corresponding to the operative case (24.5 °C). 

The amount of displacement of the image plane produced by the PD-plate, also 
called diversity, is an input parameter in the PD-code. One of our concerns has 
always been how critical the value assigned to the diversity is, or in other words, 
how much robust is the PD-code in rendering good results versus possible errors 
in the value of the diversity. This concern arises from the possibility of having 
during the flight unexpected peaks in the temperature fluctuation, out of the range 
predicted by the thermal model. 

To clarify this question, we have designed the following numerical experiment: 
in the inversion process we will take for the diversity the expected nominal operative 
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Experiment 4 
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FIGURE 2.26 — Experiment 4: PD reconstruction for different amounts of aberration. Seven 
different contributions from low-order aberrations are adopted: rms-LOA = 0, A/12, A/7, A/5, 
A/4, A/3, A/2, and for each of them three possible cases of polishing errors are considered: rms- 
ripple = 0, 2A/60, 2A/28. The exam of the PD-speckle reconstructed images from 30 realizations 
allows us to check the robustness of the IMaX WFE calibration. The IMaX case is emphasized in 
yellow and the so-called true object representing the free-from-aberrations and -diffraction image 
is located in the upper right corner. All the images are displayed by employing a common gray 
scale ranging from the minimum to the maximum intensity values of the true object. Specific 
differences characterizing every Experiment are summarized at the bottom of the figure. 
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FIGURE 2.27 — Comparative plot for the numerical Experiment 4. To understand the meaning of 
the different symbols and labels see section §2.5.3: Key notes for the interpretation of Figures 2.17 
to 2.28. 
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FIGURE 2.28 — Comparative plot of Zernike coefficients for the numerical Experiment 4. Bullets 
linked by dotted lines represent the variation of the mean Zernike coefficients retrieved from the 
PD-inversions of 30 realizations. Thick solid lines represent the variation of the original coefficients 
used as code inputs to simulate the degradation caused by low-order aberration (LOA) terms. Only 
the coefficients indexed from 4 to 11 are included in the figure. 
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value 8.51 mm, whereas in the degradation process, we will assign to the diversity 
a different value thus simulating an unexpected temperature peak which we are not 
aware of. 

In the first realization of the Experiment we have selected for the degradation 
process a rather large diversity value (12 mm) as compared with the nominal one 
used in the inversion (8.51 mm). Always within our standard range of analysis (i.e. 
for rms-ripple < 2/60 and rms-LOA < 1/4 waves), the inversion produces over- 
restoration: rms(v) ranging from 17.5 to 20.7 %, mean(|i t — i r \) from 2.6 to 5.5 % 
and rms(ij — i r ) from 3.6 to 7.1 %. Also the retrieved Zernike coefficients, partic- 
ularly «4, differ substantially from the reference values. A curious result is that the 
values obtained for are larger than those corresponding to the Experiment 3 by a 
roughly constant amount of about 0.1 waves, corresponding to a defocus distance 
of ~ 3 mm, similar to the diversity error (3.49 mm). It seems like if the error in the 
diversity would have misled the PD-algorithm in such a way that this error is now 
ascribed to the image in nominal focus position, in the form of defocus aberration. 
This could explain the resulting over- restoration mentioned above. Since this is an 
extreme experiment considering diversity values far away from the IMaX case, we 
do not display the results. 

In a second and more realistic numerical experiment we have selected a diversity 
of 9 mm for the degradation process whereas we have kept the nominal value for the 
inversion. Thus, the error introduced, 9.00 - 8.51 = 0.49 mm, is greater by a factor 
of 10 than the expected peak variation of the diversity (± 0.05 mm) according to 
the IMaX thermal model (Table 2.6). Even so, the PD-inversion gives very similar 
values to those in Experiment 3 (compare rows (a) to (f) in Figures 2.24 and 2.27). 
However, small discrepancies can be detected in the retrieved Zernike coefficients, 
from the comparison of Figures 2.25 and 2.28. In particular, the bullet-dotted line 
for Q!4 in the latter figure is shifted upwards by ~ 0.015 waves (~ 0.44 mm of defo- 
cus distance) with respect to the corresponding line in the former figure. Again we 
find that the error in the knowledge of the real diversity is transferred to the budget 
of real aberrations in the system, in the form of a defocus error. Nevertheless, in 
this case, the impact on the restorations can be consider as negligible since we get 
very similar results to those obtained in Experiment 3. 

The repercussion of the discrepancies between the results of Experiments 1, 2, 
3 and 4, are very subtle as to be easily detected by comparing the Figures 2.17, 
2.20, 2.23 and 2.26. 
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[arcsec" 1 ] 

FIGURE 2.29 — MTF shape for the case of IMaX. Left panel: Diffraction limited MFT (solid 
line), MTFfo otpr ,„t (clotted line) and MTF resulting from the total error contribution (dashed line). 
Right panel: 3D representation of the latter. 

Case of IMaX 

According to the aberrations estimated from laboratory measurements and realistic 
models, the case of IMaX fits into Experiment 3 for rms-ripple < 2/60 and rms- 
LOA = 1/5 waves. This implies a loss of contrast of ~ 5 % of the real value and 
rms(«t — i r ) ~ 3 % (see Figure 2.24). Figure 2.29 {left panel) shows the diffrac- 
tion limited MTF {solid line) for SUNRISE. The dotted line is the MTFf ootp „ n t 
representing the pixel integration effect in the CCD. The dashed line is a section 
of the MTF resulting from the total error contribution estimated for the IMaX im- 
ages including also the CCD pixel integration effect. The central depression in the 
MTF profile is very sensitive to the significance of LOA in such a way that it can 
almost reach a zero level at some particular spatial frequencies (see right panel in 
Figure 2.29). These singularities complicate the construction of the noise filter (see 
section §1.3.2) to be applied in the restoration process. 

Figure 2.30 is a comparative plot of the Zernike coefficients used for degrading 
the images (from index 4 to 36) and those retrieved from the inversion process 
(from index 4 to 25). 

2.6 Conclusions 

In this chapter, a method for in-flight calibration of aberrations in IMaX has been 
proposed. The method is based on the Phase Diversity (PD) inversion procedure. 
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IMaX case 
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FIGURE 2.30 — Zernike coefficients for the IMaX case. The figure shows for comparison the 
coefficients used to degrade the image (input) and those retrieved from the PD-inversion (output). 



The robustness of the method has been tested by numerical experiments sim- 
ulating different aberration components. The sources of aberration have been 
modelled and added in subsequent experiments. The repercussion of every new 
added ingredient in the final result from the inversion has been evaluated. 

The PD-code does not accurately reproduce the shape of the wavefront errors 
but provides reliable OTFs for subsequent satisfactory restorations. 

Small-scale irregularities in the wavefront error are not detectable by the PD- 
code where a finite (rather low) number of Zernike terms is used. This limitation 
mainly produces stray-light all over the restored image and consequently a loss of 
contrast. This effect is tolerable within certain margins, and fixes constraints to 
the polishing quality in the main mirror of the SUNRISE'S telescope and to the 
inhomogeneities in the etalon of IMaX. 

The inhomogeneities in the transmission function of the etalon are converted 
by the inversion process into some extra errors in the resulting wavefront which 
partially compensate for the loss of contrast caused by the high-order aberrations 
terms non-sensed by the PD-code. 

Experiment 3 (see Table 2.7) includes all the error sources assumed as signifi- 
cant for IMaX, and the results of its inversions can be considered reliable enough 
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as to validate the method proposed for the calibration of errors in the images of 
IMaX, at least for aberrations lying in the range: rms-ripple < 2/60 and rms-LOA 
< 1/4 waves. The case of IMaX, with the expected values of rms-ripple < 2/60 and 
rms-LOA ~ 1/5 waves, is included in this range of aberrations and, consequently, 
the calibration method is applicable to IMaX. 

For the case of IMaX, the residual errors in the proposed calibration method 
induce, in turn, errors in the subsequent restoration: mean(|ij — i r \) < 2.5 %, 
rms(«t — i r ) < 3.3 %, and a loss of rms contrast < 5 %. 

The amount of defocus (also named diversity) produced by the PD-plate is a 
critical parameter in order for the PD-algorithm to work optimally. We have got 
acceptable results for errors in the diversity of ~ 5 %. 



Part II 

Study of proper motions in solar 
active regions 
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3 

Solar active regions 

The second part of this thesis, devoted to the study of solar active regions, 
is introduced in this chapter based on the revision of literature published on 
this topic. Information about solar active regions is firstly presented starting from 
the historical point of view, secondly going into the theoretical and observational 
framework and finally reviewing the state-of-the-art in the research on the dynamics 
and evolution of photospheric structures around sunspots. 

3.1 Introduction 

The importance of the Sun for humankind is more than justified as well as the 
research on solar physics in which the Sun appears to be a fabulous laboratory 
to study multitude of phenomena and fundamental processes in astrophysics. By 
the analysis of the radiation (photons) coming from the different layers of the 
Sun (see Figure 3.1) scientists have been able to figure out a large amount of 
information about the underlying physics behind our closest star. Telescopes and 
solar instrumentation have played an important role since the very first ground- 
based observations up to the more recent space telescopes and future missions 
trying to untangle the behavior of the Sun. The Sun as a variable star experiences 
a changing magnetic activity framed in what is called the solar cycle. The main 
cycle is composed by two consecutive 11 years cycles and represents a reversal of the 
global magnetic field. The solar cycle is evidenced by the variation of the number 
of sunspots, number of solar flares and solar irradiance among other parameters. 
Magnetism then dominates the Sun's behavior and can be studied in more detail 
than in any other stellar body. The vast majority of light leaving the Sun and 
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reaching us (approx. 8 minutes later) comes from the photosphere, a thin layer 
of about 500 km of thickness (sort of 0.04% of the solar radius) that is actually 
the densest part of the solar atmosphere, composed by convective cells known as 
solar granulation. Underneath the photosphere, convection is taking place carrying 
the energy released by thermonuclear fusion in the solar core. Plasma motions in 
the underlying convection zone are therefore driving the dynamical behavior of the 
photosphere and reveal through the granulation pattern. 

FIGURE 3.1 — Cartoon repre- 
senting different layers in the 
Sun: Photosphere (Ph, yellow) , 
Chromosphere (Ch, orange) and 
Corona (Cr, blue). The plot 
shows the stratification of temper- 
ature (black line) and mass den- 
sity (dashed line) when moving ra- 
dially outwards the Sun (height). 
The temperature minimum occurs 
at a height of 500 km above the 
deepest layer observable with op- 
tical telescopes. The temperature 
minimum marks the upper limit of 
the photosphere. The Transition 
Region (red) is characterized by a 
huge increment of temperature in 
just a few hundredth kilometers. 



3.2 Sunspots 

Sunspots are the most evident manifestation of the solar activity and of the solar 
cycle. They are related to the majority of violent phenomena occurring in the solar 
external layers and moreover a clear sign of the magnetism effects in the dynamics 
and structure of our closest star. Solar active regions (sunspots and plages) are 
labeled for identification by consecutive numbers as they appear on the visible solar 
disk. One of the most widely used identifiers is the NOAA (National Oceanic and 
Atmospheric Administration) catalogue. Many important facts about the sunspots 
have been studied for many years and we have learned things about its nature, 
origin and evolution, but there are still fundamental questions we are willing to 
clarify. One of the important issues about sunspots was to discover that they were 
the first astronomical object known to harbor a magnetic field (Hale, 1908). 
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Sunspots are nowadays interpreted as complex structures having strong mag- 
netic fields that inhibit the plasma convection reducing therefore the efficiency in 
the convective energy transport. The temperature of sunspots is for this reason 
lower than that in the surrounding photosphere. They are seen as black entities 
when comparing it with the outside brightness but actually the temperature inside 
sunspots goes up to 5000 K. 

Study of sunspots is important because they are connected to intense and explo- 
sive phenomena in the chromosphere and solar corona. Coronal loops, for instance, 
are magnetic lines anchored to sunspots that after mutual interaction and recon- 
nection liberate strong amounts of energy in the form of the so-called flares and 
coronal mass ejections (CME). The released energy travels millions of kilometers 
before interacting with the Earth's atmosphere and affecting the terrestrial com- 
munications and orbiting satellites. Clarifying the sunspots internal structure and 
their physical processes is demanding to understand and predict the solar activity 
and its influence on Earth. 



3.3 Historical overview 

The first observations of sunspots on the solar surface brought another view re- 
garding the incorruptibility of heavens so that at the beginning many astronomers 
rejected to believe that sunspots were actually attached to the Sun. J. Fabricius 
was the first to write a book on sunspots in 1611 (Fabricius, MDCXI) and indepen- 
dently C. Scheiner claimed to observe sunspots the same year associating them to 
orbiting objects (Scheiner, 1623-30). Scheiner is also regarded as the first who evi- 
dences velocity fields in the solar photosphere noting that sunspots near the equator 
traverse the solar disc faster than those at higher latitudes. Galileo associated the 
sunspots to a sort of thick clouds sharing the solar rotation, as dense that they were 
able to block the Sun's light. He was the first that used the terminology of solar 
latitude and longitude to locate sunspots on the Sun. Galileo's accurate observa- 
tions using a rather good telescope for the time being, allowed him to look at many 
details on the solar surface. He was also the first to discover brighter zones which 
he called faculae. But it was not until Keppler that sunspots were indeed put right 
on the solar surface (Keppler, 1955). Keppler had first postulated that the Sun 
should rotate in less than three months but later on reduced the rotating period 
to just one day. After proving by means of geometrical arguments that Scheiner 
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was wrong when considering the behavior of sunspots at the solar limb, Keppler 
settled the question, being convinced that spots are on the Sun as we all know now. 

Sunspots observations continued in the XVII century although there is a lapse 
from 1645 to 1715 where a very few sunspots were observed. This fact is known 
as the Maunder minimum and corresponds to a time period in which the Earth 
suffered from very low temperatures and cold winters that are associated to the 
reduction of the solar activity. Nowadays, the Sun is still affecting us directly and 
indirectly and the Sun's changing magnetic behavior has an influence on the solar 
neighborhood climate. 

More solar properties and parameters, as the solar axis inclination and differ- 
ential rotation rate were inferred by Carrington (1859) and confirmed by several 
studies afterwards. The dynamics of the plasma in the solar photosphere was found 
to show several convective patterns at different spatial scales. In addition to the 
granular pattern (with typical diameters of 150-2500 km), the mesogranular one at 
intermediate scales (5-10 Mm) and the supergranular one at large scales (over 20 
Mm) organize the plasma in a sort of hierarchical convective structure. There is 
still an open debate concerning the nature, origin and structure of these convective 
patterns. 

For a more extensive and detailed historical review on observations of sunspots 
see Casanovas (1997) and references therein. 

3.4 Structure of sunspots 

The responsible for the origin and structure of sunspots is believed to be the deep- 
seated toroidal (azimuthally oriented) magnetic flux in the solar interior (Schussler 
et al, 2002). This magnetic flux is thought to be situated between the convective 
zone and the radiative interior and increases by means of a dynamo action up to 
the formation of a kind of instability. The eruption of flux leads to the formation 
of fi-shaped magnetic loops coming up to higher layers, therefore reaching the 
convective zone and emerging through the solar surface as active regions (Zwaan, 
1992). Summarizing this idea, what we have is that the magnetic field generated 
underneath the convective zone in the Sun, emerges to the solar surface forming 
sunspot groups and active regions and spreads around its surface through convec- 
tive mechanisms. 
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Diverse models for sunspots have been proposed through the years such as the 
cluster model and the monolithic model. On the one hand the cluster (spaghetti) 
model developed by Parker (1979) in which the magnetic field divides into many 
separate flux tubes in the first few Mm below the surface. On the other hand the 
monolithic model by Cowling (1957) establishes that the magnetic field underneath 
the solar surface is confined to a single flux tube and the plasma flows crossing it 
are prohibited. 

Sunspots, as complex magnetic structures embedded in a convective plasma, 
show many active and changing features on multiple scales. Convective flows and 
large-scale plasma circulation plays an important role in the dynamics and evolution 
of solar active regions (see, e.g., Schrijver & Zwann, 2000) 

Important physical processes are present in sunspots, some of them at spatial- 
scales beyond the resolution limit of the current solar telescopes. In the last years 
and using the latest restoration techniques and advances in adaptive optics, the 
complex magnetic and dynamic structure of sunspots have been unveiled in some 
detail. Nevertheless, there are many important questions concerning them that 
remain unanswered yet. 

3.4.1 Umbra 

The umbra is the central dark area that corresponds to the coolest part of a sunspot. 
The temperature inside the umbra is low (~ 3500-5000 K) due to the presence of 
a strong magnetic field that inhibits the convective energy transport underneath 
(Biermann, 1941). The magnetic field orientation inside the umbra is vertical 
respect to the solar surface and becomes slightly inclined when approaching the 
umbra-penumbral boundary. Larger sunspots with bigger umbrae exhibit a more 
intense magnetic field strength and the range inside umbrae goes from 2000 to 
3500 Gauss in average. The energy radiation in the umbra is roughly 20% of the 
photospheric value. 

Although the umbra seems to be completely dark and with non internal struc- 
ture, when observed in high-resolution, a lot of fine features in its interior can be 
detected, for instance the umbral dots which are bright dots with sizes of about 
0.2-0.5 arc seconds and hotter than the umbral background. These features are 
still under study to clarify and confirm whether they are related to convective phe- 
nomena transporting some sort of flow. Light bridges are other structures inside 
the umbra that look like luminous streaks or narrow bands with different ages and 
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sizes. 

3.4.2 Penumbra 

The penumbra is the filamentary and deeply-rooted (thick) structure in a sunspot 
that forms the ring-shaped part surrounding the umbra, where a significant part 
of the total magnetic flux is carried (Schmidt, 1991; Thomas & Weiss, 1992). 
The penumbra is seen in continuum images as radial alternating bright and dark 
filaments (e.g. Collados, del Toro Iniesta & Vazquez, 1988) and exhibits an in- 
termediate brightness between the umbra and the quiet photosphere (see also the 
reviews by del Toro Iniesta, 2001; Solanki, 2003). Figure 3.2 shows some high 
resolution images taken with the SST in La Palma, where we can see the intricate 
filamentary structure of the penumbra and the fine features inside and around pores. 

Magnetic field is in general stronger in the inner part closest to the umbra 
and decreases in the outer frontier and its energy radiation is approximately 75% 
the photospheric value. The picture of the fine structure of the penumbra has 
not always been the same but evolved during the last years. Nowadays it is com- 
monly accepted that the penumbral magnetic field is uncombed so that at least 
two different inclinations of the magnetic field coexist on a small scale in a sort of 
interlocking structure (see Figure 3.3). 

The main two distinct inclinations correspond to the component where the 
magnetic field is tilted by ~ 40-50 degrees (with respect to the vertical to the 
solar surface) with a strength of ~2000 G, and the component in which the mag- 
netic field is almost horizontal and weaker (Lites et al, 1993). Vertical magnetic 
inclinations were first associated to bright filaments whereas the more horizontal 
ones to dark filaments though nowadays this picture has changed significantly since 
the discovery of two narrow lateral brightenings located at each side of the central 
obscuration, as will be commented below. 

Recent works using high-resolution images have improved the measurements of 
the field inclination in the penumbral filaments (Bello et al, 2005; Langhans et al, 
2005). Results from full-Stokes profiles have confirmed and firmly established that 
variations of inclination and field strength are indeed present in sunspot penum- 
brae, and moreover they are anti-correlated: the less inclined components (known 
as spines) have stronger magnetic field than the more horizontal ones (know as 
intraspines, Lites et al, 1993). In the intraspines, some flux tubes (collection of 
multiple magnetic field lines) overtake the external penumbral border (Sainz Dalda 
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& Martinez Pillet, 2005) whereas some other immerse below the surface due to 
the turbulent effect of the granular convection pushing them downwards (Westen- 
dorp Plaza et al . , 1997; Thomas et a 1 , 2002). We will come back again to this topic. 

Although many unknowns have been unveiled, there are still many fundamental 
and unanswered questions such as why sunspots have a penumbra and how it is 
maintained for a few days before it disappears. A lot of internal structure has been 
resolved in penumbral filaments and they seem to be composed by a dark core 
accompanied by a bright halo at both sides. This new features named dark-cored 
penumbral filaments, recently discovered by Scharmer et al. (2002), have brought 
great interest as they can unveil the very fundamental influence of magnetocon- 
vection at small scales occurring in sunspot penumbrae. At the beginning of the 
filament, the dark core is not present and the lateral brightenings merge into a 
larger bright structure known as penumbral grains (PG) that seems to move radi- 
ally in the penumbra. 

Several authors have worked on this new observational evidence to investigate 
on the fine structure of the penumbra and nowadays there some models of the 
penumbral fine structure trying to explain all diverse phenomena. Firsly, the most 
mentioned and confronted one, the uncombed model proposed by Solanki & Mon- 
tavon (2003) in which the penumbra is seen as a collection of horizontal flux tubes 
embedded in an almost vertical and strong background magnetic field. Numer- 
ical simulations of the temporal evolution of flux tubes (the moving-tube model 
of Schlichenmaier, Jahn & Schmidt, 1998a, b) have modelled them starting at the 
magnetopause, then rising adiabatically as a result of radiative heating, by magnetic 
buoyancy, and eventually being bent into the horizontal where they cool down due 
to radiative losses while continue outward from the center of the spot. Recently, 
some new observational and theoretical works support the uncombed model 1 (e.g. 
Muller et al, 2002; Bellot Rubio et al., 2003; Borrero et al, 2005; Bellot Rubio 
et al, 2007; Jurcak & Bellot Rubio, 2008; Ruiz Cobo & Bellot Rubio, 2008). A 
second model has been proposed by Spruit & Scharmer (2006), the so-called gappy 
penumbral model (see also Scharmer & Spruit, 2006). These authors interpret the 
penumbral filaments as due to convection in field-free, radially aligned gaps below 
the surface of the penumbra and argue that this solves the discrepancy between the 
large heat flux and the low velocities present in penumbra or, in other words, that 
their model can explain the heating of the penumbra. Apart from these two models, 
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alternative ones have come out on the scene based on Micro-Structured Magnetic 
Atmosphere (MISMA, Sanchez Almeida, 1998, 2005). These models assume that 
the penumbra is formed by optically thin magnetic fibrils a few km in diameter and 
moreover that every resolution element contains a bunch of unorganized field lines 
with random strengths and inclinations. In general, all models have some problems 
(see Bellot Rubio, 2007) and the debate is not sealed yet. 

New features and fine-scale structures have been recently observed thanks to 
the improvements in the real-time and post-facto correction of the images. It is 
likely that inside dark cores there is more unseen fine structure. This new structures 
will add more ingredients to this already complicated scenario. 

The penumbra not only shows an intricate magnetic structure but also a com- 
plex dynamics. The PGs for instance, are moving in two privileged directions, 
radially towards the umbra in the inner penumbra with velocities of the order of 
500 m s _1 (Muller, 1976) and towards the surrounding photosphere in the outer 
penumbra (Sobotka et al., 1999; Sobotka & Sutterlin, 2001). The line dividing 
both distinct directions of movement is placed at a distance from the outer penum- 
bral border of about one third of the total penumbra longitude. Some of the PGs 
cross the penumbra-photosphere boundary and continue to move in the granulation 
as small bright features (Bonet et al., 2004). 

Apart from the radial proper motions of PGs, Marquez, Sanchez-Almeida &i 
Bonet (2006) find, by using local correlation tracking techniques, that radial mo- 
tions diverge away from bright filaments to converge toward dark filaments, leading 
to an association between downflows and dark features that could be a sign of 
convection in penumbra. This result is reinforced by Sanchez Almeida et al. (2007) 
on the basis of spectroscopic measurements, with high-angular resolution, of the 
Evershed effect (see next section). 



3.4.3 Evershed flow 

The Evershed Flow (hereafter EF) is associated to an observational effect in the 
penumbra registered as a global wavelength shift for spectral lines forming in the 
penumbra of sunspots. There are actually two kinds of EF registered: 1) the nor- 
mal EF at photospheric level where blueshifts are observed in the penumbra on the 
sunspot side close to the solar disc center whereas redshifts are detected on the 
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FIGURE 3.2 — Image showing a highly- resolved sunspot penumbra in a sunspot: AR 10030, 
observed the 15 July 2002 with the 1-meter SST in la Palma (left) and some fine-scale structures 
(right) placed inside and around solar pores. Courtesy Institute of Solar Physics, Stockholm. 

sunspot side towards the solar limb. It corresponds to mainly radial and horizontal 
flows traveling in the penumbra of the sunspots towards the external photosphere; 
2) at chromospheric level, the so-called inverse EF, where the material flows to- 
wards the umbra, contrary to the normal EF and much less studied so far. From 
here forth, we will focus on the normal EF. 

This effect was first observed by Evershed in 1909 (Evershed, 1909) at Ko- 
daikanal Observatory in India and is interpreted as a Doppler shifting caused by the 
radial and predominantly horizontal flow along the penumbra. The EF has been 
interpreted as a siphon flow along magnetic flux tubes (Meyer & Schmidt, 1968; 
Thomas, 1988). 

There have been several studies trying to explain the properties and nature of 
the EF in the last forty years (e.g. Haugen, 1969; Maltby, 1975; Moore, 1981; Di- 
aletis, Mein & Alissandrakis, 1985; Dere, Schmieder & Alissandrakis, 1990). For a 
very recent review the reader is referred to the work by Bellot Rubio (2007). 

In more recent studies and helped by high-resolution observations, it has been 
confirmed that there is a strong correlation between the EF and the horizontal 
magnetic fields in the penumbra. It is the expected behavior as resulted from mag- 
netohydrodynamical considerations that require for the flow to be aligned respect 
to the magnetic field lines (Bellot Rubio et a I . , 2003). Rimmele (1995) found that 
the EF is confined in thin channels which elevate upon the continuum level over 



92 



CHAPTER 3. Solar active regions 




FIGURE 3.3 — Cartoon showing the interlocking structure of the magnetic field in the filamentary 
penumbra of a sunspot. Radial bright filaments where the magnetic field is inclined (~ 40 ° respect 
to the horizontal direction in the external penumbra) coexist with the dark filaments having an 
almost horizontal magnetic field. Courtesy John H. Thomas. 

the vast majority of the penumbra. The last result has been afterwards confirmed 
in different observations. It has also been found that some of these channels and 
their associated field lines re-enter the solar surface in the outermost part of the 
penumbra and beyond (Westendorp Plaza et al., 1997; Schlichenmaier & Schmidt, 
2000). Some recent results indicate that the EF does not end abruptly at the 
white-light sunspot boundary, but that at least part of it continues beyond the 
visible boundary and forms an extensive magnetic canopy (e.g. Solanki, Montavon 
& Livingston, 1994; Rezaei et al, 2006). The existence of radial outflow in the 
sunspot canopy is still an ongoing debate. 



3.5 Photosphere surrounding sunspots 

Sunspots are immersed in the convective plasma as mentioned above, and this 
plasma is observed mainly as granulation. Granulation is the first pattern that can 
be identified when looking at the solar photosphere and covers the vast majority of 
its surface. Figure 3.4 shows the underneath convective pattern mechanism forming 
the granulation that we see on the photospheric layer. In the inter-granular lanes 
we can detect the presence of magnetic flux coming from below, sometimes visible 
as bright and elongated structures as shown in the high-spatial resolution image in 
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Figure 3.4. 

The granulation convective pattern surrounding sunspots is perturbed by the 
presence of magnetic elements that move radially outwards through an annular cell 
called the "moat". Next sections comment on the physical properties of these 
magnetic structures and on the plasma flows around sunspots. 



3.5.1 Moving Magnetic Features 

The Moving Magnetic Features (MMFs) are bright features that correspond to small 
magnetic elements of mixed polarity traveling radially outwards while immersed 
in the granulation surrounding sunspots (Sheeley, 1972; Harvey & Harvey, 1973; 
Hagenaar & Shine, 2005, for a recent review). These features move with velocities 
of some hundredths of meters per second up to 3 km s _1 and can be generally 
classified in three main groups, as follows: 

• Type I: Bipolar pairs of magnetic elements moving together outwards the 
sunspot inside the moat, with velocities in a range between 0.5 a 1 km s" 1 . 
They appear just outside the penumbra along directions that are extensions 
of the dark penumbral filaments. These bipolar pairs are interpreted as being 
footpoints of a loop emerging to the solar surface belonging to a penumbral 




FIGURE 3.4 — Sketch showing the 
underneath convective pattern mech- 
anism in a solar granulation region. 
The black arrows plotted over the 
granulation image stand for the diver- 
gence directions of the exploding gran- 
ule whereas the short black arrows be- 
low show the position of the global up 
and downflows. Magnetic lines (in or- 
ange) are also delineated and corre- 
spond to small flux tubes seen as bright 
structures in the inter-granular areas 
when reaching the photospheric level. 
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FIGURE 3.5 — Diverse interpretations for the bipolar Type I MMFs 

flux tube in which there is a strong vertical convection. The loop is then swept 
away by the moat flow. Several loops can emerge in different positions and 
different times along the same flux tube and this mechanism is seen when 
successive MMF Type I tend to follow similar trajectories across the moat 
region. As long as the flux tube is attached to the sunspot, this type of MMFs 
associated to the tube are not responsible for the sunspot decay. Figure 3.5 
shows the diverse interpretations of bipolar type I MMFs (Yurchyshyn et al, 
2001; Wilson, 1973; Harvey & Harvey, 1973). 

• Type II: Consist of unipolar magnetic elements with the same polarity as the 
parent sunspot that are shed from the sunspot border, moving across the 
moat with velocities of 2 and 3 km s _1 . They are interpreted to be flux 
tubes that have been separated from the main flux and are dragged by the 
moat. The loss of flux represented by this features is a possible mechanism 
involved in the sunspot decay (Hagenaar & Shine, 2005). Nevertheless, 
further investigations should be done to know whether or not this type of 
MMFs is more common in rapidly decaying sunspots. 

• Type III: The last type of MMFs are also unipolar features but in this case 



3.5. Photosphere surrounding sunspots 



95 



with magnetic polarity opposite to that of the parent sunspot, having high 
velocities in a range similar to that of the Type II. They are more transitory 
than the first two types and travel outward a few thousand km before disap- 
pearing, and are thought to be produced when a flux tube section emerges 
to the surface generating a rapid horizontal displacement of the footpoints 
in the moat. 

3.5.2 Moat flow 

The existence of an annular region around sunspots was first reported by Sheeley 
(1969) who termed it moat. Radially oriented outflows in moats were found by 
Sheeley & Bhatnagar (1971). 

The sunspot moat defines an organized horizontal flow pattern that ends quite 
abruptly at a distance that can be comparable to supergranular sizes or even larger. 
Sheeley (1972) reported Doppler speed of these outflows 0.5 - 1.0 km s _1 for dis- 
tances 10-20 Mm beyond the outer boundary of the sunspot. Brickhouse & Labonte 
(1988) performed a study of the properties of moats and concluded that the moat 
radius is roughly proportional to the penumbral radius. 

The characteristics and temporal evolution of moats derived at two heights in 
the atmosphere for a large sample of sunspots with different sizes, shapes and evo- 
lutionary stages have been recently studied by Sobotka & Roudier (2007). These 
authors establish a relation between sunspot types and their occurrence, areas, and 
horizontal velocities of the moats in the photosphere and in the transition region 
and found that moats do not show substantial changes for periods of about 12 
hours. Table 3.1 lists the averaged horizontal velocities inside and outside the 
moat for young and old spots and Figures 3.6 and 3.7 show plots (for young and 
old spots) of the maximum moat radius versus the sunspot radius, and the average 
azimuthal distribution of the the moat area and velocities respectively. 

Bonet et al. (2005) have found that MMFs in the moats are mostly located 
in radially oriented channels between the local divergent motions of mesogranular 
scale. This motions confine the MMFs within the channels. They also found that 
most MMFs move passively with the same average speeds as the neighbouring 
granules, that is, they are carried along by the same large-scale moat flows as the 
granules. However, a small fraction of them (6%) move away from the spot more 
rapidly (they reach velocities grater than 1.4 km s^ 1 ) than the neighbouring gran- 
ules which can be interpreted as the consequence of interactions with fast moving 
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TABLE 3.1 — Averaged horizontal velocities [m s 1 ] inside and outside moat. Taken from 
Sobotka & Roudier (2007). 



Position 


Passband 


Young spots 


Old spots 


Inside moat .... 


white light 
1600 A 


380 ± 30 
450 ± 70 


410 ± 30 
450 ± 40 


Outside moat ... 


white light 
' 1600 A 


200 ± 40 
200 ± 40 


240 ± 40 
210 ± 50 



small granular fronts or fragments. 

New findings (Sainz Dalda & Martmez Pillet, 2005) have shown that the penum- 
bral filaments extend beyond the sunspot boundary entering the region dominated 
by the moat flow where the MMF activity is detected. Thus, the temporal av- 
erage of magnetograms has unveiled the existence of moat filaments: horizontal, 
filamentary structures coming from the penumbra and reaching the photospheric 
network as an extension of the penumbral filaments. Moreover, some MMF have 
been found (Sainz Dalda & Martmez Pillet, 2005; Ravindra, 2006; Cabrera Solana 
et a I . , 2006) starting just inside the outer penumbral boundary in its way out from 
the sunspot (see Figures 3.8 and 3.9). Kubo et al. (2007) have recently established 
a relationship between the vertical component (spines) of the magnetic field in the 
so-called uncombed structure of the penumbra (section §3.4.2) and MMFs (with 
the same polarity as the sunspot) observed in the moat regions. 
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FIGURE 3.6 — Plot of the maximum moat radius versus sunspot radius. Young spots are marked 
by "+" and old spots by "{>"■ Linear fits are shown for young spots (dashed line) and old ones 
(dash-dot). Solid line correspond to values R ma x = r s . Taken from Sobotka & Roudier (2007). 
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FIGURE 3.7 — Plot of the average azimuthal distribution of the moat area (left) and average 
azimuthal distribution of velocities inside moats and outside (right) for old spots (solid line) and 
young sunspots (dotted line). Taken from Sobotka & Roudier (2007). 
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FIGURE 3.8 — Temporal evolution of the center-side penumbra and moat region as studied by 
Cabrera Solana et al. (2006). The figure shows (from top to bottom): Continuum intensity at 
1565 nm, magnetogram signal from Fe I 630.25 nm, logarithm of the total linear polarization of 
Fe I 1564.85 nm, Stokes V zero-crossing velocity of Fe 1565.28 nm, and Stokes V zero-crossing 
velocity of Fe I 630.25 nm. Negative velocities are blue shifts. Pink and blue contours delimit 
features A and B respectively. Red lines indicate the inner and outer penumbral boundaries. Tick 
marks in x-axis represent 1". The black arrow points to the disc center and the time scale starts 
at June 30, 9:34 UT. 



FIGURE 3.9 — Average of the absolute value (left) and true sign magnetogram (right) showing 
the presence of moat filaments all around the sunspot. Black arrows point out some examples. 
The mean umbral and penumbral boundaries are overlaid. Taken from Sainz Dalda & Martinez 
Pillet (2005). 




Proper motions in a complex active 

region 



his chapter is the first one devoted to the analysis of solar observations, showing 



I the results from the application of the image restoration process commented 
in the former part of the thesis. The study is based on ground-based observations 
of a complex ^-configuration solar active region in which we performed a detailed 
analysis of proper motions at photospheric level describing the dynamics of the 
plasma in the FOV. Some results presented in this chapter are an extensive version 
of Vargas Dommguez et al. (2007). 

4.1 Introduction 

Sunspots that contain umbrae of opposite magnetic polarity within a single penum- 
brae are known as ^-configuration sunspots. They are complex active regions that 
usually evolve rapidly and are likely to produce energetic solar flares (Zirin & Wang, 
1989). The present chapter is based on an excellent ~79 minutes time series of a 
sunspot group with a ^-configuration. The long duration of the series, besides the 
high and stable quality throughout the entire period, substantially improved after 
image reconstruction, make of this material an excellent data set to study morphol- 
ogy and dynamical behavior of sunspots and their surroundings. We are interested 
here in the relation between the existence of penumbra in this active region and 
the dynamics of the convective plasma surrounding it. The proper motions inside 
penumbrae are also analyzed and a statistical study between the different regions 
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over the FOV is performed. 
4.2 Observations 

The first observing campaign I was involved in at the Swedish 1 m Solar Telescope 
(SST 1 , Scharmer et al., 2003a) corresponding to the so-called Spanish Time 2 , was 
held in April 2004. During this campaign I was particularly interested on becom- 
ing familiar with the instrumentation, software and observing techniques that have 
been in constant evolution and development in recent years. Moreover, the ob- 
served images were used to test the restoration algorithms (see chapter §1). 

The second observing campaign was held at the SST in July 2005 and corre- 
sponded to the so-called International Time Program (ITP) 3 . In this framework, 
several Japanese scientists came to the Canary Observatories 4 for a joint pro- 
gramme of simultaneous observations using various solar telescopes: the Vacuum 
Tower Telescope (VTT), the Solar Swedish Telescope (SST), the Dutch Open Tele- 
scope (DOT), and the satellites Transition Region and Coronal Explorer (TRACE) 
and Solar and Heliospheric Observatory (SOHO). The observations were focused on 
targets displaying solar active regions and, in particular, regions inside and around 
sunspots where various magnetic features (e.g. MMFs, see section §3.4) are widely 
observed. The SST and the DOT were devoted to imaging, pursuing observations 
of the photosphere and chromosphere with high -spatial and -temporal resolutions. 
The main objective of these coordinated observations was to figure out what kind 
of magnetic features in the photosphere/chromosphere are responsible for quasi- 
steady heating and enhanced activities well observed in the active corona. 

4.2.1 Instrumentation and setup 

The work and results presented in all this chapter are based on the observations 
done at the SST (during the ITP collaboration, see Figure 4.1), using the adaptive 

1 The SST is operated by the Institute for Solar Physics (ISP) of the Royal Swedish Academy 
of Sciences at the Roque de los Muchachos Observatory in La Palma, Canary Islands. 

2 The Spanish Time corresponds to the 20% of the total observing time at the Observatories 
of the Canary Islands, assigned to Spain by the institutions the telescopes belongs to, by virtue of 
the International Agreements on Cooperation in Astrophysics. 

^Period corresponding to the 5% of the total observing time (see previous footnote), aimed to 
promote astrophysical research in cooperation with various telescopes and institutions. 

4 Observatorio del Teide (OT) in Tenerife and Observatorio Roque de los Muchachos (ORM) 
in La Palma. 
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FIGURE 4.1 — Left figure: Sketch showing the design of the Swedish 1-m Solar Telescope (SST) 
installed at the Roque de los Muchachos Observatory (more information available at the website 
http://www.solarphysics.kva.se/). The diagram plots the primary focus (A) where the rays coming 
from the lens converge, the Schupmann system (6) to correct for chromatic aberration, and the 
deformable mirror of the adaptive optics (C). Courtesy Institute of Solar Physics, Stockholm. 
Right figure: Photography of the SST taken by the author during the 2005 campaign. 
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SST 


CHARACTERISTICS 


Type 


Evacuated refractor 


Mount 


Turret design 


Aperture 


97 cm 


Primary focus length (fa) 


20.35 m 


Secondary focus length (fa) 


45.8 m 


Wavelength (A) range 


0.35 fim - 0.8 fim 


Image scale in fa 


4.5"/mm 


Post- focus instruments 


Correlation Tracker (CT) 




Adaptive Optics (AO) 




High resolution CCDs 




SOUP Tunable Filter 




Spectrograph 


Capabilities 


Low-order AO 




(up to 35-modes 




Karhunen-Loeve) 




Real-time frame selection 



TABLE 4.1 — Main parameters and characteristics of the SST. 



optics (AO) system. The AO system has a bi-morph mirror equipped with 37 elec- 
trodes and a Shack-Hartmann wavefront sensor with 37 microlenses. A thin-film 
filter (dichroic filter) separates the light in two beams, the blue and the red beams. 
The blue beam for imaging and the red one to feed a longitudinal magnetograph 
(the SOUP Tunable Filter that will be described in detail later on). Table 4.1 enu- 
merates the main characteristics of the telescope. 

In order to apply the restoration techniques described in chapter §1 we recorded 
two objects: images in G-band 5 and G-cont 6 . The setup included then, for the 
blue beam, two channels taking simultaneous images. In the first channel, G-cont 
(A436.3 nm, FWHM=1.15 nm), two Kodak Megaplus 1.6 CCD cameras with a 
10-bits dynamical range and 1536 x 1024 pixels were employed, as corresponding 
to the PD focus-defocus image-pair. The defocus image was generated by moving 
backwards the camera a certain distance to generate a larger optical path in com- 
parison to the focus one, hence inducing a known amount of defocus. The second 



5 Spectral band of molecular lines at A430.5 nm. 

6 Relatively clear continuum, close to the G-band, at A436.3 nm. 
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FILTER 


G-band 


G-cont 


Ca II H 


A [A] 


4305.6 


4363 


3969 


AAfwhm [A] 


10.8 


11.5 


3 


Exposure time [ms] 


~13 


-13 


~80 


CCD 


1536 x 1024 


1536 x 1024 


2048 x 2048 


Field of view ["] 


63 x 42 


63 x 42 


84 x 84 


Pixel size ["] 


0.041 


0.041 


0.041 



TABLE 4.2 — Parameters for the blue channel during the 2005 campaign. 



channel was devoted to G-band (A430.5 nm, FWHM=1.08 nm). The pixel size of 
the cameras was 0'.'041 square. 

Apart from the two above described, a third independent channel (asynchronous 
respect to G-band and G-cont) was acquiring reference images in the Ca II H line 
(A396.9 nm, FWHM=0.3 nm). Table 4.2 groups the main characteristics for the 
blue channel. 

Figure 4.2 shows the setup with all the optical elements used in the configuration 
we mounted for the campaign. The optical paths and instruments employed by the 
adaptive optics system are drawn in red in the figure. On the right-hand side, the 
blue beam includes the PD image-pair (G-cont), the G-band channel (synchronized 
with G-cont by means of a chopper) and the independent Ca II H channel. 

The Solar Optical Universal Polarimeter (SOUP, Title et al, 1986) tunable filter, 
owned by the Lockheed Martin Solar and Astrophysics Laboratory (LMSAL), was 
placed along the red beam. This instrument is essentially a temperature compen- 
sated Lyot filter. It operates by rotating a series of linear polarizers and 1/2-wave 
retarders placed between calcite crystal optics. The selected wavelength is adjusted 
by tuning the servo motors for each rotating element to a previously calibrated 
position. Two bandpass modes, wide (>100 mA) and narrow (< 100 mA), are 
selectable for any given line tuning. The complete configuration of the SOUP also 
includes two Liquid Crystal Variable Retarders (LCVRs) for measuring the polar- 
ization state of the incoming light. The SOUP filter is used for observations in the 
Ha spectral line that forms in the solar chromosphere and to obtain magnetograms 
using the magnetically sensitive line Fe I 6302 A. 

A beam splitter in front of the SOUP divides the light beam to obtain simul- 
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SST - La Pal ma Setup June/July 2005 
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FlGURE 4.2 — Setup mounted at the SST for the July 2005 observing campaign. The optical 
paths and instruments employed by the adaptive optics system are drawn in red. On the right-hand 
side the blue beam includes the PD image-pair (G-cont), the G-band channel (synchronized with 
G-cont by means of a chopper) and the independent Ca H channel. 
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FIGURE 4.3 — The Solar Optical Universal Polarimeter (SOUP) tunable filter (left) and the 
science high-resolution cameras (right) on the SST optical bench. 



taneous images in the continuum near the selected bandwidth. These images, 
understood as additional objects, will help us to restore the narrow band images 
(by MOMFBD) and to obtain high-resolution magnetograms. 

Figure 4.3 shows the SOUP (left) and the science high-resolution cameras of 
the blue beam (right) as placed on the optical bench at the SST 7 . 

4.2.2 Data acquisition 

During the period 3-13 July 2005, images were acquired as listed in Table 4.3. 
The weather conditions at the observatory through all the campaign were quite 
variable obtaining rather medium quality images though for the observing run on 
July 9 we were able to record very good images (close to the diffraction limit after 
restoration) due to an excellent seeing. 

Our targets of interest were sunspots in different evolutionary stages close to 
the solar disc center both, leaders and followers. Within all the images recorded, 
a post-facto pre-selection was made in order to get rid of the data below a certain 
quality threshold. The quality of the images is quantified by averaging the local 
gradients in absolute value inside a granulation window within the FOV. The quality 
was measured by using the G-band images as reference. Figure 4.4 shows a plot 



7 The complete description of the instrumentation and the manuals of operation of the SST 
can be found in the website http://www.solarphysics.kva.se/LaPalma/ 
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TT1VTF, TTT 


APTTVF RFdTOlST 


OTT ATJTY 


03.07.05 


8:38-11:47 


10781 


medium 


04.07.05 


16:45-17:31 


10781 (follower) 


poor 


05.07.05 


7:58-9:15 


10781 (follower) 


medium 


05.07.05 


14:16-14:44 


10781 (follower) 


good 


05.07.05 


16:23-17:33 


10781 (leader) 


medium 


06.07.05 


7:50-8:57 


10786 (/eacfer) 


good 


07.07.05 


8:55-10:12 


10786 (leader) 


medium 


07.07.05 


10:36-12:58 


10786 (east of the leader) 


good 


07.07.05 


13:40-14:25 


10786 (leader) 


medium 


08.07.05 


8:00-8:45 


10786 


medium 


09.07.05 


7:44-9:30 


10786 


excellent 


10.07.05 


8:18-9:07 


10786 


medium 


10.07.05 


16:18-16:47 


10786 (north region) 


good 


11.07.05 


7:44-10:01 


10789 (follower) 


medium 


12.07.05 


8:30-9:53 


10789 


medium 


13.07.05 


8:36-9:42 


10789 (follower) 


good 



TABLE 4.3 — Log of observations for the 2005 International Time Program (ITP) campaign. 
The quality of the images is quantified by averaging the absolute value of the local gradients inside 
a window located in a granulation region within the FOV. The quality values have been translated 
into a more qualitative scale varying from poor quality images up to excellent quality ones. The 
quality was measured using the G-band images as reference. 
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3000 — 
2000 ^ 



8.0 8.5 9.0 

Time [hours, decimal fraction] 

FIGURE 4.4 — Plot of the quality (in arbitrary units) of the G-band images during the entire 
observation on 9 July 2005. The quality scale is such that 4000 indicates good quality images 
and there are astonished good quality images having values over 7000. At the end of the recorded 
images (at ~ 9:30 UT) the quality went down to values around 3000 which forced us to stop the 
observations due to the bad seeing conditions after that. 



of the quality for the G-band images through the entire observation on July 9, 2005. 

Among all the acquired data, we decided to work with the best quality images 
that corresponded to the NOAA Active Region 10786 observed on July 9, 2005. 
This complex region, displaying a (^-configuration, was placed at heliocentric posi- 
tion /i =0.9 as shown in the SOHO full-disk continuum image in Figure 4.5. The 
sequences span for more than one hour, from 7:47 UT to 9:06 UT, following the 
evolution of the sunspots. 

4.3 Data processing 

Real-time corrections with the AO of the telescope and further post-processing 
techniques rendered up images near the diffraction limit. The main steps performed 
in the processing of the images are sequentially detailed in the next sections. 

4.3.1 Flatfield and darkfield corrections 

Standard procedures for flatfielding and dark-current subtraction were firstly ap- 
plied. During the observations, apart from the science images, we also acquired the 
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FIGURE 4.5 — Solar active region 
NOAA 10786 observed on July 9 2005 
(small black box). Image taken by 
the SOHO satellite. The white box 
indicates the SOHO high-resolution 
FOV and the disc center is marked 
with a black dot. 



corresponding flat and dark images 8 . The corrections were made by applying the 
following expression 



where i corr is the final image after flatfield and dark current corrections, i raw is 
the observed raw image, d and / are the mean dark and flat images, respectively, 
and indexes (i, j) stand for the pixel position in the image matrix. Special care has 
been taken during these corrections to eliminate dark and hot pixels, and also the 
black borders in the flat images, as a crucial step to ensure the images will not be 
affected by spurious artifacts when applying the restoration code. 

4.3.2 Restorations 

Post-processing for image restoration was performed by employing the Multi-Object 
Multi- Frame Blind Deconvolution (MOMFBD) described in section §1.5. In the 
blue beam, the observing strategy to apply this technique consisted in taking G- 
band images and simultaneous G-cont phase diversity image-pairs. Sets of about 
18 images per channel (i.e. 3 x 18 images) within a temporal cadence of 10.0517 

8 To improve the SNR, every flat and dark image resulted from real-time averages of 100 single 
exposures. The flats were taken while the telescope, out-of-focus, scanned the central area of the 
solar disc at fast speed. 



corr 



I 



f(i,j) - d(i,j) 



(4.1) 
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G-BAND G-CONT 




FIGURE 4.6 — Example of the final images in G-band (left panels) and G-cont (right panels) 
after applying the restoration process. Displayed are the worst (upper panels) and the best (middle 
panels) quality images out of the 18 images in a particular set among the 472 sets. Using the 
MOMFBD method, two restored images are generated, one for the G-band (lower left panel) and 
another one for G-cont (lower right panel). The total FOV in the original images is 63x42 
whereas for the final ones reduces to 60x40 ". 
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CHANNEL 


IMAGES Time, UT 


G-band 

G-cont (focus) 
G-cont (defocus) 


221-8818 7:47:45 - 9:09:15 
115-8712 7:47:45 - 9:09:15 
228-8825 7:47:45 - 9:09:15 



TABLE 4.4 — Images information in the restoration process. Tabulated in the second column 
are the number of the initial and final images for each time series in every channel. 

s, were combined to produce a pair of simultaneous G-band and G-cont restored 
images. Table 4.4 shows the information about the images employed in the restora- 
tion process. 

Following this procedure, two time series of reconstructed images were pro- 
duced: 472 images for each channel, G-band and G-cont, that correspond to 79 
min of continuous solar observation. Figure 4.6 displays an example of the restora- 
tions pursued, for one out of the 472 image groups processed. The total FOV in 
the original non-restored images was 63x42 " whereas the final FOV in the final 
restored images was 60x40 ". 

4.3.3 De-stretching and subsonic filtering 

Once the images were restored, they were, firstly, de-rotated to compensate for 
diurnal field rotation observed in the focal plane, induced by the alt-azimuthal con- 
figuration of the telescope and, secondly, aligned by cross correlation. 

Additional image processing included the so-called de-stretching. The stretch- 
ing or distortion effect on the images, commented in section §1.1, consist on dif- 
ferential local displacements over the FOV. They are due to differential deviations 
from the original tilts in the wavefronts coming from different point sources. At- 
mospheric turbulence is the responsible for this problem. The effect is specially 
strong during the morning time when the line-of-sight (hereafter LOS) to the Sun 
is more inclined hence inducing a greater optical thickness in the light path and 
consequently, a larger atmospheric turbulence contribution. 

Thus, the stretching is a local effect manifested with different strength at dif- 
ferent positions over the FOV and can be corrected at medium- and large-spatial 
scales by applying local correlation techniques. The local correlation tracking (LCT) 
technique was originally designed for the removal of distortions in image sequences 
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(November, 1986) and later on, was used to infer the proper motions of the gran- 
ulation in a time series of white-light images (November & Simon, 1988). 

Due to the strong distortion present in our time series, the de-stretching al- 
gorithm was applied repeatedly in four subsequent steps, to completely eliminate 
residual effects. In every step, the size of the correlation tracking window was mod- 
ified from greater to smaller sizes to ensure corrections at different spatial scales. 

Subsonic filtering is the last stage in the corrections applied to the images and 
deals with the cleaning process pursued to eliminate high frequency oscillations, 
such as the p-modes and residual jitter stemming from de-stretching, in the fre- 
quency k — w domain (Title et al, 1986). A filter with a phase-velocity threshold 
of 4 km s" 1 was employed with satisfactory results (see Figure 4.7 center). The 
sharp cutoff in the subsonic filter was replaced by a smooth transition from 1 to 0. 

An apodization in the temporal dimension (8% of the total length at each 
end) was also applied to minimize border effects in the filtering process. The final 
product were two movies (G-band and G-cont) 9 of 428 frames each, spanning over 
71 min with a cadence of 10.0517 s, and covering a FOV of 57'/8 x 34'.'4. The time 
series reached the high-resolution level shown in Figure 4.7 {lower central panel). 
This quality is pretty much suitable for subsequent scientific analysis. Figure 4.8 
shows the superb quality achieved in the best image of the time series. 



4.3.4 Maps of displacements 

The LCT method has been commonly used for determining velocity fields in time- 
resolved image data and was first applied by November &i Simon (1988) to measure 
proper-motions in solar granulation. Since then, different authors have used LCT 
techniques at many different spatial scales to study the solar dynamics, such as 
long-term evolution and flows in supergranulation (Shine, Simon & Hurlburt, 2000; 
De Rosa, Duvall & Toomre, 2000a). 

The method is based on a simple but powerful idea that works on the principle 
of finding local concordances between two frames (images). A correlation window 

9 See the movies at the web site: http://www.iac.es/galeria/svargas/movies.html 
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FIGURE 4.7 — Steps in the correction of a time series of images. The correction of stretching 
eliminates the distortion in the images (upper panel) caused by atmospheric turbulence. The 
subsonic filtering eliminates the intensity fluctuations produced by solar high-frequency oscillations 
and limits the phase-velocity of the solar structures to a maximum of 4 km s _1 (middle panel). 
After the described corrections, the time series is ready for its physical analysis (lower panel). 




FIGURE 4.8 — G-band best quality image, taken at the SST on 9 July 2005 during the ITP 
campaign, after the data processing including the MOMFBD restoration. 
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selects corresponding sub-areas in both images (the patches to be correlated). The 
selected sub-areas include those structures we want to track, such as solar gran- 
ules that evolve in time. The displacement vector is given by the difference of the 
coordinates of the sub-area centers when the best match between them is found. 
Mathematically this idea is expressed in terms of a correlation function between 
two instants. 

The selection of the sub-areas is performed by means of a window which is usu- 
ally set to be a Gaussian function with a Full-Width at Half Maximum (FWHM) 
adapted to the desired size for the correlation patch and depending on the spatial- 
scale of the structures we are about to track. 

It is well documented in the literature that LCT-techniques in general produce 
some systematic errors in the determination of displacements. The nature of the in- 
terpolation algorithm to fix the position of the correlation maximum with sub-pixel 
accuracy may play a role in the results. The LCT-method typically underestimate 
the velocity fields. Tests performed by the authors of the LCT-code we are employ- 
ing in this work lead to conclude that in the worst cases this underestimation can 
amount to 20-30 % (Yi, 1992; Molowny-Horas, 1994). 

Figure 4.9, taken from November & Simon (1988), displays a plot of the pre- 
cision of the displacement measurements by means of LCT and using a quadratic 
interpolation formulae to determine the correlation maximum. These authors have 
used solar granulation images recorded on photographic film support, digitized at 
12 positions with increments of 0.1 pixel. Then, they have computed displacement 
maps by comparison with the unshifted image applying their LCT algorithm. Fig- 
ure 4.9 shows the ratio between the computed and the real displacements. The 
former being the average of displacements over the area of each displacement map 
with error bars representing ±lcr. 

4.3.5 Magnetograms 

In the red beam, the SOUP tunable birefringent filter working at Fe I 6302A in 
combination with a modulator composed by liquid crystal retarders, produced lon- 
gitudinal magnetograms as mentioned above. The SOUP filter recorded images in 
the two flanks of the line at ±50 mA from its core. The modulator switched be- 
tween two states to obtain 10 left and 10 right circular polarization states (LCP and 
RCP) at each line position. The total time required to record all four-type images 
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D . 5 1,0 1.5 



FIGURE 4.9 — Precision of the measurements of displacement by LCT as shown by November 
& Simon (1988). A quadratic interpolation formulae has been used to determine the correlation 
maximum. See the text for more information. 



was 33 s. A beamsplitter in front of the SOUP filter deflected 10% of the light 
to obtain broad-band phase diversity image-pairs, that combined with the narrow- 
band images of the SOUP allowed MOMFBD reconstructions. The MOMFBD 
restoration was therefore performed from subsets of 40 SOUP images including all 
polarization states and line positions, and 40 simultaneously recorded wide-band 
PD image-pairs, that is, a total of 120 images to produce 5 restored ones: I R l (f P , 
IlcP' IrcP' Ilcp anc ' Icont- Magnetograms (M) and also Dopplergrams (D) can 
be generated from the saved restored images by means of the following expressions 



2 



D = i 
2 



jblue 
1 RCP 



jblue 
1 LCP 



jred 
1 RCP 



jred 
1 LCP 



jblue I jblue 

1 RCP ~r 1 LCP 



jred _i_ jred 
1 RCP ~r 1 LCP 



' jblue jred jblue jred 

1 RCP 1 RCP , 1 LCP 1 LCP 

jblue _i_ jred jblue i jred 

1 RCP 1 RCP 1 LCP 1 LCP 



(4.2) 
(4.3) 



The Doppler velocity v is given by: v = f c D, where f c is a constant to convert 
Doppler signal to Doppler velocities and is estimated to be ~ 10 km s _1 by com- 
paring the Doppler signal with the mean LOS Doppler velocity of granules (~ 1 
km s _1 ). The FOV of the restored SOUP images is 90" x 60" (with a pixel size of 
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FIGURE 4.10 — Magnetogram corresponding to the region under study obtained with the SOUP 
filter in the SST on 9 July 2005. The <5-spot neutral lines in Figure 4.11 were obtained from this 
frame. The coordinates are expressed in Mm. 

0'.'061 pixel -1 ). Figure 4.10 displays one of the resulting magnetograms having a 
spatial resolution of about 072. 

4.4 Results 

4.4.1 General description of proper motions 

The G4oand series was employed to study the proper motions of the structures in the 
FOV by using a LCT algorithm implemented by Molowny-Horas & Yi (1994). For 
our analysis, we have preferably chosen a Gaussian tracking window with FWHM 
0'.'78 being half of the typical granular size, which is suitable for tracking mainly 
granules. 

With this code we obtained two maps of horizontal displacements or proper 
motions per time step (i.e. the horizontal velocity components v x and v y ), which 
were averaged in time. Averages of horizontal velocity components were performed 
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over 5-min and 71-min intervals. 

Figure 4.11 shows the resulting flow map averaged over the whole series (71 
min). Moat flows are seen in the spot on the lower left of the FOV but more 
prominently in the spot on the right side in the upper and lower penumbral regions 
(moats flows will be treated in more detail in section 4.4.2). Several centers of 
diverging horizontal motions are present in the entire FOV. Some of them surround 
sunspots (lower left and right sides of the figure) but the most conspicuous ones are 
seen in the lesser magnetized area (lower part of the FOV between the spots), dis- 
playing greater velocities and more symmetrical shapes. These velocity structures, 
related to recurrent expansion and splitting of granules, are commonly associated 
with mesogranulation (Roudier et al., 2003; Roudier & Muller, 2004; Bonet et al., 
2005). 

The white lines in Figure 4.11 delineate the LOS neutral lines inferred from a 
SOUP magnetogram (see Figure 4.10). The regions crossed by these lines corre- 
spond naturally to horizontal fields, similar to those found in penumbrae. However, 
we caution that these horizontal fields can have different flows as those that occur 
in normal penumbra (such as shear flows, as found by Deng et al. (2006) or the 
supersonic nozzles observed by Martinez Pillet et al. (1994). Close inspection also 
shows that areas near neutral lines have a filamentary appearance that is much less 
distinct than observed in normal penumbrae extending radially from umbral cores. 

We measured a mean speed in the granulation outside the moats of ~0.35 km 
s _1 . This small value reflects the presence of high magnetic signal (dense pop- 
ulation of pores and G-band bright points all over the field), however, it is even 
smaller than the reported, also for magnetized granulation, by other authors that 
employ tracking windows of similar widths (Berger et al., 1998; Bonet et al., 2005). 

The region in between the left and right hand side sunspots in the active region 
under study is indeed a plage region of abnormal granules full of magnetic activity as 
revealed by Ishikawa et al. (2007). Figure 4.12 shows the results from these authors 
working with the central region of the FOV (in an upside-down position respect to 
the one shown in Figure 4.11). The G-band image shows a large amount of bright 
structures that are selected and marked in white (see Figure 4.12 center). These 
G-band bright points are ubiquitously distributed over the whole plage region and 
point out the intense magnetic activity in this so-called magnetized granulation. 
The corresponding magnetogram (see Figure 4.12 right) also reveals the magnetic 
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FIGURE 4.11 — Map of horizontal velocities in the entire FOV (71 min average). The black 
contours outline the borders of umbrae, penumbrae and pores. The white lines delineate the LOS 
neutral lines inferred from a magnetogram. The coordinates are expressed in Mm. The length 
of the black bar at coordinates (0,0) corresponds to 4 km s" 1 . The background represents the 
average image of the G-band series. 



4.4. Results 



121 




" X(«n»eei) " X(anssee9) 



FIGURE 4.12 — Left: G-band image displaying the central part of the 5-configuration active 
region (turned upside-down respect to the one in Figure 4.11 and scales in arc seconds). Center: 
G-band bright points are selected and marked in white. Right: magnetogram (M) of the same 
region in gray scale ranging from -2795 (black) to 1000 (white) Gauss. 

polarities behavior and intense activity in the same region. 

Inside penumbrae the proper motions were studied to infer the systematic di- 
rections of the flows. Figure 4.13 shows, with enlarged velocity scale, the proper 
motions inside penumbra using a tracking window FWHM= $.'78 and averaging 
over 71 minutes. These displacements do not necessarily correspond to real mo- 
tions of plasma but could reflect displacements of brightness. We are aware that 
the tracking window used allows us to determine only averaged displacements of 
the penumbral fine features, however, the description of the proper motions of 
this features is out of the scope of the present work. Figure 4.13 reveals mean 
flows at spatial resolution of ~ l"(mean velocity ~ 0.4 km s -1 ), moving toward 
the umbra in the inner penumbra and toward the surrounding granulation in the 
outer penumbra. The filamentary like structures in grey, mapping small horizon- 
tal velocities (< 0.13 km s" 1 ), outline the frontiers separating both kinds of flow. 
This clear division of the velocities field in two quite differentiable tendencies has 
already been mentioned by some authors (see Bonet et al., 2004) but not in a sys- 
tematically way inside sunspots in a complex active region as presented in this work. 

In order to perform a statistical study of horizontal velocities, a sequence of 14 
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FIGURE 4.13 — Proper motions inside penumbrae (FWHM= 0'.'78, 71 min average). The figure 
displays the close-to-zero horizontal speed (vh < 0.13 km s _1 ) as gray traces inside white-colored 
penumbrae. These lines delimit the changes in the privileged flow directions: 1) flows towards the 
umbra in the inner penumbra, and 2) towards the surrounding granulation in the outer penumbra. 
Umbrae and pores are delimited by black contours. The coordinates are expressed in Mm. The 
length of the black bar at coordinates (0,0) corresponds to 4 km s _1 . 
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REGION 



mean 



rms 



Penumbrae 



0.40 0.35 



Magnetic Gran. 0.35 0.24 
Moats 0.67 0.32 



TABLE 4.5 — Statistics of the horizontal speeds Vh [km s _1 ] for various regions within the FOV 
of a complex solar active region. The calculations are done using the velocity magnitudes derived 
from 14 maps (5 min averaged, FWHM 0'.'78) covering 70 minutes of solar evolution. 



velocity maps was obtained averaging over 5-min windows along the whole series. 
As expected, the resulting maps are a little more noisy but with a wider range of 
variation of the speeds (i.e. velocity magnitudes) as compared with the case of 
averaging over 71 min. Figure 4.14 shows the histogram of the speeds as measured 
in this sequence of maps, for different regions within the FOV: penumbral areas, 
granulation outside the moats and moat regions (see also Table 4.4.1 for mean 
and rms values). Note that the moat histogram in the Figure 4.14 shows a cutoff 
imposed by the lower threshold selected to define moat velocities. 

The statistical analysis of the results (Figure 4.14 and Table 4.5) renders a 
value for the mean speed in the moats of ~ 0.67 km s _1 , similar to that reported 
by Bonet et al. (2005), that is almost twice the mean value in the rest of the 
granulation field. The confinement of such significantly higher speeds in restricted 
areas close to the sunspots makes moats a sharply differentiated structure in the 
photospheric granulation. 

We pursued the same statistical calculations for velocity maps computed with 
a tracking window with FWHM tf!2. The results revealed that when reducing the 
size of the tracking window, the measured mean and rms speed values increment 
since we are also including small structures, as bright points, moving faster than 
the general flows observed in the granulation pattern evolution. In our experiment, 
for instance, the difference between the mean and rms values is about twofold when 
going from FWHM 0"78 to 0'.'2. 



124 



CHAPTER 4. Proper motions in a complex active region 



Histograms of horizontal speeds 




a 1 2 3 4 



V h [km.s' , l 



FIGURE 4.14 — Histogram of horizontal speeds Wj, (from 14 maps, 5 minutes averaged each) 
for a 5-configuration active region. A local correlation tracking technique has been employed to 
derive the horizontal velocity components using a Gaussian tracking window of FMWH= 0'.'78. 
Table 4.5 lists the statistics of velocity magnitudes in every region. 
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4.4.2 Large-scale horizontal granular flows 

To study the proper motions surrounding the sunspots and in particular the be- 
havior of the large-scale flows in the FOV, a velocity threshold can be selected to 
only plot a given range. The map in Figure 4.15 shows only those velocities in the 
granulation field with magnitude above 0.4 km s _1 . Strong radial outflows (moats) 
are evident surrounding the sunspots and 0.4 km s _1 appears to be a character- 
istic speed value to define the moats. The black contours in the figure outline 
the area of the moats. It was drawn by hand following visual criteria of proximity 
to the sunspot and avoiding a few strong exploding events at the lower part of 
the FOV that contaminate the large-scale moat flow. Empty areas in the moat 
region in the velocity maps are a result of the chosen value for the speed threshold 
and correspond to the interaction with exploding events that occur inside the moat. 

From inspection of Figure 4.15, it becomes evident that the moats are closely 
associated with the presence of a penumbra. The examples in the lower left, lower 
right and upper right corners are quite revealing. The right pointed velocity vectors 
on the lower left corner correspond to the moat flow of a penumbra that is not 
visible in the FOV of this figure, but can be seen in an image covering a wider FOV, 
of the same active region in Figure 4.16. Moats are missing in the sunspot sides 
with no penumbral structure. This is readily seen in the right side of the umbral 
core at coordinates (9,10), but also in the left side of the one at (25,18). None of 
the pores in the middle of the FOV is associated with any moat-like flow. More 
importantly, moats are also absent in granulation regions that are located side by 
side with a penumbra, but in the direction perpendicular to that marked by the 
penumbral filaments. A clear example is the region around coordinates (9,6-7) in 
Figure 4.15. As a normal (meaning here far from the neutral line) penumbra, it is 
a candidate to develop a moat flow, but none is seen in this region. The moat flow 
only appears along the direction delineated by the penumbral filaments. Moats, 
thus, appear to be a natural extension of the flow along the direction of penumbral 
filaments and do not exist as outflows in the transverse direction. 

4.4.3 Neutral lines affecting the flow behavior 

As commented in section §4.3.5 we produced magnetograms of the FOV by means 
of the SOUP filter. The high-resolution magnetogram in Figure 4.10 shows the two 
polarities. The positive polarities (in white) dominate the left and right areas of 
the FOV and correspond widely to the polarities of the respective sunspots. The 
central area of the image is dominated by a negative polarity (in black). Two neu- 
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FIGURE 4.15 — Map of the horizontal velocities with magnitude > 0.4 km s _1 (71 min average). 
The contrast within penumbrae has been enhanced by removal of a spatial running mean of the 
original image. Strong radial outflows (moats) are evident surrounding filamentary penumbrae. 
They are confined by the black contours. These moats are found to be closely associated with 
the existence of penumbra. Note that moat region surrounding coordinates (1,5) corresponds to 
a penumbra located to the left of the image and not visible in our FOV. The coordinates are 
expressed in Mm. The length of the black bar at coordinates (0,0) corresponds to 4 km s _1 . 
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tral lines can be then identified separating opposite polarities. These neutral lines 
are also overplotted in the magnetogram, dopplergram and white-light image taken 
by the DOT during the same observing run but with a wider FOV (see Figure 4.16). 

It is seen that when the neutral lines run close to a penumbral border and par- 
allel to it, the expected behavior of the outflows in the granulation region side by 
side with that penumbra changes somehow and the moats do not appear in those 
regions. The penumbral region in Figure 4.15, centered at coordinates (5,20), does 
not show clear evidence of a moat flow (although localized outflows are seen above 
and below it). This region corresponds to a strongly sheared neutral line (see Fig- 
ure 4.11) commonly seen in 5 spots (Deng et al., 2006). In the SOUP magnetogram 
(Figure 4.10) the spot centered on the left side of the FOV (white polarity) does 
not show a radially outward directed penumbra, but penumbral filaments that run 
parallel to the neutral line. This is a configuration commonly referred as sheared. 
This configuration is also inferred from the filaments observed in the G-band frames. 
We suspect that both conditions, the presence of a neutral line and the absence 
of radially oriented filaments, are responsible for the absence of a clearly developed 
moat flow. 

The nature of neutral lines with sheared configurations in 5-spots is not well 
understood. But we stress that if they correspond to relatively shallow structures, 
they would also block the upward propagation of heat and should be able to generate 
a moat flow. Such a moat flow is not observed in this sunspot. 

4.4.4 Motions of centers of divergence 

Following November & Simon (1988), (see also Marquez, Sanchez-Almeida & 
Bonet, 2006), the divergences and the vertical velocities can be calculated from 
the horizontal velocities, according to the expression 

v z (x,y) = h m V • Vh(x,y), (4.4) 

where h m =150 km stands for the scale height of the flux of mass and V • Vh{x, y) 
the divergence of the horizontal velocity field. 

For this purpose, we generated four horizontal velocity maps (Gaussian track- 
ing window of FMWH= (f!78) averaging over time windows of 5, 15, 40 and 70 
min, and starting from the beginning of the time series. From these maps we ob- 
tained the divergences and the corresponding vertical velocities. Since we were not 
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FIGURE 4.16 — Magnetogram (upper panel), Dopplergram (middle panel) and white light image 
(lower panel) for the 5-configu ration sunspot group, observed by P. Siitterlin with the DOT on 9 
July 2007. Neutral lines are identified in the magnetogram (white lines) and overplotted in the 
Dopplergram (white lines) and in the white light image (black lines). 
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AVERAGED: 5 MIN [km s" 1 ] AVERAGED: 15 MIN [km s" 1 ] 

min max rms mean min max rms mean 

-7.77 12.41 1.78 0.04 -6.30 10.64 1.46 0.03 




AVERAGED: 40 MIN [km s" 1 ] AVERAGED: 70 MIN [km s" 1 ] 

min max rms mean min max rms mean 

-4.99 7.80 1.35 0.01 -4.23 6.15 1.22 0.01 




FIGURE 4.17 — Maps of vertical velocities averaged over different time windows starting from 
the beginning of the time series (5, 15, 40 and 70 minutes, FMWH= 0'.'78). In order to identify 
general trends in the vertical flows, the maps have been smoothed by convolving with a Gaussian 
running window of FWHM=l'.'47. Yellow represents close-to-zero velocities. Extreme values of 
divergence can be easily identified: exploding granules with positive divergences (in dark red) at 
coordinates (30,0) and powerful sinks with negative divergence (in blue) at the upper right part 
of the image. Regardless the average applied, from 5 up to 70 minutes, the exploding events and 
sinks are recurrently seen and clearly defined within the FOV. The statistics of vertical velocities 
in the FOV is presented for every averaged map. Black contours outline umbrae, penumbrae and 
pores. Note that the color-scales (the representation in pseudo-color is the opposite to the one 
used for Doppler shifts) are slightly different for every map since the values for the minimum and 
maximum velocity change from one case to the other as shown by the statistics tabulated at the 
top of each map. The coordinates are expressed in Mm. 
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FIGURE 4.18 — Histogram of vertical velocities in the FOV, calculated from the 5-min averaged 
velocities maps. 



interested in describing details but only general trends in the granulation vertical 
flows, the resulting v z (x,y) maps were smoothed by convolving with a Gaussian 
running window of FWHM = l'. / 47, so that the identification of coherent structures 
was facilitated. Figure 4.17 displays the resulting four vertical velocity maps as im- 
ages representing the velocity magnitudes by using a pseudocolor scale. Structures 
corresponding to positive divergences (i.e. recurrent exploding granules) and sinks 
with negative divergences are coherently distributed all over the FOV for all the 
different temporal averages though an evident evolution can be noticed between 
these averages. The figure shows extreme cases for both positive divergences at 
coordinates (30,0) in dark red, and sinks in the upper right part of the image in blue. 

For the particular case of 5-min horizontal velocity averages we computed 14 
of these maps covering the whole time series (70 min), thus deriving 14 vertical 
velocity maps similar to the one in the upper left corner in Figure 4.17. Figure 4.18 
shows the histogram of vertical velocities as calculated from these 14 maps, that 
shows a quasi-symmetric behavior around the mean value (v ~ 0). 

The divergence structures displace in time dragged by flows at scales larger than 
the mesogranular ones. In order to follow the evolution of the centers of divergence 
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along the time series, we employed the same 14 vertical velocity maps described 
in the previous paragraph 10 . Then, by applying the LCT again but this time upon 
the sequence of 14 images of divergence, we can determine how the centers of 
divergence displace in time. Because of the larger size of the structures we want 
to track, the correlation tracking window was set to FWHM 3". 

Figure 4.19 displays the resulting map corresponding to the averaged velocities 
of the centers of divergence during the whole time series (70 min). The veloci- 
ties inside sunspots and pores are set to zero since we are only interested in the 
granulation proper motions. Figure 4.19 reveals the smooth character of the flow 
transporting the centers of divergence; the flows are more organized and uniform as 
compared with those in Figure 4.11 but they still preserve signatures of mesogran- 
ular flows. The well-organized moat flows are also clearly defined in the figure as 
they were in Figure 4.11. Note that the scales employed to represent the velocity 
vectors in Figures 4.19 and 4.11 are quite different so that the velocities displayed 
in the former are much lower than those represented in Figure 4.11. 



4.4.5 Study of convective cells 

Using the map of averaged velocities calculated in section §4.3.4 from the 428 
images conforming the time series (Figure 4.11), we can study the evolution (dis- 
placement) of passive corks homogeneously distributed in the FOV at time equal 
zero, one per pixel. By applying the mentioned averaged velocity field to the spread 
corks we can follow the displacement of every single cork in time steps of 10.0517 s. 
With this procedure, a sequence of images showing the locations of corks in every 
moment can be generated. Figure 4.20 displays the position of corks (in white) 
within the FOV after 0, 22, 45, 67, 89, 111, 134 and 156 minutes respectively. It 
is important to differentiate between the sequence of 428 solar images conforming 
the time series and the sequence of images showing the evolution of the corks. 

This evolution can be extended to times longer than the total duration of the 
series (71 min). This is done by assuming the same map of velocities affecting the 
corks, as shown in the last 4 panels of Figure 4.20. 



10 For the purpose we are interested here, the distinction between vertical velocity maps and 
divergence maps is irrelevant since both quantities only differ by a constant factor h m . 
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FIGURE 4.19 — Map of velocities of the centers of divergence, calculated from a sequence of 
fourteen 5-min averaged maps of divergence structures, along the whole time series (71 min). 
The velocity magnitudes are in a range from to 1.18 km s _1 . The length of the black bar at 
coordinates (0,0) corresponds to 1 km s _1 . 
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Some properties concerning the displacement of corks in the FOV can be inferred 
from Figure 4.20. Regions having stronger horizontal velocities are quicker emptied 
and correspond to the large-scale flows (moats) presented in Figure 4.15. After 22 
minutes there are evidences of mesogranulation and exploding granules affecting the 
displacements of the corks. Already in minute 45 the mesogranular flows dominate 
the FOV, these mesogranules appear as elongated structures when located inside the 
moat flows due to the dragging effect. After 67 minutes when the mean lifetime of 
mesogranules has been reached or even overtaken, the mesogranular structures tend 
to disappear, except in the upper central part of the FOV where the intense magnetic 
field confine and limits the horizontal proper motions of the plasma. A considerable 
number of corks are distributed along dark penumbral filaments as pointed out by 
Marquez, Sanchez-Almeida & Bonet (2006). In the more quiet granulation regions 
(lower central area of the FOV) and in the last panels of Figure 4.20 the remaining 
bright structures trace a pattern resembling supergranular structures which mean 
lifetime is superior to the duration of the present sequence. 



4.5 Conclusions and discussion 

Time series of G-band and G-cont images of a sunspot group with ^-configuration 
(NOAA active region 10786), spanning over 71 min, were observed and corrected 
for atmospheric and instrumental degradation. Proper motions in the entire FOV 
have been measured in the G-band series by means of local correlation tracking 
techniques. In this chapter, we concentrate in the results obtained in the identified 
moat flows (or lack thereof) of several structures within the FOV. The main con- 
clusions from the analysis can be summarized as follows: 

1) We have identified inside the penumbra frontiers separating mean longitudi- 
nal flows (parallel to the direction of the filaments) at spatial resolution ~1" (mean 
velocity about 0.4 km s^ 1 ), moving toward the umbra in the inner penumbra and 
toward the surrounding granulation in the outer penumbra. 

2) We have detected strong (mean speeds of 0.67 km s _1 with an rms of 0.32 
km s _1 ) outflows streaming from penumbrae radially oriented from an umbral core, 
the so-called sunspot moats. Our velocity values are comparable to those reported 
by i.e, Sheeley (1972), Vrabec (1974), Harvey & Harvey (1973), Hagenaar & Shine 
(2005) and Bonet et al. (2005), which describe the outflow velocities in the range 
0.5 - 1 km s _1 . Umbral cores sides with no penumbra do not display moat flows. 
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FIGURE 4.20 — Evolution of artificial passive corks moving under the influence of the ve- 
locities averaged over the whole time series (71 minutes), displayed in Figure 4.11. In 
one of the lower corners of each panel is written down the elapsed time (in minutes) from 
the beginning. A movie with the animation of the corks is available in the web site: 
http:/ / www. iac. es/galeria/svargas/corks.html 
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3) Furthermore, the moats are also absent on penumbral sides parallel to the 
direction defined by the penumbral filaments. They are not found in directions 
transverse to them. A special case is sheared penumbral configurations tangential 
to the umbral core. No moat flow is found there either. 

This evidence is clearly suggestive of a link between the moat flow and flows 
aligned with the penumbral filaments. The possible connection with the Evershed 
flow is inescapable. Although it can be argued that a physical relationship between 
these two flows is not firmly established in this chapter, we believe the evidence 
is clear-cut and that a statistical analysis of a larger number of regions, as done 
here, can establish this connection more solidly. This is one of the goals in the next 
chapter, where we extend the sample of solar active regions in order to consolidate 
the main conclusions of the present chapter. 

4.5.1 Implications of the results 

The results in this chapter should be put in the context of the recent findings 
by Sainz Dalda & Martinez Pillet (2005); see also Ravindra (2006) mentioned in 
chapter §3. These authors find that the penumbral filaments extend beyond the 
photometric sunspot boundary and cross the region dominated by the moat flow. 
This region is also the location where the MMF activity is normally found. Indeed, 
some of these MMFs are seen by these authors to start inside the penumbra. All 
these results point to an origin linking the moat flow and MMF activity with the 
well-known Evershed flow. Cabrera Solana et al. (2006) also suggest that the Ev- 
ershed clouds inside the penumbrae propagate to the moat, becoming MMFs once 
they leave the sunspot. It is indeed somewhat paradoxical that while the final fate 
of the Evershed flow has remained unknown for decades, independent explanations 
and physical scenarios have been proposed for the generation of the moat flow (that 
starts exactly where the Evershed flow is seen to vanish). 

Similarly, the results obtained from local helioseismology near sunspots, in par- 
ticular those related to the f-mode (Gizon et al., 2000), are also quite relevant to 
the results presented here. The interpretation of these results in terms of convective 
cells surrounding sunspots (e.g., Bovelet & Wiehr, 2003) may be different if the 
Evershed flow turns out to be the major process that injects mass into the moats 
surrounding sunspots. 



5 

Moats flows surrounding sunspots 

This chapter is an extended and more detailed version of the paper by Vargas 
Dommguez et al. (2008) published in the The Astrophysical Journal (ApJ). The 
analysis is performed employing ground-based observations of solar active regions, 
and complements the ideas already presented in the previous chapter concerning 
the plasma motions at photospheric level. 

5.1 Introduction 

The sunspot moat, as previously commented in this thesis, consists of an organized 
horizontal radial flow pattern surrounding the sunspots, which ends quite abruptly 
at a distance that can be comparable to supergranular sizes or even larger. 

In the last chapter the moat flow in a complex active region was studied and 
described in detail. As a result, some conclusions were established relating these 
large outflows and certain physical characteristics of the sunspots: a) No outflow 
was detected in the granulation next to umbral cores that lack penumbrae, b) 
Outflows were only found in the granulation regions adjacent to the penumbrae 
in the direction following the penumbral filaments, c) Granulation regions located 
next to penumbral sides parallel to the direction of the filaments show no moat flow. 

The aim of the present chapter is to extend the study done in chapter §4 to a 
larger sample of cases, that is, to establish whether this moat-penumbra relation 
is systematically found in other active regions and how the granulation convective 
pattern surrounding sunspots behaves. 
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In doing so, we have used several time series of images observed at high-spatial 
resolution. A total of 7 different sunspot series have been processed and analyzed. 
The sample includes sunspots with different penumbral configurations, varying from 
well-developed penumbrae to rudimentary penumbral morphologies. 

5.2 Observations and data processing 

The observations were obtained with the Swedish 1-m Solar Telescope (SST, Scharmer 
et al., 2003a) on La Palma between 2003 and 2006. All observations benefited from 
the use of the SST adaptive optics system (Scharmer et al., 2003b) that minimized 
the degrading effects of seeing. Image post-processing techniques were applied to 
increase the homogeneity in the quality of the time series and to enhance the image 
quality over the whole FOV. For the 2003 data set we applied Multi-Frame Blind 
Deconvolution (MFBD) using the implementation developed by Lofdahl (2002). 
The MFBD code was succeeded by the Multi-Object Multi-Frame Blind Decon- 
volution code (MOMFBD, Van Noort, Rouppe van der Voort & Lofdahl (2005)) 
which employs multiple objects and phase-diversity (see section 1.5 for more details 
on this method). MOMFBD was applied to the data sets after 2003. For all the 
time series we present in this chapter, the seeing conditions were generally very 
good and sometimes excellent, being able to record for more than 40 min in all 
cases. A large fraction of the restored images in the different time series approached 
the diffraction limit of the telescope. 

After image restoration, we applied standard techniques to the time series in- 
cluding correction for the diurnal field rotation, rigid alignment and de-stretching 
to correct for seeing-induced image warping. A p-modes filtering was also applied 
to the series (phase-velocity threshold of 4 km s" 1 ). 

Table 5.1 presents details on the different active region targets which are also 
shown in Figure 5.1. The observations were performed in regions at different posi- 
tions on the solar disc as listed in the columns 4 and 5 in Table 5.1, including one of 
them (S2) very close to the limb. Table 5.2 gives some details on the restored time 
series. The duration of the time series is in all cases longer than 40 minutes so that 
we can use them to get a reliable average of the displacements of the structures. 
Below we provide more detailed information on the different data sets, including 
the restoration process applied in each case. 
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FIGURE 5.1 — Solar active regions, from sunspot SI to S7, showing the different penumbral 
configurations analyzed in this chapter. The spatial scale is not the same in all the sunspots. 
The last two panels show a representation of the solar hemispheres including the positions of 
the observed sunspots. The dashed arcs determine the locations for ji =0.88 ( the outer one in 
the northern hemisphere), y, =0.99 (the inner in the northern hemisphere), and /i =0.55 (in the 
southern hemisphere) respectively. 
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Table 5.1 — Sum 


ipots Sample 




Name 


Active region 
NOAA 


Observing date 


Heliographic coordinates 


/" 


SI 


10440 


2003, Aug 18 


(S10, Wl) 


0.96 


S2 


10608 


2004, May 10 


(S5, E57) 


0.55 


S3 


10662 


2004, Aug 20 


(N13, E7) 


0.99 


S4 


10662 


2004, Aug 21 


(Nil, W5) 


0.99 


S5 


10789 


2005, Jul 13 


(N17, W32) 


0.88 


S6 


10813 


2005, Oct 04 


(S7, E37) 


0.78 


S7 


10893 


2006, Jun 10 


(Nl, E17) 


0.95 



5.2.1 SI: AR440, 22-Aug-2003 

An interference filter centered at the continuum at 436.4 nm (FWHM 1.1 nm) was 
used. In the remainder, this filter will be referred to as the "G-cont" filter. Every 
25 s, the three highest contrast images were selected and stored to disk. These 
three images were used for MFBD processing. This sunspot is shown in Figure 5.1 
(SI). 

5.2.2 S2: AR608, 10-May-2004 

An interference filter centered at the G-band at 430.5 nm (FWHM 1.2 nm) was 
used. In the remainder, this filter will be referred to as "G-band" filter. This 
observation involved two cameras set-up as a phase-diversity pair: one in focus and 
one slightly out of focus. 20 exposures from each camera, acquired over the course 
of 7.5 seconds, were used for MOMFBD processing. During processing of the time 
series, a few bad quality images, suffering from too much blurring, were dropped. 
The exposure time was 11 ms. This sunspot is shown in Figure 5.1 (S2). 

5.2.3 S3: AR662, 20-Aug-2004 

This time series was observed with the G-cont filter. The MOMFBD processing 
was performed using 10 exposures (acquired during 7 s) from a phase-diversity 
pair of cameras. For the central part of the FOV, centered on the sunspot, the 
MOMFBD processing was extended to involve a G-band phase-diversity pair of 
cameras. After MOMFBD restoration, the best quality image from 3 subsequent 
images was selected for further processing of the time series. The exposure time 
was 11.3 ms. This sunspot is shown in Figure 5.1 (S3). 
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5.2.4 S4: AR662, 21-Aug-2004 

This time series was processed in a manner similar to S3. After MOMFBD restora- 
tions, the worst 14% of the images were dropped from the time series. This sunspot 
is shown in Figure 5.1 (S4). 

5.2.5 S5: AR789, 13-Jul-2005 

For this time series, MFBD processing was applied to 18 subsequent G-band images, 
acquired during 10 s. Images that were too blurred and when the AO system was 
not actively compensating, were dropped from the MFBD sets. This sunspot is 
shown in Figure 5.1 (S5). 

5.2.6 S6: AR813, 4-Oct-2005 

This time series was observed with an interference filter centered on the spectral 
region between Ca K and Ca H at 395.4 nm (FWHM 1.0 nm). For the MOMFBD 
processing, 20 exposures (acquired during 8.3 s) from a phase-diversity pair of 
cameras were used. In addition, simultaneous exposures from two other cameras 
were used for the restoration: one equipped with a narrow-band filter centered in 
the Ca H wing at 396.5 nm and one equipped with a narrow-band filter centered on 
the Ca H core (396.8 nm). The exposure time was 11 ms. This sunspot is shown 
in Figure 5.1 (S6). 

5.2.7 S7: AR893, 10-Jun-2006 

This time series was observed with an interference filter centered at 630.2 nm 
(FWHM 0.8 nm). This is the pre-filter for the Lockheed SOUP filter. MOMFBD 
processing was applied to 400 exposures obtained during 19 s. The restorations 
included images from a phase-diversity pair of cameras with the wide-band filter, 
and narrow-band SOUP exposures in the Fe I 630.2 nm spectral line. This sunspot 
is shown in Figure 5.1 (S7). 

5.3 Data analysis and results 

The high quality, stability and long duration of the restored sunspot series, enable 
us to study the dynamics of the plasma around the sunspots. We have computed 
proper-motion velocity fields (horizontal velocities) employing the LCT technique 
described in the previous chapter. A Gaussian tracking window of FWHM l'.'O 
suitable for mainly tracking granules has been used for all the series. The method 
produces a sequence of map-pairs each one describing, within a short time interval 
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Table 5.2— 


Restored sunspots time series 


Name 


N. images 


Cadence [s] 


Duration [min:s] 


SI .... 


128 


24.7 


52:47 


S2 .... 


376 


7.5 


47:10 


S3 .... 


144 


19.8 


47:20 


S4 .... 


344 


8.0 


45:57 


S5 .... 


240 


10.1 


40:12 


S6 .... 


556 


8.7 


80:30 


S7 .... 


124 


19.7 


40:26 



(tens of milliseconds), the horizontal velocity components, v x and v y , along the 
FOV. These maps of components have been averaged over 5 and 10 min periods, 
and also over the respective total duration of every sunspot series (more than 40 
min in all cases). From the averaged maps of velocity components, the distribution 
of velocity magnitudes over the FOV (flow map) has been easily derived. 

As described in section §4.3.4 , it is well documented in the literature (November 
& Simon, 1988) that LCT-techniques in general produce some systematic errors 
in the determination of displacements. Their significance depends also on several 
factors like the width of the tracking window in relation to the size of the elements 
used as motion tracers, and the presence of large intensity gradients in the images 
(see section §4.3.4 for more details). Even so, this drawback does not change the 
main conclusions of the present work since we are not interested in fixing absolute 
values of velocities but rather in the detection of large-scale regularly organized 
flows around the sunspots in comparison with those typically found in a normal 
granulation field. In both cases we use the solar granules as the tracers of motion 
and the same size for the tracking window so that in case of some bias it will affect 
to both velocity fields similarly. 

5.3.1 De-projection of horizontal velocities 

The sample of sunspots studied in this paper includes a variety of heliocentric po- 
sitions on the solar disc as shown in Table 5.1 (with fi = cos#, where 6 is the 
heliocentric angle). Away from the solar disc center, the measured proper motions 
are in fact projections of the real horizontal velocities in the sunspot plane onto 
the plane perpendicular to the LOS. Thus, to evaluate the real horizontal velocities 
we have de-projected the velocities obtained from the LCT technique to a plane 
tangent to the solar surface at the center of the sunspot. 
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To that aim, let us consider in Figure 5.2 the orthogonal coordinate system 
SX,SY,SZ (the Sunspot System, SS) with the SZ-axis perpendicular to the plane 
tangent to the solar surface at the sunspot location and the SX-axis tangent to 
the meridian of the solar sphere at this location. Also included in the drawing is 
the coordinate system X,Y,Z (the Observing System, OS) with a common origin, 
the Z-axis coinciding with the LOS and the X-axis pointing to the axis crossing the 
Sun center in a direction parallel to the LOS. Note that the axis, X, SX, Z, and SZ 
are coplanar and SY and Y are collinear. The real horizontal velocities, v(v, 4>), are 
contained in the (SX, SY) plane whereas we observe their projections v'(v',4>'), in 
the (X,Y) plane. 

Calculations based on Figure 5.2 lead to the following relationships (equations 
5.1 and 5.2) between the magnitudes v and v', and the azimuths <f> and cj)', that 
will allow us to de-project our observations. 

v' 2 = (v sin (f>) 2 + (v cos <f> cos 9) 2 = v 2 (sin 2 <f> + cos 2 4> cos 2 6). (5-1) 



, vsincf) tancj) 

tan ^ = 2 a = o- v 5 - 2 ) 

v cos <p cos cos 

Since the X- and Y- axes do not necessarily coincide with our natural reference 
system (i.e. the edges of the CCD frame), the measurement of the projected az- 
imuths, <j)' ', requires the knowledge of the orientation of our FOV with respect to 
the X-axis that points to the solar disc center. 

5.3.2 Masking moats 

As mentioned above, for every sunspot series, the map of horizontal velocities aver- 
aged over the whole time series is evaluated. In order to coherently detect the moats 
and compare their statistics (section §5.3.6) in all active regions of our sample, we 
have de-projected the observed velocity vectors onto the sunspot plane (SX,SY) as 
described in section §5.3.1. 

Overlaying the observed images, we construct maps showing in the granulation 
field only those projected velocity vectors having de-projected magnitudes above 
a certain threshold. In these maps, extensive organized outflows coming out from 
the sunspots can easily be identified so that we can outline masks delimiting in a 
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LOS 




FIGURE 5.2 — Sketch showing two orthogonal coordinate systems for the projection analysis 
of the horizontal velocities: 1) The Sunspot System (SX.SY.SZ) with the Z-axis perpendicular 
to the solar surface at the sunspot location, and 2) the Observing System (X,Y,Z), where the 
observations are described, with the Z-axis coinciding with the LOS. 
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qualitative way the area of the moats for the purpose of statistical calculations and 
graphical representation. 

We have fixed the mentioned threshold to ~ 0.3 km s _1 as the best compro- 
mise to define the limits of the moats. This value should not be understood as an 
absolute (universal) velocity threshold to define moats in general. As we have men- 
tioned before, the particular method used here to measure the displacements may 
underestimate the velocity magnitudes. Moreover, outside the moat one can easily 
find velocities larger than that. The point is that this particular threshold makes 
in our case more evident the existence of a well organized radial outflows around 
the sunspots in comparison with the rest of the FOV. Lower values (even zero) for 
the threshold do not extend the areas of organized flows around the sunspots but 
produce maps with a very dense and noisy (exploding granules everywhere) rep- 
resentation of arrows where the outline of the frontiers of the moats results more 
difficult. 

The whole masking process for one of the time series is illustrated throughout 
the panels in Figure 5.3 and consist of: 

• a) Selecting the FOV we want to analyze. Hereafter we will use the spatial- 
scale in arc seconds. 

• b) Creating, by thresholding, a binary mask for the sunspots (including umbra 
and penumbrae) and pores. 

• c) Computing the proper motions of structures by applying LCT and plotting 
the flow map for the whole FOV in granulation regions. 

• d) De-projecting the velocities according to the equations 5.1 y 5.2. The 
de-projected magnitudes are represented in gray-scale. Bright/dark areas 
correspond to large/small velocity magnitudes. The velocity units are in 

km s _1 . 

• e) Creating a binary mask for the de-projected velocity magnitudes according 
to the following rule: pixels with de-projected velocity magnitudes greater 
than 0.3 kms -1 are set to "1" and the rest to "0". 

• f) Applying the mask in e) to the flow map in panel c), so that we get 
only those projected velocities that have a de-proyected magnitude above 
the threshold of 0.3 kms -1 . Over these new flow map we identified visually 
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FIGURE 5.3 — Steps of the process applied to generate the flow map showing the large horizontal 
outflows surrounding the sunspots. a) FOV. b) Sunspot mask, c) Flow map. d) De-projected 
velocities magnitudes. The coordinates are expressed in arc seconds. 
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FIGURE 5.3 (cont.) — e) Mask of de-projected velocities with magnitudes over a certain 
threshold, f) Flow map of projected velocities with de-projected magnitudes over the threshold, 
g) Moat mask, h) Flow map of moat flows. The coordinates are expressed in arc seconds. The 
lenght of the black bar at coordinates (0,0) in panels f) and h) represents 1.5 kms -1 for the 
projected velocities. 
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the flows coming out from the sunspot that correspond to the moat flows 
described before in this work. These moat flows are easily recognized sur- 
rounding the sunspot as large-scale outflows and hence they can be outlined 
without difficulty. 

• g) Creating a binary mask of moats from the mask generated in e), with the 
help of the moat boundaries outlined from the visual inspection in f). 

• h) Plotting the final flow map showing the moat flows. We apply the moat 
mask generated in g) onto the flow map of velocities in panel c). 

5.3.3 Moat flows around the sunspots 

Our sunspots sample includes different penumbral configurations as a key factor 
to establish the moat-penumbra relation in a more robust way. Figures 5.4 to 5.8 
show within the moats, only those velocity vectors with de-projected magnitudes 
above 0.3 kms -1 after applying the whole masking procedure described above. In 
Figure 5.9 we extend this representation to the entire granulation field showing that 
large velocities are also present outside the moat. These velocities are generally 
grouped and associated with exploding granules. 

First we focus on granulation regions that display moat flows. Close inspection 
of Figures 5.4 to 5.9, and the upper panel of Figure 5.12 reveals that the velocity 
vectors in the moats are oriented following the direction of the penumbral filaments. 
We observe that in all cases the moat flow direction lines up with penumbral 
filaments that are oriented radially with respect to the sunspot center. Neverthe- 
less, the completeness of the velocity vectors (density of arrows) in the flow maps 
depends on the threshold previously imposed, as expected. This can be seen in 
Figure 5.7 where the left-side penumbra does not seem to be strictly associated 
with large flows in the granulation region. However, this penumbral part has in fact 
an associated moat flow similar to the other penumbral regions that is not visible 
in the representation since the magnitude of the velocities is slightly lower than the 
threshold of 0.3 kms -1 . 



5.3.4 Lacking moat flows 

Next we consider granulation regions close to the sunspots that do not display moat 
flows. All sunspots of our sample, except S3, have parts of the umbral core sides 
that are free from penumbrae, having direct contact with granulation regions. In 
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all these granulation regions we do not detect systematic large-scale outflows that 
fulfill our criteria for moat flows. This supports the conclusion that umbral core 
boundaries with no penumbra do not display moat flows. 

Finally, there are also granulation regions in the vicinity of penumbrae that lack 
significant moat flows. To study these cases in more detail we mark peculiar regions 
of interest in Figures 5.4, 5.6, 5.7 and 5.8 with rectangular white boxes. These 
regions are represented as close-ups in Figure 5.10. In the first three panels (T), 
the penumbral filaments display significant curvature such that they are not radially 
oriented with respect to the sunspot center but extend in a direction tangential to 
the sunspot border (as marked with black lines in the overview images Figures 5.4, 
5.6 and 5.8). The lower right panel (R) in Figure 5.10 shows a penumbra extending 
radially from the umbra. The surrounding photosphere exhibits moat flows only in 
the direction of the penumbral filaments (also see Figure 5.7). These four exam- 
ples lead to the conclusion that moat flows are not found in directions transverse 
to the penumbral filaments. 

In a few cases we found small pores in the vicinity of penumbrae in a region 
where we would expect moat flows. Such cases are seen around coordinates (19,29) 
in Figure 5.4 and coordinates (29,11) in Figure 5.8. For these regions the measured 
proper motions are less reliable since we have a FWHM=l'/0 tracking window act- 
ing on a region of only a few arc seconds (~ 2-3) between the pore and penumbra. 
Another possibility is that the pores themselves are somehow blocking the large 
outflows changing the expected behavior. In the case presented in Figure 5.8 the 
most plausible explanation for the absence of moats is the presence of a neutral 
line close to coordinates (29,11) as we shall describe in section §5.3.5. 

The findings described above suggest a link between the moat flows and the 
Evershed flows in penumbrae. We come back to this relation in section §5.4. 

The debate regarding the existence of a moat flow around umbral cores and 
individual pores is still ongoing. In a work by Deng et al. (2007), these authors 
found that the dividing line between radial inward and outward proper motions in 
the inner and outer penumbra, respectively, survived the decay phase, suggesting 
that the moat flow is still detectable after the penumbra disappeared. These issues 
will be deeper mentioned in the next chapter contrasting the results presented in 
this thesis and some more examples in the literature. 
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FIGURE 5.4 — Sunspot SI: Map of the horizontal velocities inside the moat with de-projected 
magnitudes > 0.3 km s _1 . The background corresponds to the average image of the whole series. 
The sunspot exhibits a well developed penumbra in the lower region. The moat flow is clearly 
found in the same region following the penumbra filamentary direction and in no other regions 
surrounding the sunspot. The coordinates are expressed in arc seconds. The black bar at (0,0) 
represents 1.5 kms -1 for the projected velocities. A peculiar region is shown with a white square. 
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FIGURE 5.5 — Sunspot S2: Map of the horizontal velocities inside the moat with de-projected 
magnitudes > 0.3 km s _1 . The background corresponds to the average image of the whole series. 
This sunspot shows different penumbrae shapes all around it, from short to long and from wide 
to narrow penumbrae. Except for the places where the umbra is adjacent to the surrounding 
granulation, we identified moat flows all around the sunspot in roughly the same configuration 
of the penumbra. The black bar at (0,0) represents 1.5 kms -1 for the projected velocities. The 
coordinates are expressed in arc seconds. 
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FIGURE 5.6 — Sunspot S4: Map of the horizontal velocities inside the moat with de-projected 
magnitudes > 0.3 kms -1 . The background corresponds to the average image of the whole series. 
The figure corresponds to a sunspot with two distinct penumbrae. One tangential to the umbral 
core and the other one coming out radially from the right part of the umbra. The moat flow 
continues as a prolongation of the last one. The black bar at (0,0) represents 1.5 kms -1 for the 
projected velocities. The coordinates are expressed in arc seconds. A peculiar region is shown with 
a white square. 
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FIGURE 5.7 — Sunspot S5: Map of the horizontal velocities inside the moat with de-projected 
magnitudes > 0.3 km s" 1 . The background corresponds to the average image of the whole 
series. This sunspot presents very distinct areas with and without penumbrae. The penumbral 
distributions follow radial directions coming out from the umbra. Moat flows are organized radially 
as the penumbra. The black bar at (0,0) represents 1.5 kms -1 for the projected velocities. The 
coordinates are expressed in arc seconds. A peculiar region is shown with a white square. 
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FIGURE 5.8 — Sunspot S6: Map of the horizontal velocities inside the moat with de-projected 
magnitudes > 0.3 kms -1 . The background corresponds to the average image of the whole series. 
This sunspot exhibits penumbrae on the right and upper sides. The penumbra on the right is 
mostly tangential to the sunspot border. In the upper penumbra we identified the moat flow as 
the prolongation of the penumbral filaments as expected. The black bar at (0,0) represents 1.5 
kms -1 for the projected velocities. The coordinates are expressed in arc seconds. A peculiar region 
is shown with a white square. 
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FIGURE 5.9 — Sunspot S7: Map of the horizontal velocities surrounding the sunspot with de- 
projected magnitudes > 0.3 kms -1 . The background corresponds to the average image of the 
whole series. The sunspot includes a well developed penumbra but also an empty region with no 
penumbra. Looking at the large flows, they completely disappear in the emptied areas with no 
prenumbra. The white contours outline the moat regions. The black bar at (0,0) represents 1.5 
kms - for the projected velocities. The coordinates are expressed in arc seconds 
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FIGURE 5.10 — From top to bottom and left to right: Close-ups for peculiar regions in sunspots 
SI, S4, S6 and S5 respectively. The first three regions (T) show a penumbra extending tangential 
to the sunspot border. The lower right panel (R) shows a radial penumbra. The coordinates are 
expressed in arc seconds. 
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5.3.5 Neutral lines affecting the flows behavior 

Only one of the sunspots we study displays a complete regular and well-developed 
penumbra completely surrounding the umbral core. Figure 5.12 shows the flow 
map calculated for this active region, plotting the horizontal velocities surrounding 
the sunspot with de-projected magnitudes > 0.3 km s -1 . Large outflows are not 
found in a part of the right side of the upper panel of Figure 5.12 (rectangular box). 

Following the findings of section §5.3.3, we would expect to find a moat flow 
in this region. The penumbral filaments are oriented radially from the umbral core 
which suggest the presence of a moat flow in the granulation region in the direct 
vicinity. Even with a lower velocity threshold, no moat flow can be discerned in 
this region. Nevertheless, when comparing with a simultaneous magnetogram (see 
lower panel in Figure 5.12), we found an inversion in magnetic polarity just outside 
the right border of the sunspot: the magnetogram displays positive polarity (in 
white) for the sunspot but negative polarity for the small magnetic elements and 
pore just outside the penumbra. The reversal of polarity (or neutral line) is confined 
to a narrow region that roughly coincides with the sunspot border. The absence of 
large outflows following this penumbra is suggestively related to the presence of this 
neutral line that might somehow be acting as a blocking agent for the moat flow 
in this region. The position of the neutral line measured by LOS magnetograms is 
generally influenced by the location of the sunspot (9 angle). In this case since the 
sunspot is very close to the disc center we can claim for a reliable determination of 
the neutral line. 

For sunspot in Figure 5.8 there is also a neutral line crossing all along the 
right-hand border of the sunspot from top to bottom around coordinates (20,28) 
to (28,11), respectively, as reported by De Pontieu et al. (2007) and shown here 
as Figure 5.11. The penumbral filaments extend from the umbra to the right and 
seems to bend, being forced to follow the direction of the neutral line as approach- 
ing to it. The sheared configuration of this penumbra is arranged so that the 
penumbral filaments end up along directions parallel to the sunspot border. Moat 
flows are then not found beyond this sort of penumbral configuration as mentioned 
in section §5.3.4. The same behavior when having neutral lines was presented in 
section §4.4.3 for a complex ^-configuration active region. For that active region 
we found a strong sheared neutral line crossing a penumbral border in a place 
where a moat flow was expected to follow the penumbral filaments direction but 
actually not detected. More observations of complex active regions with neutral 
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FIGURE 5.11 — Images taken in Ha wide-band (left) and Ha line-center (right) of NOAA AR 
10783 taken at the SST on 4-Oct-2005 as displayed by De Pontieu et al. (2007). The FOV shown 
in both images is 90 degrees counterclockwise rotated with respect to the one in Figure 5.8. Tick 
marks are in arc seconds. The white and black contours in the left panel outline the positive 
and negative magnetic flux from a full-disk MDI magnetogram. The white arrow indicates the 
direction of the disc center. 

lines present in the vicinity of penumbrae are needed to firmly establish the relation 
between the absence of moat flows and magnetic neutral lines. 



5.3.6 Statistics of velocity fields in moats vs. quiet granulation 

In this section the statistics of time-averaged horizontal velocity fields in moats is 
performed. For the sake of comparison, the velocity statistics in quiet granulation 
areas has also been computed. To that aim, boxes (~ 9x9 arcsec 2 ) in regions of 
less-magnetized (quiet) granulation and far from the moat flows have been man- 
ually selected in every FOV. Table 5.3 summarizes the statistical properties of the 
de-projected magnitudes of horizontal velocity fields within both moat masks and 
selected granular boxes. Here we again remark that because of the possible un- 
derestimation of the displacements by using LCT, we intend to perform not an 
absolute but a comparative statistical study between moat regions, where we find 
radially organized outflows, and quiet granulation regions far from the sunspots. 
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For each sunspot of our sample three velocity magnitude maps result from av- 
eraging over different time periods, namely: 5 min, 10 min and the whole duration 
of the corresponding time series . Table 5.3 shows that in most cases the mean 
velocity magnitude in moats is greater than in granulation. In a few cases and only 
for averages in short time periods (5, 10 min), we find similar values in both or 
even slightly lower velocities in moats. However, when averaging over a long time 
period, the difference of mean velocities in moats and granulation is in all cases 
positive and more conspicuous than for short time period averages (see column 10 
in the Table). A similar behavior is obtained for the rms parameter (see column 9 
in the Table). The maximum velocity values are also systematically larger in moats 
than in granulation. 

The described statistical behavior is expected since for short time averaging 
periods, of the order of the granulation lifetime (5-16 min, see Hirzberger et al. 
(1999) and references there in), the proper motions of the granulation structures 
compete in magnitude with the velocity of large-scale flows. However, averaging 
velocity components over long time periods results in local velocity cancellations 
in short-lived structures while steady motions at large spatial scales prevail. So 
moats can be considered as a long-term and large-scale phenomenon where the 
mean velocity exceeds that of quiet granulation by about 30%. 

The significant fluctuations in A(rms) and A(mean) (columns 9 and 10 in 
Table 5.3) for the velocity fields averaged over the whole time series, deserve a 
comment and possibly future improved measurements. They could be ascribed to 
morphological differences in the various cases considered or they could be related to 
the particular evolutionary state or magnetic field strength in the various sunspots 
of our sample. But we cannot forget that the de-projection formulas (5.1) and 
(5.2) include products, divisions and sums of products of trigonometric functions. 
Accordingly, the accuracy in the measurement of 9 and (j)' is crucial for the deter- 
mination of v . Thus, two sources of inaccuracy are: a) the assumption of constant 
heliocentric angle 9 all along the sunspot area (this has important impact in re- 
gions far from the solar disc center), and b) the determination in our FOV of the 
direction pointing to the solar disc center which defines the X-axis of the Observing 

1 Actually, when averaging over 5 min (or 10 min) we obtain more than one map since we then 
compute maps every 5 min (or 10 min) up to complete the total duration of the series (i.e when 
having 40 min of duration for a time series, we would finally obtain 8 maps of 5 min or 4 maps of 
10 min). This enable us to improve the statistics of the velocity magnitudes and reduce the noise 
in our calculations. 
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Coordinate System, i.e. the angular origin for <f>' . 

Figure 5.13 shows the histograms of de-projected velocity magnitudes for av- 
erages over the whole time series for sunspots SI to S7. Red lines correspond 
to velocities inside the moats and the blue ones to velocities in quiet granulation 
boxes. In all cases we find the same general trend in both histograms. The red 
histograms are globally shifted toward the right with respect to the blue ones. 

The left wings of the histograms for granulation lie systematically above those 
of the histograms for moats, and for SI, S2, S3, S4 and S6, both histograms nearly 
coincide at the farthest left end. This means that within the moats standard veloc- 
ity values for quiet granulation are still present although with a lower weight. Such 
is the case for the fragments of exploding granules that additionally are swept by 
large-scale moat flows. 

At some point from left to right, both histograms cross each other so that the 
right wing of the moat histogram surpasses the corresponding quiet granulation 
wing and extends to larger velocity values (> 0.6 kms -1 ). This confirms the 
predominance of large velocities in moats. The intersection point of both histograms 
corresponds to ~ 0.3 km s _1 in almost all cases (a bit higher in SI, S2 and S4). 
Interestingly, this value coincides with the threshold empirically selected a priori 
(section §5.3.2) to clearly distinguish the frontiers of the moats, which confirms 
the goodness of the limit-value (0.3 km s" 1 ) chosen to detect moats in our data 
set. 
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FIGURE 5.12 — Sunspot S3: Map of the horizontal velocities with de-projected magnitudes > 
0.3 kins -1 (upper panel). The sunspot in the figure is surrounded all around by a well-developed 
penumbra. Large outflows are not found in a certain part of the right-hand side of the spot. When 
looking at the magnetogram for the same FOV (lower panel) we found a neutral line which is 
suspiciously acting as a blocking agent for the flow coming out from the spot. The black box is 
a common region in both figures. The coordinates are expressed in arc seconds. The black bar at 
(0,0) in the flow map represents 1.5 kms -1 for the projected velocities. 
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TABLE 5.3 — Statistics of de-projected horizontal velocity magnitudes [m s 1 ] 



MOAT GRANULATION 



DURATION A(rms)t A(mean)t 



SUNSPOT 


(mi n utes) 


Max 


rms 


Mean 


Max 


rms 


Mean 


(%) 


(%) 


SI 


5 


1048 


174 


333 


947 


157 


315 


10.8 


5.7 




10 


947 


162 


315 


842 


157 


314 


3.2 


0.3 




53 


639 


137 


301 


596 


123 


250 


11.4 


20.4 


S2 


5 


1187 


183 


329 


898 


162 


307 


13.0 


7.1 




10 


1278 


182 


327 


862 


155 


300 


17.4 


9.0 




47 


921 


150 


303 


678 


115 


249 


30.4 


21.7 


S3 


5 


1036 


146 


290 


870 


147 


282 


-0.7 


2.8 




10 


868 


146 


281 


671 


137 


289 


6.6 


-2.8 




47 


689 


126 


265 


451 


101 


220 


24.8 


20.5 


S4 


5 


1060 


198 


411 


863 


164 


362 


20.7 


13.5 




10 


954 


180 


379 


780 


154 


348 


16.9 


8.9 




45 


834 


179 


365 


655 


118 


262 


51.7 


39.3 


S5 


5 


999 


178 


355 


968 


172 


307 


3.5 


15.6 




10 


972 


187 


372 


1017 


159 


291 


17.6 


27.8 




40 


802 


141 


355 


545 


107 


220 


31.8 


61.4 


S6 


5 


699 


129 


257 


683 


130 


278 


-0.8 


-7.6 




10 


673 


131 


248 


723 


136 


277 


-3.7 


-10.5 




40 


622 


125 


243 


526 


104 


210 


20.2 


15.7 


S7 


5 


872 


158 


300 


759 


142 


226 


11.3 


32.7 




10 


785 


159 


320 


828 


149 


227 


6.7 


41.0 




40 


692 


131 


285 


593 


121 


196 


8.3 


45.4 



f A(rms) and A(mean) stand for increments of the rms and mean values in moats with respect 
to quiet granulation. 




FIGURE 5.13 — Histograms of the velocity magnitudes averaging over the whole time series 
(more than 40 minutes in every case) for sunspots SI to S7. The values for moats (red line) and 
quiet granulation (blue line) are plotted. The vertical and horizontal axes represent the percentages 
and the velocity magnitudes [kms -1 ], respectively. 
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5.4 Conclusions and discussion 

Seven time series of sunspots were restored from instrumental and atmospheric 
aberrations using MFBD and MOMFBD techniques. The great quality of the 
images allowed us to study proper motions of granules outside the sunspots and 
measure their time-averaged velocities which describe the general trends of motions 
(not the same as the instantaneous velocities). 

The results from chapter §4 are extended in this chapter to a larger sample 
of active regions and systematically confirmed the same findings: a) moat flows 
are oriented following the direction of the penumbral filaments; b) in granulation 
regions found adjacent to an irregular penumbral side parallel to the penumbral fil- 
aments, moats are absent, or in other words, moats do not develop in the direction 
transverse to the penumbral filaments. Note that if the moat flows were originated 
by the blockage of the heat flux from below by the penumbra, one would expect 
moat flows directed along, but also, transverse to the direction of the penum- 
bral filaments; c) umbral core sides with no penumbra do not display moat flows. 
Moreover we include in our sample a case in which a neutral line extends along 
a penumbral border where we would expect a moat flow continuation. For this 
sunspot we do not find any moat flow following the direction of the penumbral 
filaments after crossing the penumbral border where we see a change in magnetic 
polarity. The same result was found in section §4.4.3 in a penumbral portion of a 
complex active region crossed by a strong sheared neutral line. 

All these results indicate a likely connection between the moat flows and flows 
aligned with penumbral filaments. In a recent work, Cabrera Solana et al. (2006) 
suggest that the Evershed clouds inside penumbrae propagate to the surrounding 
moat and then become MMFs after crossing the sunspot border. The MMFs dis- 
placements trace very well-defined paths that can actually be clearly seen when 
averaging magnetograms in time. Some of these MMFs are seen to start inside 
the penumbra (Sainz Dalda & Martinez Pillet, 2005). We refer the reader to sec- 
tion S3. 5. 2 for more information about these new studies. 



In this chapter we also complement the results of section §4.5 with a statistical 
analysis describing in a quantitative way the differences between velocity fields in 
moat flows and in less-magnetized solar granulation in the nearby. In contrast to 
the case of granulation, moat flows are well-organized, steady and large-scale mo- 
tions. For averages made over more than 40 minutes, the mean velocity in moats 
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(0.3 km exceeds that of quiet granulation (~ 0.23 kms -1 ) by about 30%, 
although we obtain a considerable dispersion in the results. Also the rms of the 
velocity magnitude is greater in moats by a similar percentage. 

The histograms of velocity magnitudes in the moats are broader than those 
in granulation (see Figure 5.13). The histograms of granulation show conspic- 
uous maxima, most of them in a range from 0.2 to 0.3 kms -1 , whereas the 
histograms of moats present a flatter top. Systematically, the left wing of the 
granulation histogram lies above that corresponding to the moat. At some point 
about 0.3 - 0.4 km s , both histograms cross each other and the right wing of the 
moat histogram extends beyond that of the granulation to larger velocity values (> 
0.6 kms -1 ). 

We have studied one case of a sunspot penumbra displaying a neutral line all 
along a sector of the penumbra (Figure 5.12). This neutral line is detected at the 
penumbral boundary but also emphasized by a large opposite polarity concentra- 
tion nearby. Interestingly, this penumbral sector shows no moat flow. We interpret 
this evidence as an indication of the Evershed flow being forced to go into deeper 
sub-photospheric layers at a faster pace than what is normally though to occur in 
penumbral regions not associated to neutral lines (Westendorp Plaza et al., 1997). 
A similar case, but in a sheared neutral line of a <5-spot, was found in section §4.4.3. 

The important questions related to the flows inside and outside sunspots are 
yet being studied. In the present work we contribute with a new sample of sunspots 
observed between 2003 and 2006 in six different observing campaigns. The data sets 
have been selected on basis of the seeing quality (sharpness, homogeneity, duration) 
and the availability of suitable targets: the presence of spots with some form of 
irregular penumbra. In all of our samples we follow the evolution of the sunspots for 
more than 40 min. Although this only represents a snapshot in the evolution of the 
sunspots through all their emerging and decaying processes, our sample includes 
sunspots in different evolutionary stages and penumbral configurations. 



6 

Flow field around solar pores 

This chapter is devoted to the study of flows around solar pores. After the last 
two chapters, meant to analyze horizontal flows surrounding different types 
of sunspots in detail, this one deals specifically with pores. We are interested in 
observing and analyzing pores because, different from sunspots, they do not have 
penumbra at all, so that our main conclusions about the relation between moats and 
penumbrae can be tested in these scenarios as explained in the following sections. 

6.1 Introduction 

Though there is increasing evidence linking the moat flows and the Evershed flow 
along the penumbral filaments, the debate regarding the existence of a moat flow 
around umbral cores and individual pores is still undergoing. In a recent work, Deng 
et al. (2007), found that the dividing line between radial inward and outward proper 
motions in the inner and outer penumbra, respectively, survived the decay phase, 
suggesting that the moat flow is still detectable after the penumbra disappeared. 
Figure 6.1 shows the active region studied by these authors. The sunspot evolution 
is shown over six days during its decaying phase. Yellow corks trace the dividing 
line between inward and outwards proper motions in the penumbra. 

Previous works (Sobotka et al., 1999; Roudier, Bonet & Sobotka, 2002; Hirzberger, 
2003) have measured horizontal proper motions in and around pores and have ob- 
served some penetrating flows at the umbral boundaries and a ringlike structure of 
positive divergence ("rosettas") around the pores which is related to a continuous 
activity of exploding granules in the granulation around them. Roudier, Bonet & 
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FIGURE 6.1 — Sunspot evolution over six days during its decaying phase taken from Deng et al. 
(2007). White contour lines mark the umbral and penumbral boundaries. The yellow corks trace 
the dividing line between the inward and the outward proper motions in the penumbra and around 
the residual pore. The red arrows in panel f point to the dividing line around the pore. 
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FIGURE 6.2 — Horizontal flows in a solar pore observed on 1995 June 29 and studied by Roudier, 
Bonet & Sobotka (2002). Left: Image of the pore observed on 1995 Jun 29 at the Swedish Vacuum 
Solar Telescope. Center: Horizontal velocity field with 0'.'3 resolution and averaged over 5 min. 
Right: Histogram of horizontal velocity magnitudes in the whole FOV (thick solid line), inside 
(dotted line) and outside the pore (thin solid line). The spatial scales in the first two panels are 
expressed in arc seconds. 

Sobotka (2002) identified a very clear inflow around pores which corresponds to 
the penetration of small granules and granular fragments from the photosphere 
into the pores, pushed by granular motions originated in the divergence centres 
around them. They conclude that the motions at the periphery of the pore are 
substantially and continuously influenced by the external plasma flows deposited 
by the exploding granules. Figure 6.2 shows the pore studied by Roudier, Bonet 
ii Sobotka (2002) and the horizontal velocity field averaged over 5 minutes where 
the centers of divergence are clearly identified forming the ringlike structure around 
the pore with mean velocities ~ 0.3 km s _1 . 

Pores are interesting to analyze since, as they do not display penumbrae (Keil, 
1999), what we actually observe is the direct interaction between the umbra (with 
a strong vertical magnetic field that inhibits convection inside it) and the convec- 
tive plasma in the surrounding photosphere, without any intermediate structure in 
between. Many observed features such as bright granules moving in the border of 
a pore (Sobotka et al., 1999) show the complex exchange taking place between the 
pore and the granulation surrounding it. 

In this chapter we aim to get velocity flow fields from high-resolution time 
series showing the evolution of solar regions having pores. Thus, we shall be able 
to analyze the horizontal flows around the solar pores under study to figure out 
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from the flow maps how the plasma behaves. The data analyzed are sets taken 
at both ground-based and space telescopes. The observations were obtained with 
the Swedish 1-m Solar Telescope and the space telescope HINODE in 2007. In the 
next two sections the observations and the processing applied to the images are 
described in some detail. 

6.2 Observations and data processing: SST Data 
6.2.1 Observations 

The campaign at the Swedish 1-m Solar Telescope was carried out during September- 
October 2007 and corresponds to an important collaboration among several Euro- 
pean and Japanese institutions. This was a long campaign (24 days) with joint ob- 
servations using not only the SST but also three other solar telescopes (DOT.VTT 
and THEMIS) at the Canary Islands Observatories in La Palma and Tenerife, re- 
spectively. Moreover, and for the very first time, coordinated observations with the 
space solar telescope HINODE (Kosugi et al., 2007) were performed as part of the 
Hinode Operation Program 14. 

The data from the SST analyzed in this chapter were recorded during a particu- 
lar observing run on 30 September 2007 and correspond to the active region NOAA 
10971. The characteristics of the telescope and some important details concerning 
the observations are sketched and listed in Figure 4.1 and Table 4.1 respectively. 
The main target was a region close to the solar disc center (/i=0.98) with some 
pores of different sizes embedded in a plage region that exhibits an intense mag- 
netic activity. 

The optical setup consisted in two channels: (blue and red). The blue chan- 
nel was devoted to imaging and the red one to feed the SOUP filter. In the blue 
channel, images in G-band and Call H were acquired at a high cadence so that 
the MFBD restoration technique (see section §1.5) could be applied to correct for 
atmospheric and instrumental aberrations. Table 6.1 enumerates the main param- 
eters of the images in the blue channel. 

Along the red beam the SOUP filter was placed, as mentioned before. Sec- 
tion §4.2.1 details the description of the instrument and its operation. A beam 
splitter positioned in front of the SOUP deflected a fraction of light to obtain 
simultaneous broad-band PD image-pairs in the continuum near the selected wave- 
length (Fel 6302 A). These images made up an additional object for the MOMFBD 
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FILTER 


G-band Call H 


A [A] 

A\f w hm [A] 
Exposure time [ms] 
CCD (pixels) 
Field of view ["] 
Pixel size ["] 


4305.6 3969 
11.5 1.1 

~10 ~10 
2048 x 2048 2048 x 2048 
69 x 69 69 x 69 
0.034 0.034 



TABLE 6.1 — Parameters for the blue channel used during the 2007 campaign. 

algorithm to jointly restore both, the broad-band and the narrow-band images (see 
section §1.5). From the restored narrow-band images we computed magnetograms 
and dopplergrams. 

Once the observing run was over, only the best quality images (above a certain 
threshold) in both channels (blue and red) were preserved for further processing 
and the rest were discarded. 



6.2.2 Data processing 

Post-processing of the images included first the standard flatfielding and dark cur- 
rent subtraction, as well as the elimination of hot and dark pixels and spurious 
borders in the flatfielded images. After these steps, the images were ready to be 
restored so that we could group them in sets of restoration. 

Blue beam 

In the blue channel, G-band and Call H images were independently restored by 
employing MFBD. For each wavelength, the image sequence was grouped in sets 
of about 80 frames acquired within time intervals of 10 seconds each. Every set 
yielded one restored image. So, from ~ 52640 single exposures in G-band and as 
many in Call H (more than 10 5 images adding G-band and Call H) we obtained a 
total of 658 restored images per wavelength. 

Next steps in the data post-processing were: compensation for diurnal field ro- 
tation, rigid alignment of the images, correction for distortion and finally subsonic 
filtering. For more details see section §4.3.3 where similar procedures are described 
for the SST data of the 2005 campaign. 
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SST Data 



Date 2007 Serie 



Time 
(UT) 



Duration N. images Cadence 



FOV 



mm 



sec 



30 Sep .. 



1 

2 



08:43-09:31 
09:36-09:56 



48 
20 



286 
118 



10 
10 



64.3x65.0 
64.3x65.0 



TABLE 6.2 — Characteristics of the time series of solar pores observed on 30 September 2007 
with the SST. For both wavelengths observed (G-band and Call H) we have two corresponding 
time series as listed in the table, with common parameters. 



Due to periods of bad seeing in which the quality did not reach the desired top 
level, some of the images were discarded and we kept only the best and longest 
consecutive sequence of images. The very final product were 4 time series (2 for 
G-band, 2 for Call H) with the characteristics listed in Table 6.2. The time gap of 
about 5 min between the two consecutive series for each wavelength was a conse- 
quence of a telescope tracking interruption. That is also why the FOV is slightly 
different in both time series. Most images forming the time series show details 
near the diffraction limit of the telescope. Figure 6.3 displays one of the G-band 
restored images (in false color) with close-ups enhancing very tiny bright features 
in the intergranular lanes. 



Red beam 

In the red beam images were restored with MOMFBD. The procedure followed for 
image restoration and for obtaining magnetograms and dopplergrams is the same as 
that applied to the observing material of the 2005 campaign (see section §4.3.5). 
Figure 6.4 shows one of the co-temporal sets of images of the emerging active 
region after restorations. 

6.3 Observations and data processing: HINODE Data 

After launched on the 22 nd September 2006 HINODE has become a very useful 
tool in the solar physics research. It has played an important role since it started 
pointing to the Sun and managed to observe many solar features in such a high 
detail by avoiding the distortion effects produced by the Earth's atmosphere. The 
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FIGURE 6.3 — G-band restored image (false color) of the emerging active region on 30 September 
2007. Small colored boxes are close-ups showing the high-resolution reached after the restoration 
process that enable us to detail very tiny structures in the bright features placed in the intergranular 
lanes. The typical size of the granules in the image is ~ 2". 
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public archive of HINODE is an organized data base where all HINODE observa- 
tions can be searched and downloaded in such an easy way. 

We were interested in HINODE observations of solar pores taken with the Solar 
Optical Telescope (SOT), in order to pursue our study of photospheric flows as 
explained in section §6.1. Next two sections describe the HINODE data we analyze 
in the present work. These data extend the sample of cases under study in the 
present work and moreover give us the possibility to compare with the results 
stemming from ground-based observations. 

6.3.1 Data on 1 June 2007 

Solar portion including an isolated and round-shaped pore observed by HINODE 
on 1 June 2007. Images were acquired in G-band with a cadence of 30 seconds 
covering a FOV of 27'.'85 x 55'.'70 close to the disk center (/i=0.87). After subsonic 
filtering we obtained two time series: one from 21:35 to 21:55 UT with 40 images 
and a second one from 22:26 to 23:33 UT with 134 images. Table 6.3 (1 Jun) 
summarizes the parameters of the time series in more detail. 

The treatment performed with HINODE data does not include any restora- 
tion process because of their excellent quality as a consequence of the absence of 
atmospheric turbulence. 

6.3.2 Data on 30 September 2007 

The data from HINODE on 30 September 2007 corresponds to the coordinated ob- 
servations described in section §6.2.1. HINODE observed in G-band the emerging 
flux region during almost the whole day from 00:14 to 17:59 UT with a few brief 
interruptions for calibrations. Table 6.3 (30 Sep) summarizes the parameters of the 
time series in more detail. Figure 6.5 (left) sketches the correspondence between 
the FOV observed by HINODE (51'/1 x 10272) and the one framed by the SST 
(69"6 x 69! 6). A large region of the SST FOV is covered by HINODE observations 
as shown in the figure. 

HINODE data presents quite a few misalignments due to tiny tracking flaws 
and also to every temporal interruption. Hence, we decided to align all the 1030 
images (18 hours) at a sub-pixel level. Taking the first image as the reference, next 
consecutive images were aligned to it by using a routine that made a first alignment 



See the website: http://solar-b.nao.ac.jp/hsc_e/darts_e.shtml 
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FIGURE 6.5 — Left: Sketch showing the correspondence between SST and HINODE observed 
FOVs on 30 September 2007 during a coordinated campaign. Black arrows denote the two or- 
thogonal coordinate systems and 6 stands for the rotation angle between them that amount to 
~ 65 degrees. Right: First image of the series taken with HINODE on 30 September 2007 at 
00:14:55 UT in G-band. The original FOV and the final one after data processing are shown for 
comparison. The spatial coordinates are in pixels. 



approach (in case of large misalignments). In a second iteration, the alignment with 
the first image was refined by correlating pairs of subsequent images and accumu- 
lating their misalignments. Once all images were perfectly aligned, the FOV was 
clipped to get rid of spurious borders resulting from the alignment process. 



Next step on the processing is the subsonic filtering of the series. This task 
resulted a bit tedious not only because of the short temporal interruptions (a total 
of 7 during the 18 hours) in the series but also because of the big size of the images 
and the large amount of them. To overcome this drawback we fragmented the total 
18 hours series in 8 data cubes of about 130 images each so that consecutive images 
in every cube were equally spaced. Subsonic filtering was applied independently to 
every cube. The final FOV (57'/l x 102'/2) is shown in Figure 6.5 {right panel) for 
the first image at 00:14:55 UT. 
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HINODE Data 



Date 2007 


Serie 


Time 

(UT) 


Duration 
[min] 


N. images 


Cadence 
[sec] 


FOV 

["] 


1 Jun ... 


1 

2 


21:35-21:55 
22:26-23:33 


20 
67 


40 
134 


30 
30 


27.9x55.7 
64.26x65.0 


30 Sep ... 


1 


00:14-17:59 


960 


1030 


60 


55.7x111.4 



TABLE 6.3 — Characteristics of the time series of solar pores observed on 1 June and 30 Septem- 
ber 2007 with HINODE. 



6.4 Data analysis and results 

6.4.1 General description of horizontal proper motions in the FOV 

The G-band series taken from SST and HINODE have been used to analyze the 
horizontal proper motions of structures in every FOV. Proper motions were mea- 
sured by using the same LCT-technique employed in previous chapters. In this 
section we show the maps of horizontal velocities calculated for the different time 
series using a Gaussian tracking window of FWHM 1.0" which is roughly half of 
the typical granular size. 

Figure 6.6 shows the flow maps computed from the two restored time series in 
the SST. The velocities were averaged over the total duration of each series. The 
underlying background in the maps is the average image of the respective series. As 
commented above, the FOV is slightly different in both cases so that, for instance, 
pore at coordinates (17,26) in the upper map of the figure is located at (24,42) in 
the lower map. 

The maps are dominated by flows coming from exploding granules taking place 
all over the FOV. The map averaging over a longer time period (48 min) is slightly 
smoother than the other one (20 min), and displays rather lower velocities. Never- 
theless, both maps (48 and 20 min average) reproduce in general the same proper 
motions all over the FOV. The top of the FOV in Figure 6.6 (upper panel) shows 
large flows caused by exploding granular events. These events are grouped at every 
upper corner of the FOV forming two structures that fit well the supergranular scale 
(see also the magnetogram in Figure 6.4). A smaller portion of these structures 
can also be identified in the lower panel of Figure 6.6. Even though a complete 
description of the proper motions around the pores in the FOV will be done in 
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section §6.4.5, a glance at the figure reveals no evidence of a moat-like pattern 
around any of the pores. The central part of the FOV where the smaller pores 
are embedded, exhibits a lower magnitude of horizontal velocities. This behavior is 
explained by the intense magnetic activity in this part of the FOV as unveiled by 
the correspondent magnetogram of the zone shown in Figure 6.4. 

Concerning HINODE data, the computed map of horizontal velocities for the 
time series on 1 June 2007 is displayed in Figure 6.7. The flow map shown in this 
figure is calculated for the longest time series (67 minutes, see Table 6.3). The 
isolated solar pore is immerse in a granular region displaying several and recurrent 
large exploding granular events all around it. The connection of the centers of 
the exploding events outlines a round-shaped contour at a distance from the pore 
border similar to the pore radius. The periphery of the pore is hereby influenced by 
plasma flows being deposited by exploding granular events. Velocity magnitudes in 
the outer part of the outlined contour are clearly larger than those in the inner part 
between the contour and the pore border. 

6.4.2 Averaging horizontal flows in different time intervals 

Here we study the influence of different time averages of velocity maps in a time 
series of solar pores. To that aim we consider averaging intervals of 5, 10, 15 and 
20 min in the best quality G-band time series recorded at the SST (the second 
series in Table 6.2). Figures 6.8 and 6.9 show the resulting mean velocity map for 
each case. As expected, the maps are smoother for longer averaging periods. Note 
that regardless we average over 5 or 20 min intervals, we find essentially the same 
general trends and exploding granular patterns all over the FOV, meaning that the 
averaging periods we are employing are smaller than or about the lifetime of the 
observed structures, and that even from a short time sequence we get a reliable 
picture of the plasma dynamics in the FOV. It is also worth mentioning that we do 
not recognize any moat-like flow around the pores for such different time averages. 

The statistics of the velocity magnitudes for every time average has been com- 
puted and the resulting histograms are shown in Figure 6.10. Although the distri- 
bution of velocity magnitudes is very similar in the four cases considered, one can 
easily detect that the histograms shift to the left as long as the averaging periods 
increase. Thus, the largest velocity magnitudes range from 1.83 to 1.97 km s _1 
and the mean values from 0.48 to 0.50 km s _1 . In both cases, smaller/larger values 
correspond to longer/shorter time averages. 
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FIGURE 6.6 — Map of horizontal velocities (FWHM l'.'O) from the restored time series taken at 
the SST on September 30, 2007. The velocities are averaged over 48 minutes (upper panel) and 
20 minutes (lower panel). The white contours outline the border of solar pores. The length of the 
black bar at coordinates (0,0) corresponds to 1.6 km s" 1 . The coordinates are expressed in arc 
sec. The background represents the average image of the G-band series in every case, respectively. 
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FIGURE 6.7 — Maps of horizontal velocities (FWHM l'.'O, 67 minutes average) for the processed 
time series taken by HINODE on June 1, 2007. The length of the black bar at coordinates (0,0) 
corresponds to 2.3 km s - . The coordinates are expressed in arc sec. The background represents 
the average image of the G-band series. 
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FIGURE 6.8 — Maps of horizontal velocities (FWHM l'.'O) for the restored time series taken at 
the SST on September 30, 2007, averaged over different time intervals: 5 minutes (upper panel) 
and 10 minutes (lower panel). The white contours outline the border of solar pores. The length of 
the black bar at coordinates (0,0) corresponds to 1.6 km s . The coordinates are expressed in arc 
sec. The background represents the average image of the G-band series in every case, respectively. 
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FIGURE 6.9 — Maps of horizontal velocities (FWHM l'.'O) for the restored time series taken at 
the SST on September 30, 2007, averaged over different time intervals: 15 minutes (upper panel) 
and 20 minutes (lower panel). The white contours outline the border of solar pores. The length of 
the black bar at coordinates (0,0) corresponds to 1.6 km s _1 . The coordinates are expressed in arc 
sec. The background represents the average image of the G-band series in every case respectively. 
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FIGURE 6.10 — Histogram of horizontal velocity magnitudes for different time intervals. A local 
correlation tracking technique has been used to derive the horizontal velocities (FWHM=1'.'0). 

6.4.3 Long-term evolution of the velocity field 

In order to investigate how the evolution of the emerging region affects the velocity 
in the FOV for long periods of time, we have used the HINODE time series lasting 
for several hours. HINODE data on 30 September 2007 corresponds to 18 hours 
of almost continuous solar observation as reported in section §6.3.2. The stability 
and long duration of the series is ideal for the analysis of proper motions in the 
emerging flux region. In particular we are interested in the evolution of the maps 
of horizontal velocities averaged over every hour. We have computed 14 of this 
maps starting from the first hour of observation 00:14 UT up to the time 14:00 
UT; during this long period of time we have the region of interest in our FOV. 
Afterwards, because of failures in the telescope tracking and a stop for calibration 
purposes, the FOV changes significantly and the main pore within the region moves 
further away to the border of the FOV. 

Images are grouped then in sets of one hour and every set is processed indepen- 
dently resulting in a map of horizontal velocities averaged over the total duration 
of the set. The majority of the sets have a duration of 60 minutes except for four 
of them with 46,45,45 and 55 min respectively. These four sets will be under- 
stood hereafter as one-hour sets for description purposes although the average to 
calculate the map of velocities is performed over the actual duration of every set, 
respectively. The 14 maps of horizontal velocities are computed using a Gaussian 
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tracking window of FWHM l'/O. Figures 6.11-6.14 show all the flow maps with 
the indication of the time period they correspond to, and the background image 
represents the averaged image of every 1-hour set. 

The main pore is located in the spatial position (37,55) in the first map (Fig- 
ures 6.11) and surrounded by some smaller pores which form altogether a sort of 
vertical and elongated structure in the figure. The collection of these pores is evolv- 
ing in time and some of them start merging and disappearing. The final picture of 
the region displayed in Figure 6.14 shows the isolated main pore with only a very 
tiny companion. 

We do not identify any signal of moat-like flow around the pores in any of the 
evolutionary stages shown in the maps sequence but continuous activity caused by 
exploding granules. Centers of divergence are systematically identified, some of 
them very close to the pore border. Proper motions displaying inward components 
are more common around the pores and any outward regular and large-scale flow, as 
corresponding to a moat-flow, is found. These results are firmly validated because 
of the long duration of the sample and reinforce the previous analysis pursued for 
the only 1-hour SST data. 

6.4.4 Distribution of horizontal speeds in the FOV 

The global speed distribution in a granulation field can be easily visualized by 
means of false-color maps showing the areas including velocity magnitudes within 
certain ranges. Figures 6.15 and 6.16 show the location (white areas) of velocity 
magnitudes in three different ranges in km s" 1 , for the SST data (the first series 
in Table 6.2): 

• Low velocity magnitudes: < 0.3. 

• Medium velocity magnitudes: 0.3 - 0.8. 

• Large velocity magnitudes: > 0.8. 

Small speeds are mainly grouped in the central part of the FOV where an in- 
tense magnetic activity is detected as evidenced by the high concentration of G-band 
bright points and faculae present in this region (see Figures 6.3 and 6.4). Around 
pores, the velocity magnitudes mainly correspond to the lower range (< 0.3 km s _1 ) 
so that they are surrounded by white areas in Figure 6.15 (upper panel). 
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HINODE 30.09.2007 / 00:14-00:59 UT HINODE 30.09.2007 / 01:00-01:59 UT 




HINODE 30.09.2007 / 02:00-02:59 UT HINODE 30.09.2007 / 03:00-03:59 UT 
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FIGURE 6.11 — Horizontal flows observed in solar pores regions with HINODE data on 30 Sep 
2007. The length of the black bar at coordinates (0,0) corresponds to 2.3 km s~ . The coordinate 
are expressed in arc sec in all the maps hereafter. 
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HINODE 30.09.2007 / 04:00-04:59 UT HINODE 30.09.2007 / 05:00-05:59 UT 




FIGURE 6.12 — Horizontal flows observed in solar pores regions with HINODE data on 30 Sep 
2007 (cont). 
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HINODE 30.09.2007 / 08:00-08:59 UT HINODE 30.09.2007 / 09:00-09:59 UT 




FIGURE 6.13 — Horizontal flows observed in solar pores regions with HINODE data on 30 Sep 
2007 (cont). 
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HINODE 30.09.2007 / 12:00-12:59 UT HINODE 30.09.2007 / 13:00-13:59 UT 




FIGURE 6.14 — Horizontal flows observed in solar pores regions with HINODE data on 30 Sep 
2007 (cont). 
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The areas mapping medium velocity magnitudes are regularly spread out all over 
the FOV except in the proximity of pores (Figure 6.15 lower panel). The shape of 
empty (non white) areas around the pores almost reproduce the pores shape. 

Large velocities (> 0.8 km s~ x ) are not homogeneously distributed in the FOV 
but mainly located in the two upper corners of Figure 6.16 (upper panel) where 
the granulation is less-magnetized (see magnetogram in Figure 6.4). These zones 
show large velocity flows that might reflect the presence of supergranular cells as 
commented in section §6.4.1. 



Figure 6.16 (lower panel) shows the mask of areas (in white) corresponding to 
very low horizontal velocities (< 0.2 km s _1 ) plotted over co-spatial and temporal 
Ha (false color) image. The correspondence between the bright structures in Ha 
and the regions of close-to-zero horizontal velocities is clearly seen. In the filamen- 
tary stripe-like structures in Ha connecting the two polarities (see magnetogram 
in Figure 6.4), there is a high density of low velocities (white areas) whereas this 
density shrinks when going to the upper part of the FOV in Figure 6.16 (lower 
panel). 

6.4.5 Velocity distribution around solar pores 

From the computed velocity fields we can do detailed analysis of the velocity distri- 
bution around the solar pores. To pursue this study we will use the region observed 
with the SST were we have a collection of pores to work with. As we have two 
time series for this region (see Table 6.2), we will employ the first one covering a 
longest time period, except for one of the pores which is out of the FOV. For this 
pore we will use the second time series. Figure 6.17 shows the FOV including all 
the pores under study which are labeled for easy identification hereafter. 

Figure 6.18 illustrates the projection of the velocities into radial and transversal 
components as a convenient way to compute inward and outward motions. The 
figure plots two points in the granulation surrounding a solar pore (red dots) with 
their corresponding velocity vectors v. The pore is centered at coordinates (x c ,y c ) 
with respect to the orthogonal coordinate system X,Y placed in the lower left corner 
of the FOV. Vector r is the position vector of a given point with respect to the pore 
center. The pore center is located at its gravity-center calculated by weighting the 
position of every point inside the pore with the inverse of its respective intensity. 
Velocity vectors in every point of the granulation surrounding the pore are projected 
into radial v r and transversal v t components. 
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FIGURE 6.15 — Location of the areas (in white) where the magnitude of horizontal velocities is 
lower than 0.3 km s _1 (upper panel) and between 0.3 - 0.8 km s~ 1 (lower panel). The coordinates 
are expressed in arc seconds. 
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FIGURE 6.16 — Location of the areas (in white) where the magnitude of horizontal velocities is 
greater than 0.8 km s~ 1 (upper panel). Location of small velocities overplotted (areas in white) 
on a co-temporal Ha restored image (lower panel). The coordinates are expressed in arc seconds. 
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FIGURE 6.17 — Emerging active region observed with the SST on 30 September 2007. The 
FOV include several pores (hereafter Canary Pores) which for the sake of simplicity have been 
labeled as shown in the figure to facilitate their identification. The coordinates are expressed in 
arc seconds. 
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Figure 6.18 — Sketch 

showing the projection 
applied to the velocities 
around a solar pore centered 
at (x c ,y c ) respect to the 
orthogonal coordinate system 
X, Y. The figure shows the 
velocity vectors v for two 
points in the granulation 
region around the pore (red 
dots). Projection of veloci- 
ties is done along both radial 
and transversal directions 
(dotted lines) obtaining the 
radial v r and transversal 
v t velocity components, 
respectively. 

By using the pore named as TENERIFE in Figure 6.17 we will illustrate the 
procedure followed to calculate velocity directions, throughout the panels in Fig- 
ure 6.19. In the following itemization, we will describe the steps corresponding to 
each panel in the figure: 

• a) Selection of the FOV that includes the pore under study. 

• b) Selection of a correlation box. Since active regions in general, and pores 
in particular, exhibit their own displacement (due to differential rotation and 
intrinsic motions) while embedded in the granulation pattern, we align the 
time series respect to an area framing the pore (correlation box) so that we 
make sure we are measuring plasma motions with respect to the pore. 

• c) Horizontal shifts for alignment correction as computed by cross-correlation. 

• d) Vertical shifts for alignment correction as computed by cross-correlation. 

• e) Map of horizontal velocities as computed from LCT. 

• f) Calculation of radial directions for short distances. The "radial directions" 
that we use as the reference to project the velocities are defined as perpendic- 
ular to the pore border. These perpendicular directions are calculated from 
the gradients of intensity in a smoothed pore-mask image. This way, we can 
also manage with non round-shaped pores. The limit for short distances is 
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defined by thresholding the intensity gradients. The threshold depends on 
the pore shape and size. 

• g) Calculation of radial directions for long distances. At large distances we 
consider all pores as round-shaped structures (as in Figure 6.18) and the 
radial directions are defined by the position vector r of a given point with 
respect to the pore center. 

• h) Determination of inward and outward motions. According to Figure 6.18, 
inward/outward motions correspond to cos/3 negative/positive, where (5 is 
the angle formed by v and the positive radial direction (outward) at each 
point of the FOV. The value of cos (3 is mapped in gray scale in panel h), 
ranging from 1 for purely radial outward velocities (in black) to -1 for purely 
inward velocities (in white). 

• i) Creation of a binary mask. A binary mask is created from the gray-scaled 
map in h). Areas in black and white correspond to velocities with positive 
(outward) and negative (inward) radial components, respectively. 

• j) Selection of a final clipped FOV. 

Once the method to describe the velocity directions has been set, we apply it to 
our sample of 7 different pores shown in Figure 6.17. The results are displayed in 
Figures 6.20 to 6.26 where we have: the averaged FOV around the pore, the veloc- 
ity field, the gray-scale representation of cos 0, and the final mask of velocities with 
inward (white) and outward (black) radial components. The analysis of all cases 
establish that the flows display a clear preference for inward directions around the 
pores. This fact is systematically found in all examples. The more regular-shaped 
pores (Figures 6.20 and 6.21) are surrounded by an also regular annular-like area 
with inward velocity components, having a mean width similar to the center-to- 
border distance in the pore. The dividing line between inward and outward motions 
is connecting the centers of divergence. 

Figures 6.27 and 6.28 show for the two regular-shaped pores, an average image 
(48 min) of the time series (upper panels) and a color-scaled image representing 
the horizontal velocity magnitudes (middle panels). Both pores exhibit a sort of 
proto-penumbral structure though not very emphasized. Signals of proto-penumbra 
appearance are also visible in other pores of our sample (e.g. see Figure 6.22). 
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FIGURE 6.19 — Steps of the process applied to generate a map describing the direction of the 
velocities in the vicinity of solar pores, a) FOV. b) Correlation box. c) Shifts to correct for vertical 
misalignment d) Shifts to correct for horizontal misalignment, e) Flow map. f) Perpendicular 
directions to the pore border; arrows point to the inward (negative) directions. 
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FIGURE 6.19 (cont.) — g) Radial directions; arrows point to the inward (negative) directions 
and the (white dot) represents the pore gravity center, h) Map of cos j3. i) Mask of inward 
(white) and outward (black) radial velocities, j) Map of velocities masked with i). 
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FIGURE 6.20 — Discrimination between inward and outward motions surrounding GRAN CA- 
NARIA pore. Black and white areas in the last panel correspond to velocities displaying outward 
and inward radial components, respectively. The length of the black bar at coordinates (0,0) in 
the upper right panel corresponds to 1.6 km s _1 . The spatial units are in arc sec. 
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SST 30.09.2007 - PROPER MOTIONS ANALYSIS AROUND TENERIFE PORE 




FIGURE 6.21 — Discrimination between inward and outward motions surrounding TENERIFE 
pore. Black and white areas in the last panel correspond to velocities displaying outward and 
inward radial components, respectively. The length of the black bar at coordinates (0,0) in the 
upper right panel corresponds to 1.6 km s _1 . The spatial units are in arc sec. 
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SST 30.09.2007 - PROPER MOTIONS ANALYSIS AROUND LA PALMA PORE 




FIGURE 6.22 — Discrimination between inward and outward motions surrounding LA PALMA 
pore. Black and white areas in the last panel correspond to velocities displaying outward and 
inward radial components, respectively. The length of the black bar at coordinates (0,0) in the 
upper right panel corresponds to 1.6 km s _1 . The spatial units are in arc sec. 
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SST 30.09.2007 - PROPER MOTIONS ANALYSIS AROUND LA GOMERA PORE 
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FIGURE 6.23 — Discrimination between inward and outward motions surrounding LA GOMERA 
pore. Black and white areas in the last panel correspond to velocities displaying outward and 
inward radial components, respectively. The length of the black bar at coordinates (0,0) in the 
upper right panel corresponds to 1.6 km s _1 . The spatial units are in arc sec. 
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SST 30.09.2007 - PROPER MOTIONS ANALYSIS AROUND EL HIERRO PORE 
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FIGURE 6.24 — Discrimination between inward and outward motions surrounding EL HIERRO 
pore. Black and white areas in the last panel correspond to velocities displaying outward and 
inward radial components, respectively. The length of the black bar at coordinates (0,0) in the 
upper right panel corresponds to 1.6 km s~ . The spatial units are in arc sec. 
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SST 30.09.2007 - PROPER MOTIONS ANALYSIS AROUND LANZAROTE PORE 




0246S 2 4 6 S 10 



FIGURE 6.25 — Discrimination between inward and outward motions surrounding LANZAROTE 
pore. Black and white areas in the last panel correspond to velocities displaying outward and inward 
radial components, respectively. The length of the black bar at coordinates (0,0) in the upper 
right panel corresponds to 1.6 km s _1 . The spatial units are in arc sec. 
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SST 30.09.2007 - PROPER MOTIONS ANALYSIS AROUND FUERTEVENTURA PORE 
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The FOV in Figures 6.27 and 6.28 is dominated by flows with velocity mag- 
nitudes up to 1 km s -1 . The centers of divergence are also visible and clearly 
identified in these images as black structures around the pores. Another important 
distinctive feature that comes out from these images is that the highest magnitudes 
of horizontal velocities are placed beyond these black structures away from the pore. 

In order to study how speeds of proper motions around the pore are changing, 
we calculate the mean velocity magnitude in a sort of ribbons (stripes) around 
the pore. A total of 12 adjacent ribbons (~ 0'.'24 wide), referred to as ring-like 
structures surrounding the pore, are used for this analysis. Figures 6.27 and 6.28 
(lower panels) plot the relation between speeds of proper motions versus distance to 
the pore border. The velocity magnitudes increment as a function of this distance: 
1) the mean value of speeds increments rapidly at distances ranging from ~ / /15 
to 0'.'4; 2) in the range Ql'A - l'.'O the variation curve is flatter; 3) in the range 
l'.'O - l'/6 we again obtain very sharp increments. Up to this distance from the 
pore border, we observe the same behavior in both pores. Nevertheless, the mean 
speed values at 1'.'6 are quite different (250 vs. 390 m s _1 ). At further distances 
(> 1'.'6) the trend in the mean velocity magnitudes differs substantially. We must 
have in mind that the flows around the pore can also be affected by the intrinsic 
characteristics of every single pore and by the contribution from other sources in 
the neighbourhood, e.g. pores in the vicinity. 

6.5 Conclusions and discussion 

We obtained time series for solar regions displaying pores from both ground-based 
and space observations. We then analyzed the horizontal displacements of the 
plasma around solar pores by means of LCT. This is the first time we have tested 
with HINODE data, our standard reduction routines for LCT and p-modes filtering. 
The long duration, stability and high-resolution of the data achieved by HINODE 
enable us to get extremely useful time series to study dynamical properties of the 
horizontal flows in the photosphere in periods much longer than those we are used 
to have in ground-based observations which are restricted by the changing seeing 
conditions. 

The LCT technique applied to the series allowed us to track the proper motions 
of structures in the regions nearby solar pores. We conclude that the flows calcu- 
lated from different solar observations are coherent among each other and show the 
determinant and overall influence of exploding events in the granulation around the 
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SST 30.09.2007 - VELOCITY MAGNITUDES 
OF PROPER MOTIONS AROUND GRAN CANARIA PORE 
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FIGURE 6.27 — Analysis of the velocity magnitudes of proper motions around GRAN CANARIA 
pore. Upper panel: Averaged FOV in G-band. Central panel: Magnitude of horizontal velocities in 
false-color representation. Lower panel: Plot of the mean velocity magnitude versus the distance 
from the pore border. The spatial units are in arc sec. 
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SST 30.09.2007 - VELOCITY MAGNITUDES 
OF PROPER MOTIONS AROUND TENERIFE PORE 
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FIGURE 6.28 — Analysis of the velocity magnitudes of proper motions around TENERIFE pore. 
Upper panel: Averaged FOV in G-band. Central panel: Magnitude of horizontal velocities in 
false-color representation. Lower panel: Plot of the mean velocity magnitude versus the distance 
from the pore border. The spatial units are in arc sec. 
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pores and in the whole FOV for every single set. Motions towards the pores in their 
nearest vicinity are the dominant characteristic we claim to observe systematically. 

The motions at the periphery of the pores are basically influenced by the exter- 
nal plasma flows deposited by the exploding events, as suggested by other authors 
in previous works (Sobotka et al., 1999; Roudier, Bonet & Sobotka, 2002). 

We interpret the dividing line between radial inward and outward motions, found 
by Deng et al. (2007) outside the residual pore (see panel f in Figure 6.1), as 
corresponding to the location of the centres of divergence of the exploding events 
around the pore. The outward motions these authors describe, which are not in the 
immediate surroundings of the pore but separated by the annular inward motion, 
would then correspond not to moat flows but to the outward flows originated in 
the regular mesh of divergence centers around the pore. 
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Conclusions 



The first part of this work deals with the problem of image degradation caused 
by atmospheric and instrumental aberrations. By using image restoration tech- 
niques we can correct the images and bring them near the diffraction limit of the 
telescope. One of the aims at working with these novel techniques is the correction 
of the instrumental aberrations that will be present in the SUNRISE mission, in par- 
ticular in the images produced by the IMaX instrument, a magnetograph designed 
and built by a Spanish consortium leaded by the Institute de Astrofi'sica de Canarias. 

Nowadays, the use of numerical restoration techniques is determinant and com- 
pulsory to obtain high-resolution images, and stands as a key tool for recent and 
future solar physics advances and developments. In that sense, the present work rep- 
resents an effort to have high-quality observations in IMaX by defining the method 
to calibrate the in-flight instrumental aberrations. The method is based on the 
Phase Diversity (PD) inversion algorithm. 

The robustness of the method has been tested by numerical experiments simu- 
lating different aberration components applied to synthetic solar images. Degrada- 
tion followed by PD-inversion has allowed to compare the goodness of the match 
between the input values for degradation and those resulting from the inversion. 
The sources of aberration have been modelled and added in subsequent experi- 
ments. The repercussion of every new added ingredient in the final result from 
the inversion has been evaluated. The PD-code does not accurately reproduce the 
shape of the wavefront errors but provides suitable OTFs for subsequent satisfac- 
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tory restorations. The results of the inversions can be considered reliable enough 
as to validate the method proposed for the calibration of errors in the images of 
IMaX, at least for aberrations lying in the range: rms-ripple < 2/60 and rms-LOA 
< 1/4 waves. 

The second part of the work is devoted to the observation and analysis of so- 
lar active regions. In particular we have studied the dynamical behaviour of the 
plasma around sunspots and solar pores at a photospheric level. We analyzed data 
from ground-based (SST) and space (HINODE) telescopes. Ground-based data 
are firstly corrected for atmospheric and instrumental aberrations with the methods 
described in the first part of the work, generating in this way time series of images 
displaying the solar evolution in active regions. 

Local correlation tracking techniques have been employed to infer proper mo- 
tion of structures within the FOV. We are mainly interested in the photospheric 
horizontal flows around sunspots and pores. A sample including various sunspots 
with different penumbral configurations, and solar pores was analyzed. 

The main findings concerning the existence of moats surrounding sunspots de- 
pending on the penumbral configurations exhibit by them are: 

• Moat flows are oriented following the direction of the penumbral filaments. 

• In granulation regions found adjacent to penumbral sides parallel to the 
penumbral filaments, moats are absent, or in other words, moats do not 
develop in the direction transverse to the penumbral filaments. 

• Umbral core sides with no penumbra do not display moat flows. 

• Around solar pores there are not signs of moat flows. In all the cases studied 
we have found evidence of exploding granules activity in the pore's surround- 
ings. 

Moreover we include in our study two examples in which a neutral line extends 
along a penumbral border where we would expect a moat flow continuation. For 
these cases we do not find any moat flow following the direction of the penumbral 
filaments after crossing the penumbral border where we see a change in magnetic 
polarity. The presence of neutral lines appears to be significantly affecting the be- 
haviour of the plasma flows, preventing the appearance of moats. 
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7.1 Future prospects 

More upcoming data from the HINODE telescope as well as from new facilities 
and missions like SUNRISE will provide long-time series of active regions with a 
constant image quality and enough spatial-resolution to establish firm confirmation 
of the evidences and findings presented in this thesis. The addition of simultane- 
ous Doppler and magnetogram data to the continuum intensity (or G-band) data 
sets will enhance our understanding of the link between the Evershed flow inside 
penumbrae and the moat flows around sunspots. Needless to say, the study of 
those stages were penumbrae are just formed or destroyed becomes of extreme im- 
portance to validate our findings. Similarly, the results that can be expected in the 
coming years from local helioseismology, in particular those related to the f-mode 
(Gizon et al., 2000), describing the flow patterns in the deeper layers near sunspots 
will be crucial for the establishment of a clear link between these two well-known 
flow patterns, Evershed and moat, that has not been appreciated yet. 

The interpretation of these results in terms of convective cells surrounding 
sunspots (e.g., Bovelet & Wiehr, 2003) may be different if the Evershed flow turns 
out to be the major process that injects mass into the moats surrounding sunspots. 

More observations of complex active regions with neutral lines present in the 
vicinity of penumbrae are needed to firmly establish the relation between the ab- 
sence of moat flows and magnetic neutral lines. 

We believe that one of the more relevant results put forward in this work sug- 
gesting that the moat flow is the continuation outside the spots of the Evershed 
flow, deserves adequate attention in the future from both an observational and a 
theoretical point of view. 



Conclusiones 



En la primera parte de este trabajo de tesis se aborda el problema de degradacion 
de imagen causada por aberraciones atmosfericas e instrumentales. Gracias a 
las nuevas tecnicas de reconstruccion de imageries podemos corregir estos efectos 
perturbadores y alcanzar un alto grado de resolucion cercano al h'mite de resolucion 
del telescopic 

Uno de los principales objetivos que se presiguen en esta tesis al trabajar con 
estas innovadoras tecnicas es el de aplicarlas en la correccion de las aberraciones in- 
strumentales en IMaX. Este instrumento es un magnetegrafo solar que forma parte 
de la mision SUNRISE (un telescopio solar de un metro de apertura que volara en 
un globo aerostatico y recogera imagenes solares continuamente durante 15 dias). 
IMaX esta desarrollado enteramente por un consorcio de instituciones espanolas, 
siendo el Institute de AstrofTsica de Canarias la institucion que lidera el proyecto. 

Actualmente las tecnicas de restauracion de imagenes se han convertido en her- 
ramientas indispensables para el avance en el estudio de la Fi'sica Solar. En este 
sentido, el presente trabajo representa un esfuerzo para tener observaciones de alta 
calidad en IMaX definiendo el metodo para calibrar las aberraciones instrumentales 
en vuelo. Este metodo se basa en el algoritmo de inversion de Diversidad de Fase 
(DF). 

Se ha validado la robustez del metodo de calibracion de aberraciones mediante 
experimentos numericos, simulando diferentes componentes de aberracion sobre im- 
agenes solares sinteticas. La degradacion de las imagenes seguida de su inversion 
con DF, ha sido el proceso utilizado para comparar la bondad del ajuste entre los 
valores de entrada al degradar y los obtenidos despues de la inversion. Las diver- 
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sas fuentes de aberracion se han modelado y anadido en sucesivos experimentos. 
Hemos evaluado la repercusion que tiene cada ingrediente anadido, sobre el resul- 
tado final de la inversion. El codigo de DF no reproduce en detalle la forma de 
los errores del frente de onda sino que proporciona OTFs adecuadas para producir 
restauraciones satisfactorias. Los resultados de las inversiones pueden considerarse 
lo suficientemente fiables como para validar el metodo propuesto para la calibracion 
de errores en las imagenes de IMaX, al menos para aberraciones dentro del rango: 
rms-ripple < 2/60 and rms-LOA < 1/4 waves. 

La segunda parte de la tesis esta dirigida al estudio y analisis de regiones solares 
activas. En particular se ha estudiado el com porta miento y dinamica del plasma a 
nivel fotosferico alrededor de manchas y poros solares. Para ello se han utilizado 
observaciones hechas en Tierra (SST) y del telescopio espacial (HINODE) recien- 
temente lanzado al espacio. A traves de series temporales se han estudiado los 
movimientos propios de estructuras en la zona granular que circunda manchas y 
poros, considerando una amplia muestra de ambos. 

Los principals hallazgos que surgen de este trabajo de tesis son: 

• Los flujos de gran escala aldededor de las manchas solares (llamados flujos 
de foso) se orientan siguiendo la direccion de los filamentos penumbrales. 

• En las regiones de granulacion adjacentes a penumbras en donde los filamen- 
tos son paralelos al borde penumbral, no se encuentran signos de flujos de 
foso. 

• Las regiones de granulacion adjacentes a nucleos umbrales (es decir sin 
penumbra), no muestran indicios de presencia de flujos de foso. Es decir, 
los flujos de foso no se producen en direccion transversal a los filamentos 
penumbrales. 

• Alrededor de los poros analizados no se encuentran signos de flujos de foso. 
En todos los casos se han encontrado evidencias de flujos provenientes de 
granulos explosivos alrededor de estos poros. 

Sumado a esto, se han analizado tambien un par de ejemplos en donde hay 
presencia de Imeas neutras que se extienden longitudinalmente a lo largo de bordes 
penumbrales en donde se esperaria observar flujos de foso. Para estos casos sin em- 
bargo no se ha detectado este tipo de flujos. La presencia de Imeas neutras parece 
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afectar de forma notable el comportamiento de los flujos de plasma impidiendo la 
aparicion de flujos de foso. 

Futuros estudios 

Datos provenientes del telescopio HINODE y de nuevos instrumentos y misiones 
como SUNRISE, permitiran tener muchos mas ejemplos de series temporales esta- 
bles, de alta resolucion espacial y temporal, sobre las cuales se podra estudiar la 
dinamica de las estructuras en las zonas activas para verificar los hallazgos presen- 
tados en este trabajo. El analisis de series de magnetogramas y dopplergramas si- 
multaneas a series en banda G, enriqueceran indudablemente nuestro entendimiento 
de la relacion entre el flujo Evershed dentro de la penumbra y el flujo de foso alrede- 
dor de las manchas solares. 

Sin duda el estudio de manchas en diferentes estados evolutivos donde la penum- 
bra ha sido recientemente formada o esta en fase de disolucion, es de vital impor- 
tancia para validar los resultados y conclusiones que se exponen en el presente 
trabajo. 

De igual forma aquellos resultados aportados por nuevos estudios de heliosis- 
mologia local, en particular los relacionados con los modos-f (Gizon et al., 2000) 
tendientes a describir el patron de flujos a diferentes profundidades bajo las man- 
chas solares y sus regiones circundantes, seran cruciales para establecer una clara 
relacion entre estos dos flujos: dentro y fuera de la penumbra. La interpretacion 
los flujos de foso en terminos de celdas convectivas de gran escala alrededor de 
las manchas (ver Bovelet & Wiehr, 2003) puede cambiar radicalmente si el flujo 
Evershed resultara ser el principal proceso que inyecta masa en las regiones moat 
alrededor de las manchas solares. 

Creemos que uno de los resultados mas relevantes que se defienden en este 
trabajo y que supone que el flujo de foso es una continuacion, fuera de la penumbra, 
del flujo Evershed, merece especial atencion en el futuro tanto desde el punto de 
vista observacional como del teorico. 
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A 



Zernike polynomials 

The wavefront aberration <p(r,9) in the pupil of an optical system is often 
parametrized by using a convenient set of basis functions. Because of their 
simple analytical forms circle polynomials of Zernike Zj(p,0) (Zernike, 1934) de- 
fined in a unit circle (0 < p < 1) are often used as basis functions to expand <fi(r, 9), 
as follows: 

4>(r,9) = Y,^Z j (p,9) (A.l) 

where: 

(r, 9) are the polar coordinates centered on the pupil. 

p = r/R is the normalized radial distance (R is the radius of the pupil). 

ctj are the coefficients of the terms in the expansion. 

For the sake of simplicity, e.g. the orthogonality properties, we use the Zernike 
polynomials as modified by Noll (1976). Table A.l list the first Zernike terms up 
to 6 th . 

Figure A.l shows the first few Zernike polynomials which correspond to the 
classical Seidel's aberrations 1 . 



1 Phi li pp Ludwig von Seidel (1821-1896) was a German mathematician. In 1857, von Seidel 
decomposed the first order monochromatic aberrations into five constituent aberrations. They are 
now commonly referred to as the five Seidel Aberrations. 
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TABLE A.l — Zernike polynomials (as formulated by Noll, 1976). 



ZERNIKE MODE RADIAL AZIMUTHAL 

DEGREE FREQUENCY 
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FIGURE A.l — Classical Seidel's aberrations corresponding to the first few Zernike polynomials 
in Table A.l. In the panels labeled with two Z's the second one differs only by a simple azimuthal 
rotation. 
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B 



Glossary of acronyms 



T 



he list of acronyms used in the thesis are: 



INSTITUTES, ORGANISATIONS AND GROUPS 



CSIRO Commonwealth Scientific and Industrial Research Organisation 

GACE Grupo de Astronorma y Ciencias del Espacio. 

IAA Instituto de Astrofi'sica de Andaluda 

I AC Instituto de Astrofi'sica de Canarias 

INTA Instituto Nacional de Tecnica Aeroespacial 

ISP Institute of Solar Physics - Royal Swedish Academy of Sciences 

LMSAL Loockhead Martin Solar Astrophysical Laboratory 

NOAA National Oceanic and Atmospheric Administration 

SPARG Solar Physics and upper-Atmosphere Research Group 

ULL University of La Laguna 



TELESCOPES AND INSTRUMENTATION 



OBSERVATORIES 



ENO European Northern Observatory 
ORM Roque de los Muchachos Observatory 
OT Teide Observatory 
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INSTRUMENTATION 



AO Adaptive Optics 

AOM Adaptive Optics Mirror 

BS Beamsplitter 

CAM Camera optics 

CCD Charge-Coupled Device 

COLL Collimator 

CT Correlation Tracker 

CWS Correlation Tracker and Wavefront Sensor 

E Etalon 

FP Fabry-Perot 

I Iris 

IF Interference Filter 

IMaX Imaging Magnetograph Experiment 

ISLiD Imaging Stabilisation and Light Distribution 

LCVR Liquid Cristal Variable Retarders 

M1,M2,M3 Folding Mirrors 

MDI Michelson Doppler Image 

MO Microscope Objective 

RB Reflecting Ball 

RE Reference Plate 

SOT Solar Optical Telescope 

SF Spatial Filter 

SOUP Solar Optical Universal Polarimeter 

SUFI SUNRISE Filter Imager 

SUPOS Polarimetric Spectrograph 

WFE Wavefront Error 

WFS Wavefront Sensor 

XRT X-rays Telescope 



GROUND-BASED TELESCOPES 



DOT 
SST 



Dutch Open Telescope 
Solar Swedish Tower 



THEMIS Heliographic Telescope for the Study of 

the Magnetism and Instabilities on the Sun 



SPACE TELESCOPES 

HINODE Solar-B Telescope 

SOHO Solar Heliospheric Observatory 

TRACE Transition Region and Coronal Explorer 

GENERAL DEFINITIONS 

FOV Field-of-view 

IDL Interactive Data Language 

ITP International Time Program 

LCT Local Correlation Tracking 

LDB Long Duration Balloon 

MTF Modulation Transfer Function 

PD Phase Diversity 

PSF Point Spread Function 

PV Peak-to-Valley 

RMS Root Mean Square 

RSS Root Sum Square 

SNR Signal-to-Noise Ratio 



VTT 



Vacuum Tower Telescope 



SOLAR 



DEFINITIONS 



BP 
CME 



Bright Points 

Coronal Mass Ejections 



G-cont 
LB 

MMF 



G continuum 
Light Bridge 



Moving Magnetic Features 
Penumbral Grains 
Urn bra I Dots 



PG 
UD 
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Scientific contributions 




he list of scientific contributions on the topic presented in this work is given 
below: 



Conference Contributions: 

[1] S. Vargas Dommguez, L.Rouppe van der Voort, J. A. Bonet, V. Martinez 
Pillet, Y. Katsukawa Moat flow as the continuation of the Evershed Flow. Initial 
Result from Hinode, 2007. 

[2] S. Vargas Dommguez, J. A. Bonet Navarro, V. Martinez Pillet. Evidence 
of an association between the presence of penumbra and strong radial outflows in 
sunspots. Proc. Second Solar Orbiter Workshop, 2006. 

[3] D. Orozco Suarez, L.R. Bellot Rubio, S. Vargas Dommguez, J. A. Bonet 
Navarro, V. Martinez Pillet. From simulations and IMaX observations to the VIM 
design. Proc. Second Solar Orbiter Workshop, 2006. 

Refereed Publications: 

[1] S. Vargas Dommguez, J. A. Bonet, V. Martinez Pillet. Characterization 
of horizontal flows around solar pores from high-resolution time series. Astrophys- 
ical Journal (in preparation). 
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[2] R.J. Rutten, P. Siitterlin, S. Vargas Dommguez. DOT tomography of 
the solar atmosphere. VIII. Nature and visibility of pseudo Ellerman bombs. As- 
tronomy & Astrophysics (in preparation). 

[3] S. Vargas Dommguez, L.Rouppe van der Voort, J. A. Bonet, V. Martinez 
Pillet, M. Van Noort, Y. Katsukawa. Moat flow in the vicinity of sunspots for var- 
ious penumbral configurations. 2008, Astrophysical Journal, 679, 900 

[4] R. Ishikawa, S. Tsuneta, Y. Kitakoshi, Y. Katsukawa, J. A. Bonet, S. Var- 
gas Dommguez, L.H.M.Rouppe van der Voort, Y. Sakamoto and T.Ebisuzaki. 
Relationships between magnetic foot points and G-band bright structures. 2007, 
Astronomy and Astrophysics, 472, 911 

[5] S. Vargas Dommguez, J. A. Bonet, V. Martinez Pillet, Y. Katsukawa, 
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Astrophysical Journal, 660L, 165 
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muy buenos, el futbol (PlayStation !!) de los lunes con los colegas, las doradas, 

las chuletadas, los paseitos, el Teide etc, etc, espero volver y volvere. Quiero 

mencionar tambien a Juanlu y recordarlo con alegna, la misma que el irradiaba. 

Ha sido sin duda un buen trayecto recorrido en este camino desde que sah" de 
casa hace casi 8 ahos, y en el, algunas personas especiales compartieron suenos 
y luchas. A mis amigos dudlenos Juan Felipe, Hector, Pablo Mateo, Santiago, 
Manuel, Dario y Gilberto que comenzaron conmigo esta aventura, agradezco el 
saber que cuento con ellos pese a las distancias y que siempre sera un placer re- 
unirnos como en aquellas epocas de Lingwood Hall y Firs Street en Dudley Dudley. 
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Steve was also an important part in the English experience, thanks mate for all 
your kindness and support, and all the best for ASTRAL. 

Una nueva vida sin duda se presenta ante mis ojos y eso se lo debo a Laura, 
quien a lo largo de estos ultimos anos ha sido la persona que con su energi'a, com- 
prension y amor ha sabido llenar los espacios vacios que tenia y que espero lo siga 
haciendo por muchos, muchos mas. Gracias a el la y a su inigualable familia, a sus 
padres Jesus y Maria Jose, a su hermano Carlos, a los yayos Rafa, Encarna y Julia a 
los ti'os y primos, he sentido un carino excepcional dificilmente igualable. Ahora ire 
allf a reunirme con ellos y a seguir disfrutando de unas buenas paellas domingueras, 
uhhmmmm ham nam. 

A mi familia en Colombia que ha sido un apoyo incondicional en toda mi vida, 
a mi abuelita, a mis tibs, tfas, primas y a cada uno de los que se ha alegrado con 
mi victorias y sentido tristeza con mis penas; ocuparan siempre un amplio espacio 
en mi corazon. 

Y el final, lo dedico al mas importante apoyo en mi vida, mis padres y mi her- 
mano. Tendria que escribir un capitulo entero o quizas una y hasta varias tesis para 
agradecerles con palabras todo su amor y su dedicacion. Los quiero inmensamente 
y cada dia los llevo en mis pensamientos sabiendo que todo este esfuerzo es por y 
para ellos. 

A Dios, a la vida y al Universo por hacerme tan feliz y darme siempre lo que he 
querido. 



Tenerife, Estocolmo, Valencia 
Santiago Vargas Dommguez 
2008 



0.020 




0.0 0,5 1,0 1.5 

Velocfty nfiagrntuda (krn/s) 



